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Abstract. Various aspects of the Cauchy problem for the Einstein equations are 
surveyed, with the emphasis on local solutions of the evolution equations. Particular 
attention is payed to giving a clear explanation of conceptual issues which arise 
in this context. The question of producing reduced systems of equations which 
are hyperbolic is examined in detail and some new results on that subject are 
presented. Relevant background from the theory of partial differential equations is 
also explained at some length. 

1 Introduction 

One of the most striking differences between the Newtonian theory of gravity 
and its successor, general relativity, is that in the latter the gravitational field 
acquires its own dynamical properties. Its time evolution is complicated even 
in the absence of matter. This contrasts with the fact that in the Newto- 
nian theory the field vanishes when no matter is present. The field equation, 
namely the Poisson equation, together with the boundary condition that the 
field vanishes at infinity, which is an essential part of the theory, combine 
to give this result. The Einstein equations, the field equations of general rel- 
ativity, allow idealized situations which represent gravitational waves in an 
otherwise empty universe, without any material sources. This reflects the 
different mathematical nature of the equations involved in these two cases. 
The Poisson equation is elliptic while the Einstein equations are essentially 
hyperbolic in nature. The meaning of the term 'essential' in this context is 
not simple and explaining it is a major theme in the following. 

In order to understand the special theoretical difficulties connected with 
the Einstein equations and what mathematical approaches may be appro- 
priate to overcome them, it is useful to compare gravitation with electro- 
magnetism. The motion of charged matter can be described within the full 
Maxwell theory. However, there is also another possibility, which is used when 
relativistic effects are small, for instance in many situations in plasma physics. 
Here dynamical matter is coupled to the electrostatic field generated by this 
matter at any given time. In this quasi-static model the electric field follows 
the sources in a passive way while in the full theory there are propagating 
degrees of freedom. As in the case of gravitation, the elliptic equation in the 
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non-relativistic theory (the Poisson equation again) is replaced by a system 
of hyperbolic equations, the Maxwell equations. 

The fact that the Maxwell equations are so tractable is due to their lin- 
earity, a convenient feature not shared by the Einstein equations. The theory 
of linear partial differential equations in general, and of linear hyperbolic 
equations in particular, are much better developed than the corresponding 
nonlinear theories. As a side remark it may be noted that the combined 
equations describing electromagnetic fields together with their sources are 
nonlinear and that in that context serious theoretical problems, such as that 
of describing radiation damping, do appear. 

Solutions of hyperbolic equations can be uniquely determined by their 
values on a suitable initial hypersurface. The Cauchy problem is the task 
of establishing a one to one correspondence between solutions and initial 
data, and studying further properties of this correspondence. The solution 
determined by a particular initial datum may be global, i.e. defined on the 
whole space where the equations are defined, or local, i.e. only defined on a 
neighbourhood of the initial hypersurface. 'Local' and 'global' could be called 
local and global in time since in the case where a preferred time coordinate 
is present that is exactly what they mean. 

From what has been said so far we see that in studying the Einstein equa- 
tions we are faced with a system of nonlinear hyperbolic equations. Among 
nonlinear hyperbolic equations in physics, those which have been studied 
most extensively are the Euler equations, and so we may hope to get some 
insights from that direction. At the same time, it is wise to be careful not to 
treat the analogy too uncritically, since the status of the Euler and Einstein 
equations is very different. The Euler equations are phenomcnological in na- 
ture and much is understood about how they arise from models on a more 
fundamental level. The Einstein equations have been thought of as represent- 
ing fundamental physics for most of their history and the recent idea that 
they arise as a formal limiting case in string theory will require, at the very 
least, a lot more work before it can offer a solid alternative to this. To return 
to the Euler equations, one of their well-known features is the formation of 
shocks. While there is no indication of a directly analogous phenomenon for 
the Einstein equations, it does draw attention to a fundamental fact. For lin- 
ear hyperbolic equations it is in general possible to solve the Cauchy problem 
globally, i.e. to show the existence of a global solution corresponding to each 
initial datum. For nonlinear hyperbolic equations this is much more difficult 
and whether it can be done or not must be decided on a case by case basis. 
The formation of shocks in solutions of the Euler equations is an example of 
the difficulties which can occur. A general theory for general equations can 
only be hoped for in the case of the local (in time) Cauchy problem. 

For the Einstein equations we must expect to encounter the problem that 
solutions of the Cauchy problem for nonlinear hyperbolic equations do not 
exist globally. In the case of the Einstein equations this issue is clouded by the 
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fact that the distinction between local and global solutions made above does 
not apply. To define the notions of local and global we used the concept of the 
space where the equations are defined. In other words we used a background 
space. As we will see in more detail later, in the case of the Einstein equations 
there is no background space; the space-time manifold is part of the solution. 
It is better in this case to talk only about local and global properties of 
solutions and not about local and global solutions. We may then loosely use 
the words 'local' and 'global' to refer to all aspects of the Cauchy problem 
which refer to local and global properties of solutions, respectively. Solutions 
of the Einstein equations present global features such as the formation of 
black holes which are peculiar to this system and which are made possible 
by the lack of a background space. 

In this article we are not concerned with general systems of nonlinear hy- 
perbolic equations, but with a particular one, which is given to us by general 
relativity. Actually, when the coupling to matter fields is taken into account, 
we do get a variety of hyperbolic systems. Nevertheless, we might hope that 
in at least some situations of interest, such as gravitational collapse, the dy- 
namics of the gravitational field would dominate the qualitative behaviour 
and let the effects of the particular matter model fade into the background. 
In any case, it is useful to retain the distinction between the local and global 
Cauchy problems. The global problem is what we want to solve, but the local 
problem is a natural first step. Our original plan was to cover both topics, 
but along the way we discovered that the first step is already so rich that 
on grounds of time and space we have relegated global questions to passing 
remarks. Along the way we stumbled over a variety of 'well-known' things 
which turned out not to be known at all, or even to be false. 

The theory of the Cauchy problem allows us to formulate and establish 
relativistic causality within general relativity. Another basic function of the 
solution of the Cauchy problem is to parametrize solutions of the field equa- 
tions in a useful way. To single out the class of solutions relevant to the 
description of a given physical situation, we can single out an appropriate 
subclass of initial data, which is often simpler to do. This does not mean 
that by identifying this class of solutions we have solved the problem. Rather 
it means that we have found a class of problems which may be of physical 
relevance and which it is therefore desirable to investigate mathematically. In 
the end the central mathematical problem is to discover some of the global 
properties of the solutions being studied. As has already been said, this ar- 
ticle is almost entirely restricted to local questions, so that we will say little 
more about this point. However we will return to it in the last section. 

The point of view of the Cauchy problem can also be used to throw light 
on various other issues. For instance, it provides a framework in which it 
can be shown that certain approximation methods used in general relativity 
really provide approximations to solutions of the Einstein equations. It can 
help to provide an analytical basis for the development of efficient and reliable 
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numerical schemes to solve the Einstein equations. The Cauchy problem can 
furnish examples which throw light on general conjectures about solutions of 
the Einstein equations. It can be used to investigate whether certain prop- 
erties of explicit solutions of the Einstein equations of importance in general 
relativity are stable to perturbations. 

The structure of the article will now be outlined. In the second section a 
number of fundamental concepts are introduced and critically reviewed. In 
particular, the splitting of the field equations into evolution equations and 
constraints is described. The rest of the article concentrates on the evolu- 
tion equations. The question of how to solve the constraints on the initial 
hypersurface is not considered further. The second section presents the basic 
elements which go into solving the local Cauchy problem. It discusses in par- 
ticular the question of gauge freedom and how to show that the constraints 
are satisfied everywhere, given that they are satisfied on the initial hyper- 
surface (propagation of the constraints). These are aspects of the question 
of hyperbolic reduction, i.e. how questions about the Cauchy problem for the 
Einstein equations can be reduced to questions about the Cauchy problem 
for hyperbolic equations. 

The solution of the Cauchy problem relies on the use of techniques from 
the theory of partial differential equations (PDE). The third section presents 
some of the relevant techniques and attempts to explain some of the im- 
portant concepts which play a role in the theory of the Cauchy problem for 
hyperbolic equations. We have chosen to concentrate on symmetric hyper- 
bolic systems rather than on other kinds of hyperbolic equations such as 
nonlinear wave equations, which could also be used as a basis for studying 
the Cauchy problem for the Einstein equations. It will be seen that symmet- 
ric hyperbolic systems provide a very flexible tool. A comparative discussion 
of various notions of hyperbolicity is also given. 

In the fourth section we return to the question of hyperbolic reduction. 
Different ways of reducing the Einstein equations in 3+1 form are presented 
and compared. The ADM equations, which were already introduced in Sect. 
||, are discussed further. A form of these equations which has proved successful 
in numerical calculations is analysed. A number of other illustrative examples 
are treated in detail. The first example is that of the Einstein-Euler system 
for a self-gravitating perfect fluid. The case of dust, which is significantly 
different from that of a fluid with non- vanishing pressure, and which leads 
to serious problems in some approaches, is successfully handled. The second 
example is a variant of the pure Cauchy problem, namely the initial boundary 
value problem. The third example is the Einstcin-Dirac system, which gives 
rise to particular problems of its own. 

The fifth section starts with a discussion of the proofs of the basic state- 
ments of local existence, uniqueness and stability for the Einstein equations, 
based on the PDE theory reviewed in the third section and a particular hy- 
perbolic reduction already introduced in Sect. @. Then the extension of these 
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results to the case where different kinds of matter are present is sketched. 
Existence, uniqueness and stability constitute together the statement that 
the Cauchy problem is well-posed. The significance of this is underlined by 
an example of a situation where the Cauchy problem is ill-posed. Finally it 
is shown how existence and uniqueness imply the inheritance of symmetries 
of the initial data by the corresponding solutions. 

The choice of topics covered in this article is heavily dependent on the 
interests of the authors. In the last section we list a number of the important 
topics which were not discussed. While this list is also influenced by personal 
taste, we hope that it will provide the interested reader with the opportunity 
to form a balanced view of the subject. 

2 Basic Observations and Concepts 

In this section we give an introductory survey of various aspects of the field 
equations. Most of them will be discussed again, in greater detail, in later 
sections. 

Gravitational fields represented by isometric space-times must be consid- 
ered as physically equivalent and therefore the field equations for the metric 
must have the property that they determine isomctry classes of solutions, not 
specific coordinate representations. This is achieved by the covariant Einstein 
equations 

Rfj.i' — g + ^9fiv — k-T^ v . 

It is often said that many of the specific features of the Einstein equations are 
related to this covariance. One may wonder, however, what is so special about 
it, since the wave equation, the Yang-Mills equations, the Euler equations, 
just to name a few examples, are also covariant. The difference lies in the 
fundamental nature of the metric field. While the examples just quoted are 
defined with respect to some background structure, namely a given Lorentz 
space (in the case of the Yang-Mills equations in four dimensions a conformal 
structure suffices), the Einstein equations are designed to determine Lorentz 
manifolds without introducing any extraneous structures. Consequently, the 
solutions of the equations themselves provide the background on which the 
equations are to be solved. 

This makes it evident that Einstein's equations can at best be quasi-linear. 
That they are in fact quasi-linear can be seen from their explicit expression. 
Writing the equation in the form 

Rfjiv — ^Qjxv = k {Tfj,v — ^ T g^ v ), (2.1) 

the principal part of the differential operator of second order which acts on 
the metric coefficients is given by the left hand side. In arbitrary coordinates 
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we have 

R - 1 n A P f d 2 9tlu d 2 gx P d 2 g w d 2 g pv \ 

2 9 { 8x* dxP + dx» dx» dx x dx» dx»dx x \ +WiluK9, 9h 

' (2 ' 2) 

where Q denotes a rational function of the metric coefficients and their first 
order derivatives. We see that the equations are quasi-linear but not better, 
i.e. they are linear in the derivatives of highest order but these derivatives 
come with coefficients which are given by the unknown itself. The Euler 
equations are also quasi-linear in this sense. However, they are studied in 
general on an independent background space-time in terms of which their 
solutions may be analysed. 



2.1 The Principal Symbol 

Consider a system of k partial differential equations of order m for an Re- 
valued unknown u defined on an open subset U of R™. Suppose that in 
coordinates x^ , /i = 1 , . . . , n on U it takes the form 

Pu= A a D a u = f. (2.3) 

\a\<m 

Here the A a denote smooth real k x fc-matrix-valued functions, f(x) is a 
smooth Revalued function, the a = (ai, . . . , a n ), with non-negative integers 
otj, are multi-indices, |a| = ct\ + . . . + a n , and D a = d"l . . . d™n- Assume, 
first, that the equations are linear, so that A a = A a {x). 

The question of the formal solvability of the equation leads to the impor- 
tant notion of a characteristic. Suppose that H = {<£ = const., ^ 0} is a 
hypersurface of U, defined in terms of a smooth function <P. A Cauchy data 



set on H for equation (2.3) consists of a set of functions ua,Ui, . . . ,M( m _i) 
on H . Let the coordinates x^ be chosen such that <P = x 1 near H . Interpret- 
ing the functions m as the derivatives d l xl u of a solution u to our equation, 
defined in a neighbourhood of H, we ask whether equation ( |2.3| ) allows the 
function 9™ u to be determined uniquely on H from the Cauchy data. Since 
all functions D a u with ot\ < m — 1, can obviously be derived from the data, 
it follows that we can solve the equation on H for <9™ u if and only if the 
matrix ^4( m <°> - >°) i s invertible at all points of H. Using for any covector £ M 
and multi-index a the notation = S^ 1 . . . , we observe that with our 
assumptions A ( - m '°>-'°') = J2\ a \= m A a D a The invertibility of the matrix 
on the right hand side is independent of the coordinates chosen to represent 



equation (2.3) and the function <P used to represent H. This follows from 
the transformation law of the coefficients of the differential equation and of 
covectors under coordinate transformations. 

We are thus led to introduce the following covariant notions. For a cov- 



ector £ at a given point x S M, we define the principal symbol of (2.3) 



The Cauchy Problem 7 



at (x,£) as the matrix (r(x,£) = J2\a\=m A a (%) £ a ( or the associated linear 
map). The reader should be warned that a definition with replaced by 
where i is the imaginary unit, is often used, since it fits well with the Fourier 
transform. The latter convention will not be used here. A hypersurface H, 
represented by a function <P as above, is called nowhere characteristic for 
( |2.3| ), if det(o~(x, dti?)) ^ for x £ H and we say that H is a characteristic 
hypersurface, or simply a characteristic for (2.3), if det (a (x, d<P(x))) = for 
x £ H. 

Given Cauchy data on a hypersurface H which is nowhere characteristic, 
we can determine a formal expansion of a possible solution u in terms of 
x 1 on H by taking formal derivatives of (2J3) and solving for d l x± u on H, 
i = m, m + 1, .... Conversely, if H is characteristic, Cauchy data cannot 
be prescribed freely on H , because the equation implies relations among the 
Cauchy data on H. We refer to these relations as the inner equations induced 
by the equations on the characteristic. 

If the rank of the matrix a(x, d<P(x)) is k — j with some positive integer j, 
there are j such relations. The principal symbol is not sufficient to describe 
the precise nature of these relations. The complete information of (2.3) is 
needed for this. 

If equation (2.3) is quasi-linear, so that the coefficients of the equation do 
not only depend on the points of M but also on the unknown and its deriva- 
tives of order less than m, we have to proceed slightly differently. Suppose 
u is a solution of the quasi-linear equation ( [2.3| ) . For given covector £ at the 
point x £ M we define the principal symbol of ( |2 . 3| ) with respect to u as the 
matrix a(x,£,) = J2\ a \=m A °( x > u ( x )> ■ ■ ■ - D& u(%)\\!3\=m-i) and use it to 
define, by the condition above, the characteristics of ( |2.3| ) with respect to u. 
Thus for quasi-linear equations the characteristics depend on the solution. 

We return to the Einstein equations. Assuming that (M, g) is a solution of 
( [j~l| ) of dimension n, we find that for given covector £ at x £ M the principal 
symbol of the operator in (2.2) defines the linear map 

k{iv > (ct • k)fu/ — — — g p {x) {k^ £\ £ p + k\ p £ p £„ — k pp t;\ £ v — k\ v £^ £ p } , 

(2-4) 

of the set of symmetric covariant tensors at x into itself. If k pv is in its kernel, 
we have 

= k pv £ p £ p + g p ^ k p x £p — k^ p £ p — k vp £ p (2-5) 

from which we see that tensors of the form fc p „ = £ p f} v +^ v rj^, with arbitrary 
covector 77, generate an n-dimensional subspace of the kernel. This subspace 
coincides with the kernel if £ p £ p 7^ 0. If £ P £ P = 0, £ p 7^ 0, equation (2.5) 
takes the form £ M 77^ + ^uVfj, = 0, with some covector r\ p . Since ^ / we 
must have r\ p = or, equivalcntly, | g px k p \ — k pv tf . The solutions, given 
by k pv = h pv +a (£ M r)*+^ r)*) with a £ R, 77* a fixed covector with £ p 77* 7^ 0, 

and h pv satisfying £ p — 0, g pv h^v = 0, span a space of dimension , 
which is strictly larger than n if n > 4. 
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It follows that any hypersurface is characteristic for the Einstein equa- 
tions. In the case of spacelike or timelike hypersurfaces n relations are im- 
plied on Cauchy data sets, in the case of null hypersurfaces there are " ^~ 3 ^ 
additional relations. The fact that the rank of the principal symbol drops fur- 
ther for covectors £ which are null has the immediate consequence that the 
inner equations induced on null hypersurfaces arc completely different from 
the inner equations induced on nowhere characteristic hypersurfaces. This is 
related to the fact that null hypersurfaces (and appropriate generalizations 
admitting caustics) may represent wave fronts, swept out by high frequency 
perturbations of the field. The propagation of the latter is governed by the 
inner equations induced on these hypersurfaces (cf. [^7j and, for a modern 
discussion, also (88) for the mathematical aspects of this statement). 

The fact that on any hypersurface the tensors of the form k^ v = £ M r\ v + 
£,i> are m the kernel of the principal symbol map at (x, £) can be deduced by 
an abstract argument which uses only the covariance and the quasi-linearity 
of the equations (cf. fH| , ]6^| ). This emphasizes again the special role of null 
hypersurfaces. Assume that cf> T is the flow of a vector field X. If we denote 
by 4>* the associated pull-back operation, we have 

#(Ric[j]) = Ric[#(j)]. 
Taking derivatives with respect to r we obtain 

£ x Ric{g] = RiCg[£x g], 

where the right hand side denotes the (at g) linearized Ricci operator applied 
to the Lie derivative Cx g^v — V M X v + V„ X^ of the metric. Considering 
this equation as a differential equation for X, the principal part is given by 
the terms of third order appearing on the right hand side. Since the equation 
is an identity, holding for all vector fields, these terms must cancel. It follows 
that the expressions £ M X v + £ v Xa are in the kernel of the principal symbol 
of the linearized and thus also of the non-linear Ricci operator. 

For the further discussion it will be convenient to use coordinates related 
to a given hypersurface. For simplicity and since it is the case of most interest 
to us, we will assume the latter to be spacelike and denote it by S. Let T be 
a non-vanishing time flow vector field transverse to S (but not necesssarily 
timelike) and t = x° a function with {t — 0} = S and < dt,T > = 1. 
We choose local coordinates x a , a — 1,2,3, on S and extend them to a 
neighbourhood of S such that < dx a , T > = 0. In these coordinates T = d t 
and the metric has the ADM representation 

g = -{adtf + h ab ([3 a dt + dx a ) ((3 b dt + dx b ), (2.6) 

where h = h a b dx a dx b represents the induced Riemannian metric on the slice 
S. To ensure the non-degeneracy of the metric we assume the lapse function 
a to be positive. We write j3 = (3^ L <9 M = (3 a d a for the shift vector field, which 
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is tangent to the time slices. We assume further that the metric h induced on 
the time slices is Riemannian and we denote its contravariant form by h ab , 
so that h ab h bc = S a c . 

The (future-directed) unit normal to the time slices is given by = 
- o — whence n p = — aS° p . We write hf 1 „ = g p v + n u for 
the orthogonal projector onto the time slices, hu V = g plJ + n p for the 4- 
dimensional representation of the interior metric on the time slices, and use 
g to perform index shifts. The second fundamental form (extrinsic curvature) 
of the time slices is given by 

y = h x h p \7\n =-h x h P C a\ =-£ h 

\p,v — ">u v A tb p 2 P> u n 2 n ' 

For tensor fields . . . v intrinsic to the time slices, i.e. having vanishing con- 
tractions with n, the /i-covariant derivative D on a fixed time slice is given 

by 

D p t„... v = h x p h* ^..h*"?^..* . 

Finally, we write n(f) = n p d p /, note that a p = n v V„ n p = D p log(a), 
and use the coefficients 7 M v p = r\ v -n h x p h v v K* p which represent the co- 
variant derivative D in the sense that the components "f a b c are the Christof- 
fel symbols of the metric h a b- The connection coefficients of g^ can then be 
written 

f/ P = ^n" n p n(\og a) + n p a v n p - a p n v n p - n M n v a p (2.7) 

+n M 7i p — /3" w n 7 * - n p — V ,„h* p-n p — j3 v ~ „ 
a a a 

+n v Xp-p ~ n p x v P ~ Xp. v n p + 7 M v p . 



2.2 The Constraints 

To isolate the geometrically relevant information contained in the Cauchy 
data, we reduce the coordinate freedom tranverse to S by assuming T to be 
a geodesic unit vector field normal to S. The metric then takes the form 

g = -dt 2 + h ab dx a dx b , (2.8) 

and the pull-back of the second fundamental to the time slices is given by 

Xab = \ d t Kb- (2.9) 



Expressions (2.8) and (2.9) suggest that the essential Cauchy data for the 
metric field are given by the induced metric h ab and the second fundamental 
form Xab on the spacelike hypersurface S. Since coordinate transformations 
are unrestricted, the n inner relations induced on general Cauchy data on a 
spacelike hypersurface 5* must be conditions on the essential Cauchy data. 
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If we contract the covector £, which corresponds in our context to the 
normal n of S, with the principal symbol of the Einstein tensor G p „ = i? p „ — 
g^ v and evaluated on fc, we get 

e ((<J-k)^-^g^g xp (a-k) x ^ =0. 

This identity indicates the combinations of the equations which contain only 
derivatives of first order in directions transverse to 5*. Indeed, if we express 
the equations 

= Z H = n»n» (Gp, +\g fl „ - kT„ v ), 

= Z? = r/ 1 hu p (G MP + A g w - k T mp ), 

on S in terms of h and %, write p — n v T pl/ , j u = — tt/ 1 h p v T^ p , and pull- 
back to S, we get the constraint equations on spacelike hypersurfaces, i.e. the 
Hamiltonian constraint 

= 2 Z H = r - X ab X ah + (Xa a ? ~ 2 A - 2 k p, (2.10) 

and the momentum constraint 

= Z b M = D aXb a -D b X a a + Kj b , (2.11) 

where r denotes the Ricci scalar and D the connection of h. 

These equations have the following geometric meaning. On a spacclike 
hypersurface S the curvature tensors R^ vp7] and upr) of g and h respectively 
and the second fundamental form Xi*" are related, irrespective of any field 
equation, by the Gauss equation 

r" vm = h" x R X .H> h^uh^ph^rj-x^p Xur, + X M « Xu P (2.12) 
and the Codazzi equation 

D„ Xu V - D v Xm = -nx R X «w h * v h* „ h* „. (2.13) 



The constraint equations follow from (2.12) and (2.13) by contractions, the 



use of the field equations, and pull-back to S. Thus the constraints represent 
the covariant condition for the isometric embeddibility of an initial data set 
(S, h a b, Xab, P,ja) into a solution of the Einstein equations. 

We note here that the constraints ( |2.1C| ) and ( |2.1l| ), which are analogues 
of the constraints of Maxwell's equations, have important physical conse- 
quences. One of the most important of these is the positivity of the mass 
which can be associated with an asymptotically flat initial data set (subject 
to reasonable conditions) p2j , [B4j . 
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2.3 The Bianchi Identities 

Before analysing the structure of the field equations further, we note some 
important identities. The Riemann tensor Riem[(7] of the metric g 7 given by 

R<* ux P = d x r p » v - d p a " „ + r x " s r p s „ - r p » s r x s v , (2.14) 

where _T p M v denotes the Christoffel symbols of g pl , , has the covariance prop- 
erty 

Rien#*(<?)] = ^(Riemtff]), 

where (f> denotes a diffeomorphism of M into itself. Two important identities 
are a direct (cf. p4[]) consequence of this, the first Bianchi identity 

R" xu P + R" 1 P xu + R" 1 u P x = 0, (2.15) 
and the second Bianchi identity 

V M W xu P + V p Ri X pu + V„ W Xm = 0. (2.16) 
The latter implies the further identities 

V li R"v\p = V\R vp -VpRv\, (2.17) 

V"V-^„-R = 0. (2.18) 

The second Bianchi identity will serve us two quite different purposes. 
Firstly, it will allow us to resolve certain problems arising from the degeneracy 
of the principal symbol considered above (it is the integrability condition 
which allows us to show the propagation of suitably chosen gauge conditions 
and the preservation of the constraints). Secondly, it will provide us with 
alternative representations of the field equations. 

2.4 The Evolution Equations 

In this section we shall discuss a few basic ideas about the evolution problem. 
Our observations about the constraints and the decomposition of 

Z pv — p,v "T" ^Qp,v ^-^/xy? 

given by 

V = Z% - n v Zf - n„ Zf + n p n v Z H , (2.19) 

with Z^ v = h p x h p v Zx p , suggest that the basic information on the evolution 
equations should be contained in 
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or any combination of it with the constraints. To obtain simple expressions 
in terms of the field h^u and the second fundamental form x^f > 

defined by 

the generalization 

Crh^ v = 2ax»v +Cph tiv , (2.20) 



of fl2.9| ), we consider the equation 

V = Z%- 1 -h iXV {h^Zl p -Z H ) (2.21) 

= - (Ar X»v - Xp.v - Dfj. D u a) + r^y + x P P X^-2 X PP Xu p 
a 



Together with (2.20) it should be regarded as an evolution equation for the 
fields hfja,, x^v 

It then appears natural to analyse the general solution of the Einstein 
equations by the following procedure: Find initial data, i e. a solution h ab , 
Xat of the constraints, on the slice S = {t = 0}. Then find the solution h ab , 
Xab of the equations 

dth ab = 2axab + £phab, (2.22) 
d t Xab = -a {r ab + Xc C Xab~2 Xac Xb c ) (2.23) 

+D a D b a + Xab + a(Xh ab + k (T ab - - T h ab )), 

equivalent to ( [2.20 ) and ( [2.21 ), which induces these data on S. The first 
step will not be considered further in this article; we shall give some relevant 
references on the problem of solving the constraint equations in Sect. |(| Here 
we want to comment on the second step, which raises several questions. 

i) What determines the functions a, (3 a ? Is it possible to prescribe them, 
at least locally near S, as arbitrary functions a = a(t,x c ) 7 (3 a — (3 a (t,x c ) of 
the coordinates, possibly with the restriction a 2 — (3 C (3 C > 0, which would 
make dt timelike? We could give a positive answer to this question, if, starting 
from a representation of the metric (^j]) in terms of some coordinate system 

with t' = x° = on S, we could always find a coordinate transformation 
t = t(x fJ/ ), x a — x a (x >1 ) with t(0, x a ) = 0, which casts the metric into the 
desired form, i.e. achieves 

1 = 9 00 =g^'(x X ')tyty, (2.24) 



a 2 {t,x c ) 

1 



a 2 (t, x c ) 



(3 a (t, x c ) = ff ua = g» v (x x ) V x a y . (2.25) 



If the left hand side of the first equation only depended on t and x^ , the 
standard theory of first order PDE's for a single unknown could be applied 
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to this equation and second equation would essentially reduce to an ODE. 
However, in general this theory does not apply, because the dependence of 
the function a on x a introduces a coupling to the second equation. 

ii) Suppose we could prescribe lapse and shift arbitrarily. Could we then 
show the existence of a (unique) solution h a b, \ab of the initial value prob- 
lem for the equations (2.22) and (2.23) (possibly coupled to some matter 
equations)? 

iii) Suppose we could answ er th e last question positively, wou ld th e re- 
sulting solution to ( [2.22 ) and ( 2.23 ) then satisfy the constraints ( 2.10 ) and 
( [2.11 ) on the slices {t = const.}? Only then would we know that the metric 
9^v, obtained from our fields h a t,, Xab, oe P c by ( |2.6| ), is a solution of the 
Einstein equations. 

We can answer question (iii) as follows. Suppose equation ( [2.21 ) is satisfied 
on a set M =] — a, c[xS, with a, c > 0, and the solution induces the given 
data o n {0} x S, which we assume to be identified with S in the obvious way. 
Since ( 2.21 ) is equivalent to Z^ v — Z H = 0, we can write on M by ( 2.19 ) 



G>„ + A g^ u — k 



+ {2 nil n v + g^}Z H . 



Taking the divergence, using the contracted Bianchi identity ( 2.1§| ), assuming 
that the matter field equations have been given such as to ensure = 0, 

and splitting into normal and tangential parts, we get the equations 



H 



tM 



tM , 



2Z H V,n p , 



^ V p -D V Z H = -Zjf V : 



-z 



M 



-Z p n,j n 



2Z 



H 



which imply subsidiary equations, satisfied by Z H and Zff , 



i 

±h at 



o 



-Ipd _ h cd 
-h ad -±h ac 3 



z H 
z™ 



where h, h a denote linear functions of Z H , Z^f . 
v = t {Z H \Z^) of the form 

Bv = 0, 



(2.26) 

Since it is a system for 



with symmetric matrices and a positive definite matrix A , it is symmetric 
hyperbolic (cf. Sect. 3.1). Moreover, it has characteristic polynomial 



det(A" £ M ) = - det(h ab ) (n" £ P ) 2 £„, 



(2.27) 



which implies that its characteristics are hypersurfaces which are timelike or 
null with respect to the metric g^ v . 

Consequently, if S' = {(f> = 0,<j> tll ^ 0} C M, with tf> € C°°(M), is 
a spacelike hypersurface, the matrix A^ <p^ is positive definite on S' . Sup- 
pose 5", S" are two spacelike hypersurfaces which intersect at their common 
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2— dimensional boundary Z and bound a compact 'lens-shaped region' in M. 
Then it follows from the discussion of symmetric hyperbolic systems in Sect. 
|3.1| that the fields Z H , Z^ 1 must vanish on S" if they vanish on S' . 

To make a precise statement about the consequences of this property, 
we need the important notion of the domain of dependence. Let us assume 
that there is given a time orientation on (M,g). If U is a closed subset of 
M we define the future (past) domain of dependence of U in M as the set of 
points x € M such that any g— non-spacelike curve in M through x which is 
inextendible in the past (future) intersects U. We denote this set by D + (U) 
(resp. D~{U)). 

It can be shown that the result about lens-shaped regions referred to above 
and the fact that the fields Z H Zf vanish on 5" imply that Z H and Zf , 
whence also Zjf , vanish on the domain of dependence D(S) = D + (S)UD~(S) 
of S in M. This shows the p reser vation of the constraints under the evolution 
defined by equations ( 2.22j ) and ( 2.23| ) and the prescribed lapse and shift. 



Questions (i), (ii) are more delicate. Let us assume that the coefficients 
gfi v ^ are rea j analytic functions for x x in an open subset V of R 4 with 
V n {x° — 0} 7^ 0, a nd th at a > an d (3 C are real analytic functions of t and 
x a . Then equations ( [2.24 ) and ( 2.25 ) can be written in the form 



t t0 > = F 1 (t,x a ,t tC ,,x a >c >,x" ), x a , , =F 2 (t,x a ,t, c >,x a tC ,,x» ), 

with functions F\ , F2 which are real analytic for (t , x a , t )C / , x a >c i , x^ ) £ R 16 x 
V. Thus, by the theorem of Cauchy-Kowalewskaya (cf. po| ), the differential 
problem considered in question (i) can be solved in a neighbourhood of the set 
{x° = 0} C V. Using the covariance of equations ( 2.24 ), ( |2.25 ), it follows that 



given a real analytic Lorentz space (M, g), an analytic spacelike hypersurface 
S in M with coordinates x a on S, and analytic functions a = a(t,x c ) > 0, 
pa _ pa(t^ x c^ there exist unique real analytic coordinates t, x a on some 
neighbourhood of S in M such that t = on S and the lapse and shift in the 
expression of g in these coordinates are given by a and (3 a . 

Assume now that the 3— manifold S, the initial data h ao , \ab solving the 
constraints on S, as well as the functions a(t,x c ) > 0, (3 a (t,x c ) are real 
analytic and assume for simplicity that T^ v — 0. Then we can derive from 
( |2.22| ) and ( [2.23 ) a differential system for u — (h a },,k a b c = h a0iC ,Xab) of 



the form dtu = H(u,t, x a ) with a function H which is real analytic where 
det(h a b) 7^ 0. Again the theorem of Cauchy-Kowalewskaya tells us that this 
system, whence also ( [2.22 ) and ( 2.23| ), has a unique real analytic solution on 



R x S near {0} x S for the data which are given on {0} x S after identifying 
the latter in the obvious way with S. By our discussion of (iii) we know that 
we thus obtain a unique analytic solution to the full Einstein equations. 

It should be noted that the solution obtained in this way depends a priori 
not only on the data but also on the chosen lapse and shift. That the latter 
do in fact only affect the coordinate representation of the solution can be seen 
as follows. Given another set of analytic functions a'(t',x c ) > 0, (3 a (t',x c ), 
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we can either find, as remarked above, a coordinate transformation t — > t', 
x a — > x a such that the given metric has the new values of lapse and shift in 



the new coordinates , or we can deal with the initial value problem for (2.22) 



and (2.23) with the new lapse and shift. However, due to the uniqueness of 
this solution it must coincide with the first solution in its new coordinate 
representation. 

It is a remarkable fact that in the course of solving the Einstein equations 
we can prescribe rather arbitrarily four functions a, (3 a (in the analytic case) 
which are considered first as functions on some abstract R 4 but which, once 
the solution has been constructed, acquire the meaning of lapse and shift for 
the coordinate expression of the metric. Since the coordinates in which they 
have this meaning are denned by a, /3 a implicitly (via the field equations), 
we refer to these functions as the gauge source functions of our procedure 
(we shall see below that, depending on the chosen equations, quite different 
objects can play the role of gauge source functions) and to the act of pre- 
scribing these functions as imposing a gauge condition. The considerations 
above show also that the manifold on which the solution is constructed must 
be regarded as part of the solution. The transition functions relating the dif- 
ferent coordinates we have considered, as well that the domains of definition 
of these coordinate systems themselves, are determined by the gauge source 
functions, the field equations, and the initial data. 

We have seen that in the case where the data and the given lapse and shift 
are real analytic, we can answer our questions in a satisfactory way. However, 
the assumption is not satisfactory. This is not meant to say anything against 
analytic solutions. In fact, most of the 'exact solutions' which are the source 
of our intuition for general relativistic phenomena are (piecewise) analytic. 
However, we should not restrict to analyticity in principle. One reason is that 
it would be in conflict with one of the basic tenets of general relativity. Given 
two non-empty open subsets U, V of a connected space-time (M, g) such 
that no point of U can be connected by a causal curve with a point of V, any 
process in U should be independent of what happens in V. However, if the 
space-time is analytic, the field in V is essentially fixed by the behaviour of g 
in U. For instance, we would not be able to study the evolution of data h a b, 
Xab on a 3— manifold S where h a b is conformally flat in some open subset of 
S but not in another one. 

Therefore, Einstein's equations should allow us to discuss the existence 
and uniqueness of solutions, and also the continuous dependence of the lat- 
ter on the data (stability), in classes of functions which are C°° or of even 
lower smoothness. In other words, Einstein's equations should imply evo- 
lution equations for which the Cauchy problem is well posed (cf. p5[ , p9| , 
[ |9l| ) . Whether an initial value problem is well posed cannot be decided on the 
level of analytic solutions and with the methods used to prove the Cauchy- 
Kowalewskaya theorem. On this level there is no basic distinction between 
initial value problems based on spacelike and those based on timelike hy- 
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persurfaces, though in the latter case the stability property is known not to 
be satisfied. Thus we are led to search for evolution equations which satisfy 
some 'hyperbolicity condition', i.e. a condition (essentially) on the algebraic 
structure of the equations which entails the well-posedness of the Cauchy 
problem. 

A number of different hyperbolicity conditions are known, all of them 
having in common that they require the equations to admit at each point a 
maximal number of 'r eal' characteristics: if the equations have a local ex- 



pression of the form (2.3), then the operator P is hyperbolic at x G U 
only if there is a covector f G T*M\ {0} such that every 2— dimensional 
plane in the cotangent space T* M containing £ intersects the conormal cone 
{det(J2^ a \ =m A a £ Q ) = 0} in k x m real lines (counting multiplicities) (cf. 
[ p7| ) . But notice that this condition alone does not ensure the well-posedness 
of the Cauchy problem. Some further remarks about different notions of hy- 



perbolicity can be found in Sect. 3.3 



To analyse the situation in the case of ( 2.22 ) and ( 2.2S ) , we solve ( 2.22 ) 



for Xab and insert into ( 2.23 ) to obtain a system of second order for h ab which 
takes the form 

\ {d 2 t h ab - (3 C d c 8t h ab - (3 C d t d c h ab + (3 C (3 d d c d d h ab } (2.28) 
cr 

-h cd (d c d d h ab + d a d b h cd - d a d c h bd - d b d c h ad ) 

2 2 

d a d b a (h c{a db) d t [i c - (3 d h c r a d b) d d f3 c ) 

= terms of lower order in h ab , a, [3 C . 

To analyse the characteristic polynomial, we have to know how a and f3 c 
are related to the solution h ab . Suppose, for simplicity, that a and (3 C are 
given functions.. Then we have to calculate for given covector ^ the 
determinant of the linear map 

Kb >-> Kb = -g^ f M £ v Kb - h cd k cd H a (, b + 2 C c K (a £ b) , 

of symmetric tensors, where we set ^ a = h ab £ b . Denoting by A'(£) the linear 
transformation which maps the independent components k ab , a < 6, onto the 
Kb, a < b, we find 

det A'(£) = -K^) 6 (<r^) 3 . 

Thus, if we consider lapse and shift as given, the conormal cone of system 
( [2.26 ) satisfies the condition required by hyperbolicity. Moreoever, the char- 



acteristics are timelike or null as one would expect for the evolution equations 
of a theory which is founded on the idea that physical processes propagate 
on or inside the light cone of the metric field. 

However, in spite of the fact that these equations have been used for a 
long time and in various contexts, and in spite of the naturality of equation 
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(2.2S) which appears to be indicated by the hyperbolicity of the subsidiary 



equations, it is not known whether the Cauchy problem for equations (2.2t) 
is well posed. It appears that we need to make use of the constraint equations 
to obtain suitable evolution equations. 

Therefore we proceed along a different route. While the form of the con- 
straint equations induced on a given hypersurface is unique, there is a huge 
freedom to modify our evolution equations. We can try to bring the principal 



part of ( 2.2S ) into a suitable form by using the constraints, by assuming lapse 
and shift to be functionals of the metric, or by subjecting them to some equa- 
tions. Then gauge source functions of quite a different nature may appear. 
Our choice in this section is motivated by the following two observations ( fj5f , 
||, cf. also H)): 

(i) Suppose S is some spacelike hypersurface of some Lorentz space (M, g) 
and x a are coordinates on some open subset U of S. Let — F p (x v ) be 
four smooth real functions defined on R 4 . Then there exist coordinates x v 
on some neighbourhood of U in M with x° = 0, x a = x a on U and such 
that the Christoffel coefficients of g in these coordinates satisfy the relations 
r»'(x u ') = Ff'ix"'), where f"' = g x ' p ' r x > p >. 

Obviously, one can construct coordinates x v on some neighbourhood of U 
by solving Cauchy problems for the semi-linear wave equations V M \7 p x" — 
— F^ (x p ) with Cauchy data on U which are consistent with our require- 
ments. When the wave equations are expressed in these coordinates the re- 
lations above result. 

(ii) The 4— dimensional Ricci tensor can be written in the form 

Rp» = -\ 9 Xp 9^,\p + V(„r v) (2.29) 

+A " n g v s g xp r p s „ + 2r s \ g sp 9x(M r v) \, 



where we set 



r v = 9vil r», v = d»r v - r p x „ r A 



Thus, if we consider the F v as given functions, the Einstein equations take 
the form of a system of wave equations for the metric coefficients. 

Before we show that these observations lead to a short and elegant argu- 
ment for the well-posedness of the initial value problem for Einstein's equa- 
tions (assuming well behaved matter equations), we indicate how the form of 
equations ( 2.29| ) relates to our previous considerations. 



From the expressions ( p.7|) we get the relations 

d t a-a, a p a =a 2 ix-n^r"), (2.30) 



d t P a - P a ,t P b = a 2 ( 7 a - D a log a - h a v r»), (2.31) 
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which indicate that prescribing the functions F v may fix the evolution of 
lapse and shift. Writing the 3— Ricci tensor similarly to ([2.29|) 



1 



Tab = - g h cd h a b,cd + -D(a7f>) + 7c " a hf d h ce -f e 1 b + 2 7 C a e h CJ h d ( a 7fc) e / , 



d h cf i 



with 



7a = ^ah 7 = ^ab ^ 7c d = h° {h ac ^ d - - /l 0( J,a)i 



(2.32) 



we obtain ( 2.28| ) in the form 
1 



— (&t hub — 2 h a b,tc (3° + h a b,cd f3 c P d ) — h cd h abtC d 



-tf D (a [ h b)c (dt P c - 13 d d d /3 C )] + 2 D [alb) -2D a D b log a 
= terms of lower order in h ab , ct, (3 C . 
Using ( |2.31 ) to replace the terms of second order in the second line, we get 



■gTK 



-2D (a {h b)c h%r^} 



(2.33) 



+ _ [ Da a D b a + X ab (d t a - C fj a) + 2 D (a a h b)c (d t p c - (3 C >d (3 d ) 
+2 (3 C ,( h b )c.t + h c (a P d ,b) P c ,d - P d ,( £/3 h b ) d ) 

+ 2 {2 Xac Xb C - Xc C Xab - 7c d a kf d h Ce 7e f b 



1 



-2 lc d R h c f h d(a lb) e f + Xh ab + K {T ab — - T h ab ) 



which can be read as a wave equation for h ab . From ( 2.30 ) and ( 2.22 ) we get 
1 



(d t - (3 C d c ) 2 a = (d t - P c d c ) x + 2 a X - 5 a 2 X n v r v 



(2.34) 



+3 a 3 (n v r v f -a(9t- f3 c d c ) {n u r v ). 

Takin g the trace of ( 2.23 ) and using ( 2.22| ) and the Hamiltonian constraint 
Q , we get 



{d t - f3 c d c )x = D a D a a-a{xabX a ~^ + K-{p + h ab T ab ) \ , (2.35) 



whence, using (2.34), the wave equation 



{dt - (3 C d c ) 2 a~D a D a a = ~a { X ab X ab - A 



(2.36) 



+ K -{p + h ab T a b) + 5 a x n v r u - 3 a 2 (n v r u ) 2 + (d t - f3 c d c ) (n v r v ) 
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From ( |2.32 ), (2.22), and the momentum constraint (2.11) follows 



(d t - I3 C d c ) 7 a = h cd f3 a >cd + aD a X + 2a K j a (2.37) 

+2 X ac D c a- X D a a-2ax cd lc a d-[3 a 
which implies together with ( 2.30 ), ( 2.31 ), and ( 2.22| ) the wave equation 



\ (d t - (3 C d c f (3 a - h cd f3 a , cd = 2anf 



(2.38) 



-4(x ac 



X h ac ) D c a-2a ( X bc - X h bc ) lb a c - /3° , c 7 C + P a ,o D c log a 
-2 a n v r v ( 7 Q - D a log a - h a v F v ) - 2 a x h a v F v 
+D a (a n v r v ) - (d t - p c d c )(h a v r v ). 



Equations ( 2.33 ), ( 2.36 ), and ( 2.38 ) form a hyperbolic system for the fields 
h a b, a, (3 a if we consider the functions as given. We note that besides 
( |2.23j ) we used the constraints as well as (2.7) to derive this system. 

We return to the evolution problem. Suppose we are given smooth data 
h a b, Xab, P, jai i- e. a solution of the constraints, on some 3— dimensional 
manifold S, which, for simplicity, we assume to be diffeomorphic to R 3 and 
endowed with global coordinates x a . Following our previous considerations, 
we set M — R x S, denote by t — x° the natural coordinate on R and extend 
the x a in the obvious way to M. We embed our initial data set into M by 
identifying S diffcomorphically with {0} x S (the need for this embedding 
shows that it is in general not useful to restrict the choice of the coordinates 
x a ). 

We now choose four smooth real functions F^(x ) on M which will be 
assigned the role of gauge source functions in the following procedure. As sug- 



gested by (2.29) and the preceeding discussion, we study the Cauchy problem 



for the reduced equations 



1 



9 Xp 9nu,\ P + V (/1 F v) + r x % 9llS g x P r p 



(2.39) 



i 



-Xg^ + k (T^ - - Tg^). 



Since we do not want to get involved in this section with details of the mat- 
ter equations, we assume the Cauchy problem for them to be well posed 
and the energy-momentum tensor to be divergence free as a consequence of 
these equations (one of the remarks following formula (4.70) shows that the 
situation can be occasionally somewhat more subtle). 

To prepare Cauchy data for (2.3E), we choose on S a smooth positive 
lapse function a and a smooth shift vector field f3 a , which determine with 
the datum h a b the ADM representation (2.6) for <? M „ on S. Using now h a b, 
Xab, Oi, (3 a in equations (2.30), (2.31) (with F v replaced by g^^F^), and 



( |2.22 ) , we obtain the corresponding datum dt g^v on S 
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It will be shown in Sect. [| that the Cauchy problem for (2.39) and the 
initial data above is well posed. Thus there exists a smooth solution g^ v to it 
on some neighbourhood M' of S. To assure uniqueness, we assume that M' 
coincides with the domain of dependence of S in M' with respect to gp V . 



With the functions calculated from g^ v , the reduced equation (2.39) 
takes the form 

Gp V + A g^ v — k Tp U = V( p D y ) - - V p D p 

with Dp = Fp — Fp. To see that our gauge conditions are preserved under 
the evolution defined by ( [2.39| ) and our gauge source functions, we show that 
Dp = on M'. 



Using the Bianchi identity ( [2.1§| ) and our assumptions on the matter equa- 
tions, we get, by taking the divergence of the equation above, the subsidiary 
equation 



V p A, 



Since we used (2.3C) and (2.31) to determine the initial data, we know that 
Dp = on S. Moreover, equations ( |2.36 )and ( 2.38 ) may be used, on the one 



hand, to calculate dt from g^ Vl but they are, on the other hand, contained 
in ( 2.39 ), with Dp replaced by Fp. This implies that d t Dp = on S. Using the 



uniqueness property for systems of wave equations discussed in Sect. 5.1 



we 



conclude that Dp vanishes on M 1 and gp V solves indeed Einstein's equations 
(U) on M'. 

We refer to the process of reducing the initial value problem for Einstein's 
equations to an initial value problem of a hyperbolic system as a hyperbolic 
reduction. The argument given here was developed for the first time in [fl6f 
with the 'harmonicity assumption' 7^, = 0. We prefer to keep the complete 
freedom to specify the gauge source functions, because some important and 
complicated problems are related to this. Our derivation of the system (2.33), 
( [2.36 ), (2.38) illustrates the intricate relations between equation (2.23), the 
constraints, and the conservation of the gauge condition in our case. Though 
other such arguments may differ in details, the overall structure of all hyper- 
bolic reduction procedures are similar. 



2.5 Assumptions and Consequences 

If one wants to investigate single solutions or general classes of solutions 
to the Einstein equations by an abstract analysis of the hyperbolic reduced 
equations, there are certain properties which need to be 'put in by hand' and 
others which are determined by the field equations. 

The properties in the first class depend largely on the type of physical 
systems which are to be modelled and on the structure of the hypersurface 
on which data are to be prescribed. The latter could be spacelike everywhere 
as in the Cauchy problem considered above. One may wish in this case to 
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consider cosmological models with compact time slices, one of which is the 
initial hypersurface, one may wish to model the field of an isolated system, 
in which case the initial hypersurface should include a domain extending to 
infinity with the field satisfying certain fall-off conditions, or one may wish 
to consider initial hypersurfaces with 'inner boundaries' whose nature may 
depend on various geometrical and physical considerations. Other cases of 



interest are the 'initial-boundary value problem'(cf. Sect. 4.3), where data 
are given on a spacelike and a timclikc hypersurface which intersect in some 
spacelike 2-surface, the 'characteristic initial value problem' where data are 
given on two null hypersurfaces which intersect in some spacelike two-surface, 
or various other combinations of hypersurfaces. Finally, we have to make 
a choice of matter model which may introduce specific problems which are 
independent of the features of Einstein's equations considered above. All these 
considerations will affect the problem of finding appropriate solutions to the 
corresponding constraint equations and the nature of this problem depends 
to a large extent on the nature of the bounding hypersurface. We note that 
the fall-off behaviour of the data on spacelike hypersurfaces can only partly 
be specified freely, other parts being determined by the constraint equations. 

There are conditions which we just assume because we cannot do better. 
For instance, one may wish to analyse solutions which violate some basic 
causality conditions (cf. f3^]) or which admit various identifications in the 
future and past of the initial hypersurface. Such violations or identifications 
introduce compatibility conditions on the data which cannot be controlled 
locally. At present there are no techniques available to analyse such situations 
in any generality. In the context of the Cauchy problem we shall only consider 
solutions (M, g) of the Einstein equations which are globally hyperbolic and 
such that the initial hypersurface S is a Cauchy hypersurface of the solution, 
i.e. every inextendible non-spacelike curve in M intersects S exactly once 

HI- 

The class of properties which need to be inferred from the structure of 
the data and the field equations includes in general everything which has to 
do with the long time evolution of the field: the development of singularities 
and horizons or of asymptotic regimes where the field becomes in any sense 
weak and possibly approximates the Minkowski field. This does not, of course, 
preclude the possibility of making assumptions on the singular or asymptotic 
behaviour of the field and analysing the consistency of these assumptions 
with the field equations. But in the end we would like to characterize the 
long time behaviour of the field in terms of the initial data. 

3 PDE Techniques 

3.1 Symmetric Hyperbolic Systems 

In order to obtain results on the existence, uniqueness and stability of solu- 
tions of the Einstein equations it is necessary to make contact with the theory 
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of partial differential equations. A good recent textbook on this theory is 




One part of the theory which can usefully be applied to the Einstein equations 
is the theory of symmetric hyperbolic systems. This will be discussed in some 
detail in the following. The aim is not to give complete proofs of the results of 
interest in general relativity, but to present some arguments which illustrate 
the essential techniques of the subject. It is important to note that the use of 
symmetric hyperbolic systems is not the only way of proving local existence 
and uniqueness theorems for the Einstein equations. The original proof ]T6| 
used second order hyperbolic equations. Other approaches use other types 
of equations such as mixed hyperbolic-elliptic systems (see, e.g. ]2f|, Chap. 
10). The basic tools used to prove existence are the same in all cases. First 
a family of approximating problems is set up and solved. Of course, in order 
that this be helpful, the approximating problems should be easier to solve 
than the original one. Then energy estimates, whose definition is discussed 
below, are used to show that the solutions of the approximating problems 
converge to a solution of the original problem in a certain limit. 

Now some aspects of the theory of symmetric hyperbolic systems will be 
discussed. We consider a system of equations for k real variables which are 
collected into a vector-valued function u. The solution will be defined on an 
appropriate subset of R x R n . (The case needed for the Einstein equations 
is n = 3.) A point of R x R n will be denoted by (t, x). The equations are of 
the form: 

A°(t, x,u)d t u + A i (t,x,u)d l u + B(t, x,u) = 

This system of equations is called symmetric hyperbolic if the matrices A 
and A 1 are symmetric and A is positive definite. This system is quasi-linear, 
which means that it is linear in its dependence on the first derivatives. The 
notion of symmetric hyperbolicity can be defined more generally, but here 
we restrict to the quasi-linear case without further comment. It is relatively 
easy to prove a uniqueness theorem for solutions of a symmetric hyperbolic 
system and so we will do this before coming to the more complicated existence 
proofs. A hypersurface S is called spacelike with respect to a solution u of 
the equation if for any 1-form n a conormal to S (i.e. vanishing on vectors 
tangent to S) the expression A a (t, x, u)n a is positive definite. This definition 
has a priori nothing to do with the usual sense of the word spacelike in general 
relativity. However, as will be seen below, the two concepts are closely related 
in some cases. Define a lens-shaped region to be an open subset G of R x R n 
with compact closure whose boundary is the union of a subset So of the 
hypersurface t = with smooth boundary dSo and a spacelike hypersurface 
Si with a boundary which coincides with 9 So- It will be shown that if G 
is a lens-shaped region with respect to solutions u\ and U2 of a symmetric 
hyperbolic system and if the two solutions agree on So then they agree on all 
of G. This statement is subject to differentiability requirements. 

Consider a symmetric hyperbolic system where A a and B are C 1 functions 
of their arguments. Let u\ and U2 be two solutions of class C 1 . Let G be a 
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lens-shaped region with respect to u\ and u 2 whose boundary is the union 
of hypersurfaces So and S% as before. There exist continuous functions M a 
and N such that: 

A a (t,x, Ul ) - A a (t,x,u 2 ) = M a (t,x,u 1 ,u 2 )(u 1 - u 2 ) 
B(t,x,ui) - B(t,x,u 2 ) = N(t : x,ui,u 2 )(ui - u 2 ) 

This sharp form of the mean value theorem is proved for instance in . It 
follows that: 

A a {u 1 )d a {u 1 - u 2 ) + [M a ( Ul ,u 2 )(d a u 2 ) + N( Ul ,u 2 )}( Ul -u 2 ) = Q 

Here the dependence of functions on t and x has not been written out explic- 
itly. This equation for u\ — u 2 can be abbreviated to 

A a (u 1 )d a (u 1 - u 2 ) = Q(ui - u 2 ) 

where Q is a continuous function of t and x. It follows that 

d a {{e- kt A a {u 1 )(u 1 -u 2 ),u 1 -u 2 )) 

= e' kt ([~kA°( Ul ) + {d a A a )( Ul ) + 2Q]( Ul ~ ua),ui - u 2 ) 

for any constant k. Now apply Stokes' theorem to this on the region G. This 
gives an equation of the form I\ = Iq + Ig where Iq, I\ and Iq are integrals 
over 5*0, S\ and G respectively. Because S\ is spacelike, I\ is non- negative. 
If k is chosen large enough the matrix P = kA°(ux) — (d a A a )(ui) — 2Q is 
uniformly positive definite on G. This means that there exists a constant 
C > such that (v, P(t, x, Ui)v) > C(v,v) for all v S R k and (t,x) in G. 
Thus if ui — u 2 is not identically zero on G, then the volume integral Ig 
can be made negative by an appropriate choice of k. If u\ and u 2 have the 
same initial data then Iq = 0. In this case Iq + Iq is negative and we get a 
contradiction. Thus in fact u\ — u 2 = 0. 

This argument shows that locally in a neighbourhood of the initial hy- 
persurface 5* the solution of a symmetric hyperbolic system at a given point 
is determined by initial data on a compact subset of S. For any point near 
enough to S is contained in a lens-shaped region. In this context we would like 
to introduce the concept of domain of dependence for a symmetric hyperbolic 
system. There are problems with conflicting terminology here, which we will 
attempt to explain now. If a given subset E of S is chosen, then the domain 
of dependence D(E) of E is the set of all points of R x R n such that the value 
of a solution of the symmetric hyperbolic system at that point is determined 
by the restriction of the initial data to E. Comparing with the definition of 
the domain of dependence for the Einstein equations given in Sect. 2.4 we 
find a situation similar to that seen already for the term 'spacelike'. The ap- 
plications of the term to the Einstein equations and to symmetric hyperbolic 
systems might seem at first to be completely unrelated, but in fact there is a 
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close relation. This fact was already alluded to in Sect. 2.4. Even for hyper- 
bolic equations there is another ambiguity in terminology. The definition for 
symmetric hyperbolic systems conflicts with terminology used in some places 
in the literature p?[ |. Often the set which is called domain of dependence 
here is referred to as the domain of determinacy of E while E is said to be a 
domain of dependence for a point p if p lies in what we here call the domain 
of dependence. In the following we will never use the terminology of p7| . 
However, we felt it better to warn the reader of the problems than to pass 
over them in silence. 

The relationship of the general definition of the domain of dependence 
for symmetric hyperbolic systems adopted in this article with the domain of 
dependence in general relativity will now be discussed in some more detail. 
As will be explained later, the harmonically reduced Einstein equations can 
be transformed to a symmetric hyperbolic system by introducing the first 
derivatives of the metric as new variables. The comparison we wish to make is 
between the domain of dependence D(S) defined by this particular symmetric 



hyperbolic system and the domain of dependence D(S) as defined in Sect. 2.4. 
The first important point is that a hypersurface is spacelike with respect to 
the symmetric hyperbolic system if and only if it is spacelike in the usual 
sense of general relativity, i.e. if the metric induced on it by the space-time 
metric is positive definite. Given this fact, we see that the notion of lens- 
shaped region used in this section coincides in the case of this particular 
symmetric hyperbolic system derived from the Einstein equations with the 



concept introduced in Sect. 2.4 



The above uniqueness statement will now be compared for illustrative 
purposes with the well-known uniquess properties of solutions of the wave 
equation in Minkowski space. Here we see a simplified form of the situation 
which has just been presented in the case of the Einstein equations. The wave 
equation, being second order, is not symmetric hyperbolic. However, it can 
be reduced to a symmetric hyperbolic system by introducing first derivatives 
of the unknown as new variables in a suitable way. It then turns out that once 
again the concept of a spacelike hypersurface, as introduced for symmetric 
hyperbolic systems, coincides with the usual (metric) notion for hypersurfaces 
in Minkowski space. It is well known that if initial data for the wave equation 
is given on the hypersurface t = in Minkowski space, and if p is a point in 
the region t > 0, then the solution at p is determined uniquely by the data 
within the intersection of the past light cone of p with the hypersurface t = 0. 
Once the wave equation has been reduced to symmetric hyperbolic form this 
statement can be deduced from the above uniqueness theorem for symmetric 
hyperbolic systems. It suffices to note the simple geometric fact that any 
point in the region t > strictly inside the past light cone of p is contained 
in a lens-shaped region with the same property. Thus the solution is uniquely 
determined inside the past light cone of p and hence, by continuity, at p. 
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The existence of the domain of dependence has the following consequence. 
If two sets u\ and of initial data on the same initial hypersurface S coincide 
on an open subset U of S 1 , then the corresponding domains of dependence 
Di(U) and D2QJ), as well as the corresponding solutions u\ and ui on them, 
coincide. In other words, on D\(U) the solution u\ is independent of the 
extension of u\ outside U. It follows in particular that there is no need to 
impose boundary conditions or fall-off conditions on u\ 'on the edge of S" or 
some periodicity condition to determine the solution locally near a point of 
S. This is referred to as the localization property of symmetric hyperbolic 
systems since it means that the theory is not dependent on knowing how 
the initial data behave globally in space. In the following data on R n will 
be considered which are periodic in the spatial coordinates. This is equiv- 
alent to replacing R n as domain of definition of the unknown by the torus 
T n obtained by identifying the Cartesian coordinates in R n modulo 2n. In 
order to prove statements about solutions of a symmetric hyperbolic system, 
something must be assumed about the regularity of the coefficients. This will 
in particular imply that A is continuous. The continuity of A and the com- 
pactness of the torus then show together that A (t, x, u) is uniformly positive 
definite on any finite closed time interval for any given continuous function 
u defined on this interval. 

The general strategy of the existence proof we will present here will now 
be discussed in more detail. This proof is essentially that given in pcj . The 
first detailed proof of an existence theorem for general quasi-linear symmetric 
hyperbolic systems was given by Kato [ |60| using the theory of semigroups. 
The proof here uses less sophisticated functional analysis, but the basic pat- 
tern of approximations controlled by energy estimates is the same in both 
cases. For simplicity we will restrict to the special case where A is the iden- 
tity matrix. Once the proof has been carried out in that case the differences 
which arise in the general case will be pointed out. Assuming now that A is 
the identity, the equation can be written in the form: 

dtu — —A l (t, x, u)diU — B(t, x, u) 

This looks superficially like an ordinary differential equation in an infinite 
dimensional space of functions of x. Unfortunately, this point of view is 
not directly helpful in proving local existence. The essential point is that 
if A l (t, x,u)di is thought of as an operator on a space of functions of finite 
differentiability, this operator is unbounded. The strategy is now to replace 
the unbounded operator by a bounded one. In fact we even go further and re- 
place the infinite dimensional space by a finite dimensional one. This is done 
by the introduction of a mollifier (smoothing operator) at appropriate places 
in the equation. The mollifier contains a parameter e. The resulting family 
of equations depending on e defines the family of approximating problems 
referred to above in this particular case. The approximating problems can be 
solved using the standard theory of ordinary differential equations. It then 
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remains to show that the solutions to the approximating problems converge 
to a solution of the original problem as e tends to zero. 

A convenient mollifier on the torus can be constructed using the Fourier 
transform. We now recall some facts concerning the Fourier transform on the 
torus. If u is a continuous complex-valued function on T n then its Fourier 
coefficients are defined by: 

u{Cj = [ u{x)e- l{xA) dx 

Here £ is an element of Z n , i.e. a sequence of n integers. In fact this formula 
makes sense for any square integrable function on the torus. It defines a lin- 
ear mapping from the space L 2 (T n ) of complex square integrable functions 
on the torus to the space of all complex- valued functions on Z n . Its image 
is the space of square summable functions L 2 (Z n ) and the L 2 norm of u is 
equal to that of u up to a constant factor. In particular the Fourier transform 
defines an isomorphism of L 2 (T n ) onto L 2 (Z n ). Thus to define an operator 
on L 2 (T n ) it suffices to define the operator on L 2 (Z n ) which corresponds to 
it under the Fourier transform. These statements follow from the elementary 
theory of Hilbert spaces once it is known that trigonometric polynomials (i.e. 
finite linear combinations of functions of the form e l ^> x *>) can be used to 
approximate all continuous functions on the torus in the sense of uniform 
convergence. This is worked out in detail for the case n = 1 in |8(|, Chap. 
4. The case of general n is not much different, once the approximation prop- 
erty of trigonometric polynomials is known in that context. This follows, for 
instance, from the Stone- Weierstrass theorem (0, p. 122). 

Let 4> be a smooth real- valued function with compact support on R n which 
is identically one in a neighbourhood of the origin and satisfies < < 1 
for all £ and (j>{— £) = For a positive real number e let </> e (£) = e~"</>(£/e). 
If we identify Z n with the set of points in R n with integer coordinates it is 
possible to define functions (f> e on Z n by restriction. The mollifier J e is defined 

J e U = (f) t U 

for any square integrable complex-valued function u on the torus. It is a 
bounded self-adjoint linear operator on L 2 (T n ). In fact || J e || = I for all e. 
The symmetry condition imposed on (f> ensures that if u is real, J e u is also 
real. Note that, for given e, J e u is only non-zero at a finite number of points, 
the number of which depends on e. It follows that J f u is a trigonometric 
polynomial and that the space of trigonometric polynomials which occurs for 
fixed e and all possible functions u is finite dimensional. As a consequence the 
image of L^(T"), the subspace of L 2 {T n ) consisting of real- valued functions, 
under the operator J e is a finite dimensional space V e . As e tends to zero the 
dimension of V e tends to infinity and this is why these spaces can in a certain 
sense approximate the infinite-dimensional space L 2 ^(T n ). 

The original problem will now be approximated by problems involving or- 
dinary differential equations on the finite-dimensional spaces V e . If the initial 



The Cauchy Problem 27 



datum for the original problem is u°, the initial datum for the approximating 
problem will be u° = J e u°. The equation to be solved in the approximating 
problem is: 

d t u e = -J e [A i (t,x,u e )J e d i u e + B(t,x,u e )} =F(u e ) (3.1) 

Here the unknown u e takes values in the vector space V e . In order to ap- 
ply the standard existence and uniqueness theory for ordinary differential 
equations to this it suffices to check that the right hand side is a smooth 
function of u e provided A 1 and B are smooth functions of their arguments. 
Let us first check that F is continuous. If u n is a sequence of elements of V t 
which converges to u then J e diU n converges to J e diU. The sequence of func- 
tions A l (t, x, u n )JfdiU n -\-B(t, x, u n ) converges uniformly to A l (t, x, u)J e diU + 
B(t, x, u). To sec that F is continous it remains to show that if v n — > v uni- 
formly J e v n — ► J e v. If v n — ► v uniformly then the convergence also holds 
in the L 2 norm. This implies that v n — > v in L 2 (Z n ). Hence 4> c v n — > (j> e v 
pointwise. From this it can be concluded that J e v n — > J f v. This completes 
the proof of the continuity of F. Differentiability can be shown in a simi- 
lar way. The main step is to show that if u and v are elements of V e then 
lim s ^o s [A % (t, x, u + sv) — A l (t, x, u)] exists (in the sense of uniform conver- 
gence). By smoothness of A 1 it does exist and is equal to D 3 A(t, x, u)v, where 
Z?3 denotes the derivative with respect to the third argument. In this way the 
existence of the first derivative of F can be demonstrated and the explicit 
expression for the derivative shows that it is continuous. Higher derivatives 
can be handled similarly. Hence F is a smooth function from V e to V e . The 
standard theory of ordinary differential equations now gives the existence of 
a unique solution of the approximating problem for each positive value of 
e. Existence is only guaranteed for a short time T e and at this point in the 
argument it is not excluded that T e could tend to zero as e — > 0. It will now 
be seen that this can be ruled out by the use of energy estimates. 

The basic energy estimate involves computing the time derivative of the 
energy functional defined by 

£ = J \u £ \ 2 dx (3.2) 

Since the functions u e are smooth and defined on a compact manifold differ- 
entiation under the integral is justified. 

d£/dt = 2 J (u e ,d t u e ) 

= -2 J (u £ ,J £ [A l J e d t u £ + B}) 

Using the fact that J e commutes with the operators 9, and is self-adjoint, 
gives 

j(u e ,J e [A l J e d lU£ }) = J (J e u e ,A%J e u e ) 
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Using Stokes theorem and the symmetry of A 1 gives 

J {J t u t ,Ad l J t u t ) = - J (diA 1 J e u e , J € u e ) - J {A i J e diU e , J e u e ) 
It follows that 

J (u e , J^A 1 J e diU e ]) =-\J {d t A l J <i u <i ,J t u t ) 
Substituting this in equation (3.2) gives: 

d t £ = J ((SSsA*) J e « e - 2B, J £ u £ ) 

Now 

\\(d i A i )J e u e -2B\\ L * < ||S! i A*|| z -||« e || £a +2||fl|| £a 
Hence it follows by the Cauchy-Schwarz inequality that 

dt£ < \\d i A i \\ Laa £ + 2\\B\\ L 2£ 1 / 2 

This is the fundamental energy estimate. Note that this computation is closely 
related to that used to prove uniqueness for solutions of symmetric hyperbolic 
systems above. 

The existence proof also requires higher order energy estimates. To obtain 
these, first differentiate the equation one or more times with respect to the 
spatial coordinates. Higher order energy functionals are defined by 



E 

|osj<s 



/ {D a u e ) 2 dx 



The square root of the energy functional £ s is a norm which defines the 
Sobolev space H s . It is because of the energy estimates that Sobolev spaces 
play such an important role in the theory of hyperbolic equations. Differen- 
tiating the equation for ( |3.l| gives 

d t (D a u e ) = -J e [A i (t,x,u e )J e d i D a u e } - J t [D a (A(t,x,u e )J e diu) 

-A(t, x, u e )J e d t D a u e ] - J e D a (B(t, x, u e )) (3.3) 

Differentiating the expression for £ s with respect to time causes the qua ntit y 
dt(D a u e ) to appear. This can be substituted for using the equation (3.3). 
Stokes theorem can be used to eliminate the highest order derivatives, as in 
the basic energy estimate. To get an inequality for £ s similar to that derived 
above for £ it remains to obtain L 2 estimates for the quantities 

D a (A l (t, x, u e )J e d % u e ) - A l {t, x, u e )J e d % D a u e 

and 

D a B{t,x,Ut) 
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This can be done by means of the Moser inequalities, which will be stated 
without proof. (For the proofs see |9(J.) 

The first Moser estimate says that there exists a constant C > such 
that for all bounded functions /, g on T 3 belonging to the Sobolev space H s 
the following inequality holds: 

\\D a (fg)\\ L 2 < C{\\f\\ L ~\\D s g\\ L . + \\D s f\\ L 4g\\ La ») 

Here s = \a\ and for a given norm ||-D S <?|| is shorthand for the maximum of 
||Z? Q g|| over all multiindices a with \a\ — s. The second estimate says that 
there exists a constant C > such that for all bounded functions /, g on T 3 
such that the first derivatives of g are bounded, / is in H s and g is in 
the following inequality holds: 

\\D a (fg) - fD a g\\ L2 < C(\\D s f\\ L , \\g\\ L ~ + \\D fU^D*- 1 g\\ L .) 

The third estimate concerns composition with nonlinear functions. Let F be 
a smooth function defined on an open subset of R k and taking values in R k . 
There exists a constant C > such that for all functions / on T 3 taking 
values in a fixed compact subset K of U and belonging to the Sobolev space 
H s and any multiindcx a with s — \a\ > 1 the following inequality holds: 

\\D a {F{f))\\ L ,<C\\Df\\ s L -2\\D s f\\ L , 

The result of using the Moser inequalities in the expression for d£ s /dt is 
an estimate of the form: 

d t £ s < C£ s 

where the constant C depends on a compact set K in which u e takes values, as 
above, and the L°° norm of the first derivatives of u e . This can be estimated 
by a C 1 function of the C 1 norm of u e . By the Sobolev embedding theorem 
(see e.g. |||), there is a constant C such that ||u e ||ci < C||M e ||_ff s for any 
s > n/2 + 1. Thus we obtain a differential inequality of the form 

d t £ s < f(£ s ) 

for a C 1 function /. It follows that the quantity £ s (t) satisfies the inequality 

£ s (t) < z(t) 

where z(t) is the solution of the differential equation dz/dt = f(z) with initial 
value £ (0) at t = 0. (A discussion of comparison arguments of this type can 
be found in pH , Chap. 3.) Since the function z remains finite on some time 
interval [0, T], the same must be true for £ s for any s > n/2 + 1 and the 
solutions u c exist on the common interval [0, T] for all e. Morover ||ti e (t)|| #s 
is bounded independently of e for all t 6 [0, T\. Putting this in the equation 
shows that ||<9tit e (t)||#s-i is also bounded. 
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The proof of existence of a solution with the given initial datum now 
follows by some functional analysis. The main point is that the boundedness 
of a sequence of functions in some norm often implies the existence of a 
sequence which is convergent in another topology. The reader who does not 
possess and does not wish to acquire a knowledge of functional analysis may 
wish to skip the rest of this paragraph. The space defined by the L°° norm 
of || it (t) 1 1 jjs is a Banach space and is also the dual of a separable Banach 
space. The family u e is bounded in this space and hence, by the Banach- 
Alaoglu theorem (see |^l), pp. 68-70) there is a sequence Ck tending to zero 
such that u tk (t) converges uniformly in the weak topology of H s (T n ) to some 
function u which is continuous in t with values in H s (T n ) with repect to the 
weak topology. The same argument can be applied to the family dtU e . As a 
result, after possibly passing to a subsequence, dtu ek (t) converges uniformly 
in the weak topology of H s ~ 1 (T n ). On the other hand, since T n is compact, 
the inclusion of H s (T n ) in H s (T m ) is compact. (This is Rellich's Lemma. 
See e.g. 47 , p. 88) Thus, by the vector- valued Ascoli theorem 29 it can 



be assumed, once again passing to a subsequence if necessary, that u ek (t) 
converges uniformly in the norm topology of H s ~ 1 (T n ). If s > n/2 + 2 then 
the Sobolev embedding theorem can be used to deduce that it converges 
uniformly in the norm of C 1 (T"). Once this is known it follows directly that 
the expression on the right hand side of the approximate equation converges 
uniformly to that on the right hand side of the exact equation. On the other 
hand dtu^ converges to dtu in the sense of distributions and so u satifies the 
original equation in the sense of distributions. The equation then implies that 
the solution is C 1 and is a classical solution. It also has the desired initial 
datum. 



Remarks 

1. The existence and uniqueness theorem whose proof has just been sketched 
says that given an initial datum uq in H s (T n ) with s > n/2 + 2 there is 
a solution for this datum which is a bounded measurable function u{t) with 
values in H s (T n ) and is such that dtu is a bounded measurable function with 
values in H s ~ 1 {T n ). With some further work it can be shown that n/2 + 2 
can be replaced by n/2 + 1 and that the bounded measurable functions are 
in fact continuous with values in the given space |jo) . 

2. The time of existence of a solution which is continuous with values in 
H s {T n ) depends only on a bound for the norm of the initial data in H s (T n ). 

3. The analogues of the higher energy estimates whose derivation was sketched 
here for the solutions of the approximating problems hold for the solution it- 
self This implies that as long as the C 1 norm of u(t) remains bounded 
the H s norm also remains bounded. This leads to the following continuation 
criterion. If u(t) is a solution on an interval [0, T) and if ||u(i)||ci is bounded 
independently of t by a constant C > then the solution exists on a longer 
time interval. For there is some T\ such that a solution corresponding to any 
data mq with ||uo||ij« < C exists on [0, Ti). Considering setting data at times 
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just preceding T, and using uniqueness to show that the solutions obtained 
fit together to give a single solution shows the existence of a solution of the 
original problem on the interval [0, T + T\). 

4. The previous remark implies the existence of C°° solutions corresponding 
to C°° data. For suppose that a datum uo of class C°° is given. Then uq 
belongs to H s (T n ) for each s. Thus there is a corresponding solution on an 
interval [0, T s ) which is continuous with values in H s . Assume that T s has 
been chosen as large as possible. It is not possible that T s — ► as s — > oo. For 
the boundedness of the solution in H s , s > n/2 + 1 implies its boundedness 
in C 1 and this in turn implies that the solution can be extended to a longer 
time interval. Thus in fact T s does not depend on s. Using the equation then 
shows that the solution is C°° . 

5. Analogous statements to all of the above can be obtained for the case 
where A is not the identity with similar techniques. The essential point is to 
modify the expressions for the energy functionals using A . For instance the 
basic energy for the approximating problems is given by £ — f (A°u e , u e ). 

6. Our treatment here differs slightly from that of Taylor |MJ by the use of 
a reduction to finite-dimensional Banach spaces. Which variant is used is a 
matter of taste. 

The statement made in the second remark is a part of what is referred to 
as Cauchy stability. This says that when the initial data is varied the corre- 
sponding solution changes in a way which depends continuously on the data. 
This assertion applies to solutions defined on a fixed common time interval 
[0, T]. The above remark shows that a common time interval of this kind can 
be found for data which are close to a given datum in H s (T n ). A related 
statement is that if the coefficients of the equation and the data depend 
smoothly on a parameter, then for a compact parameter interval the corre- 
sponding solutions for different parameter values exist on a common time 
interval [0, T] and their restrictions to this time interval depend smoothly on 
the parameter. 

3.2 Symmetric Hyperbolic Systems on Manifolds 

In the last section it was shown how existence and uniqueness statements 
can be obtained for solutions of symmetric hyperbolic systems with data 
on a torus. It was also indicated that these imply results for more general 
settings by means of the domain of dependence. In this section some more 
details concerning these points will be provided. First it is necessary to define 
what is meant by a symmetric hyperbolic system on a manifold M. (In the 
context of the last section M — T 3 x /.) In general this will be an equation 
for sections of a fibre bundle E over M which, for simplicity, we will take 
to be a vector bundle. The easiest definition is that a symmetric hyperbolic 
system is an equation of the form P(u) — 0, where P is a nonlinear differential 
operator on M which in any local coordinate system satisfies the definition 
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given previously. For the general invariant definition of a differential operator 
on sections of a bundle the reader is referred to |ffl[ . The object defined 
in local coordinates by er(£) = A a t; a is called the principal symbol of the 



differential operator (cf. Sect. 2.1). Here A a is supposed to be evaluated on a 
given section u. The principal symbol is defined invariantly by the differential 
operator. Let L(V) be the bundle of linear mappings from V to itself. Then 
the invariantly defined principal symbol is a section of the pullback of the 
bundle L(V) to the cotangent bundle T*M. Given a local trivialization of 
V it can be identified locally with a k x k matrix-valued function on the 
cotangent bundle. Here k is the dimension of the fibre of V, i.e. the number 
of unknowns in the local coordinate representation. 

In order to define the condition of symmetry in the definition of symmetric 
hyperbolic systems in an abstract context, it is necessary to introduce a 
Riemannian metric on the bundle V. This metric is not visible in the local 
coordinate representation, since in that context this role is played by the flat 
metric whose components in the given coordinates are given by the Kronecker 
delta. The symmetry condition is expressed in terms of the principal symbol 
and the chosen inner product as: 

(a(£)v,w) = (v,a(£)w) 

for all sections v, w of V and all covectors £. The positivity of A is replaced by 
the condition that for some covector £ the quadratic form associated to <r(£) 
via the metric on V (by lowering the index) is positive definite. A hypersurface 
such that all its non-vanishing conormal vectors have this property is called 
spacelike. 

Suppose now we have a symmetric hyperbolic system on a manifold M, 
a submanifold S of M and initial data such that S is spacelike. The aim is 
to show the existence of a solution on an open neighbourhood U of S in M 
with the prescribed data and that this is the unique solution on U with this 
property. Choose a covering of S by open sets U a with the property that for 
each a the closure of this set is contained in a chart domain over which the 
bundle V can be trivialized globally. The problem of solving the equations 
on some open region in M with data the restriction of the original initial 
data to U a can be solved by means of the results already discussed. For the 
local coordinates can be used to identify U a with an open subset U' a of a 
torus T n and the initial data can be extended smoothly to the whole torus. 
The equation itself can also be transferred and extended. Corresponding to 
the extended data there exists a solution of the equations on T" x I for 
some interval I. The coordinates can be used to transfer (a restriction of) 
the solution to an open subset W a of M which is an open neighbourhood 
of U a . The domain of dependence can then be used to show that there is 
an open neighbourhood where the solutions on the different W a agree on 
the intersections. More specifically, on an intersection of this kind the two 
in principle different solutions can be expressed in terms of the same local 
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coordinates and trivialization. By construction the coordinate representation 
of the data is the same in both cases. Thus the local uniqueness theorem 
for symmetric hyperbolic systems implies that both solutions are equal in 
a neighbourhood of the initial hypersurface. Thus the solutions on different 
local patches fit together to define the desired solution. Moreover the domain 
of dependence argument shows that it is unique on its domain of definition. 

3.3 Other Notions of Hyperbolicity 

There are two rather different aspects of the concept of hyperbolicity. The 
one is an intuitive idea of what good properties a system of equations should 
have in order to qualify for the name hyperbolic. The archetypal hyperbolic 
equation is the wave equation and so the desired properties are generalizations 
of properties of solutions of the wave equation. The first property is that the 
system should have a well-posed initial value problem. This means that there 
should exist solutions corresponding to appropriate initial data, that these 
should be uniquely determined by the data, and that they should depend 
continuously on the data in a suitable sense. The second property is the 
existence of a finite domain of dependence. This means that the value of 
solutions at a given point close to the initial hypersurface should depend 
only on data on a compact subset of the initial hypersurface. 

The first aspect of hyperbolicity which has just been presented is a list of 
wishes. The second aspect is concerned with the question, how these wishes 
can be fulfilled. The idea is to develop criteria for hyperbolicity. For equations 
satisfying one of these criteria, the desired properties can be proved once and 
for all. Then all the user who wants to check the hyperbolicity of a given 
system has to do is to check the criteria, which are generally more or less 
algebraic in nature. An example is symmetric hyperbolicity, which we have 
already discussed in detail. 

Symmetric hyperbolicity can be applied to very many problems in gen- 
eral relativity, but there are also other notions of hyperbolicity which have 
been applied and cases where no proof using symmetric hyperbolic systems is 
available up to now. We will give an overview of the situation, without going 
into details. The concepts which will be discussed are strict hyperbolicity, 
strong hyperbolicity and Leray hyperbolicity. This multiplicity of definitions 
is a consequence of the fact that there is no one ideal criterion which covers 
all cases. 

There is unfortunately to our knowledge no fully general treatment of 
these matters in the mathematical literature. A general theory of this kind 
cannot be given in the context of this article, and we will refrain from stating 
any general theorems in this section. We will simply indicate some of the 
relations between the different types of hyperbolicity and give some references 
where more details can be found. 

The discussion will use some ideas from the theory of pseudodifferential 
operators. For introductions to this theory see psf, [|B9|, 0]. One important 
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tool is an operator A which defines an isomorphism from functions in the 
Sobolev space H s to functions in the Sobolev space Powers of this 

operator can be used to associate functions with a given degree of differentia- 
bility (in the sense of Sobolev spaces) with functions with a different degree of 
differentiability. The operator A is non-local. On the torus it could be chosen 
to be the operator (1 — Z\) 1 / 2 , defined in an appropriate way. 

The concept of strict hyperbolicity is easily defined. Suppose that we have 
a system of N equations of order m. Suppose there is a covector r such that 
for any covector £ not proportional to r the equation deta(x, r + A£) has 
Nm distinct real solutions A. In this case we say that all characteristics are 
real and distinct and that the system is strictly hyperbolic. The weakness of 
this criterion is due to the fact that it is so often not satisfied for systems of 
physical interest. For example, the system for two functions u and v given by 
the wave equation for each of them is not strictly hyperbolic. The wave equa- 
tion itself is strictly hyperbolic but simply writing it twice side by side leads 
to characteristics which are no longer distinct, so that strict hyperbolicity is 
lost. 

Given a strictly hyperbolic system of order m with unknown it, let Ui — 
A~ % d\u for % — 1, 2, .... m — 1. Then the following system of equations is 
obtained: 



dfUi = Aui+i; < i < m — 1 
d t u m - 1 =A- m+ \d?u) 



where it is understood that d™u should be expressed in terms of uq, . . . , u m -i 
by means of the original system. This system is first order in a sense which is 
hopefully obvious intuitively and which can be made precise using the theory 
of pseudodifferential operators. This is not a system of differential equations, 
due to the nonlocality of A, but rather a system of pseudodifferential equa- 
tions. Suppose we write it as dtv + A l (v)diV + B(v) = 0, where v is an 
abbreviation for uq, . . . ,u m -i. The idea is now to find a pseudodifferential 
operator A°(v) of order zero depending on the unknown so that A°(v)A 1 (v) 
is symmetric for all v and i and A has suitable positivity properties. If an A 
of this kind can be found, then multiplying the equation by it gives something 
which looks like a pseudodifferential analogue of a symmetric hyperbolic sys- 
tem. This symmetrization of the reduced equation can in fact be carried out 
in the case of strictly hyperbolic systems and the resulting first order pseu- 
dodifferential equation be treated by a generalization of the methods used 
for symmetric hyperbolic systems , |9^] . This gives an existence theorem, 
but does not directly give information about the domain of dependence. This 
information can be obtained afterwards by a different method [§8| . 

The algebra involved in symmetrizing the reduction of a strictly hyper- 
bolic system has not been detailed. Instead we will present some general 
information about symmetrization of first order systems using pseudodiffcr- 
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ential operators. Consider then the equation 

dtU + A^ujdiU = 

If there exists a function A°(u) such that A is positive definite and A A 1 
is symmetric for each i then multiplying the equation by A°(u) gives the 
equation: 

A"(u)d t u + A°{u)A i (u)d i u = 

which is symmetric hyperbolic. This can be generalized as follows. Each dif- 
ferential operator has a symbol which is polynomial in £. In the theory of 
pseudodifferential operators, operators are associated to symbols which no 
longer have a polynomial dependence on £. Within this theory the general 
statement (which must of course be limited in order to become strictly true) 
is that any algebra which can be done with symbols can be mimicked by op- 
erators. In terms of symbols the problem of symmetrization which has been 
stated above can be formulated as follows. Given symbols (&,«), which 
are matrix-valued functions, find a positive definite symmetric real matrix- 
valued function a°(x,u) such that a°(x, it)^a l (x, u) is symmetric. When the 
problem has been formulated in this way a generalization becomes obvious. 
Why not allow a to depend on £? In the context of pseudodifferential op- 
erators this can be done. In order for this to be useful a must be of order 
zero, which means that it must be bounded in £ for each fixed (x, u) and 
that its derivatives with respect to £ and x satisfy similar conditions which 
will not be given here. In fact it is enough to achieve the symmetrization of 
the symbol for £ of unit length, since the £ dependence of the procedure is 
essentially homogeneous in £ anyway. It is also important that it be possible 
to choose a in such a way that it depends smoothly on £, since the theory of 
pseudodifferential operators requires sufficiently smooth symbols. To prevent 
confusion in the notation, a will from now on be denoted by r while we write 
a(x, £, u) = a l (x, u)&. 

A criterion which ensures that a first order system of PDE can be sym- 
metrized in the generalized sense just discussed is that of strong hyperbolicity 
[ p3| . It is supposed that all characteristics are real, that they are of constant 
multiplicity, and that the symbol is everywhere diagonalizable. Furthermore, 
it is assumed that this diagonalization can be carried out in a way which 
depends smoothly on x, £ and u. This means that there is a symbol b(x, £, it) 
which satifies 

b^ 1 (x, £, u)a(x, £, u)b(x, £, u) = d(x, £, u) 

where d is diagonal. Let r(x,£,u) = [b(x, £, u)b T (x, £, u)]~ 1 where the su- 
perscript T denotes the transpose. Then r has the desired properties. For 
r = fe _1 (& _1 ) T is symmetric and positive definite while ra = (b~ 1 ) T db~ 1 is 
also symmetric. The one part of this criterion which is not purely algebraic 
is the smoothness condition. It is, however, often not hard to check once it is 
known how to verify the other conditions. 
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Another related notion of hyperbolicity is Leray hyperbolicity ]66[ . It al- 
lows operators with principal parts of different orders in the same system. We 
will discuss this in terms of using powers of the operator A to adjust the or- 
ders of the equations. This is somewhat different from the approach of Leray 
but is closely related to it. Suppose we have a system of PDE with unknown 
u. Let us split the vector u into a sequence of vectors u\, . . . ,Ul. Then the 
differential equation which we may write schematically as P(u) = can be 
rewritten as Pj(ui) = 0, where i and j run from I to L. Let Vi = A s ^Ui mul- 
tiply the equations Pj(ui) = with A l ^\ The result is a system of equations 
of the form Qj{vi) = 0, where Qj(Vi) = A t( -^P 3 {A-^Vi). The order of Q 3 in 
its dependence on V{ is the order of Pj in its dependence on Ui plus t(J) — s(i). 
Now we would like to choose t(J) and s(i) in such a way that all operators Qi 
have the same order (say one) in their dependence on the corresponding vari- 
able Vi and lower order (say zero) in their dependence on Vj for i ^ j. Adding 
the same amount to all indices simultaneously is irrelevant. The indices are 
only determined up to a common additive constant. The system is Leray hy- 
perbolic if this can be achieved by suitable choices of the decomposition and 
the indices, if the operators Pi corresponding to the given decomposition, 
considered with respect to their dependence on Ui, are strictly hyperbolic, 
and if the characteristics of these strictly hyperbolic operators satisfy a cer- 
tain condition. (We say more on this condition later.) It is now plausible 
that this decomposition can be combined with the symmetrization of strictly 
hyperbolic systems to obtain a system of first order pseudodifferential equa- 
tions which admit a pseudodifferential symmetrization. We do however stress 
that as far as we know the only place where the existence proof for Leray 
hyperbolic systems is written down in the literature is in the original lecture 
notes of Leray |36| and that the details of the plausibility argument presented 
here have not been worked out. It has the advantage that it gives an intuitive 
interpretation of the meaning of the indices s(i) and t(j). All the components 
of v in the solution will have the same degree of differentiability. If this is H k 
then the variables Ui will have differentiability H k ~ s ^\ 

We now comment on the condition on the characteristics mentioned in 
the discussion of Leray hyperbolicity. For a strictly hyperbolic system there 
is a notion of spacelike covectors analogous to that for symmetric hyperbolic 
systems. The position of these spacelike covectors is closely related to the 
position of the characteristics. The required condition is that there should 
be vectors which are simultaneously spacelike for the L strictly hyperbolic 
systems occurring in the definition. 

These general ideas concerning Leray hyperbolic systems will now be il- 
lustrated by the example of the Einstein-dust system. This example was first 
treated by Choquet-Bruhat JlTj . A discussion of this and other examples can 
also be found in the books of Lichnerowicz |37), Jss| . We will have more to 
say about this system in Sects. 4.2 and 5.4. The Einstein-dust equations are 
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V a {pU a )=0 
U a V a UP = 

The problem with these equations from the point of view of symmetric hy- 
perbolic systems is that while the equation for the evolution of p contains 
derivatives of U a the equation for the evolution of U a does not contain p 
This could be got around if the derivatives of U a in the evolution equation 
for p were considered as lower order terms. This is only possible if the first 
derivatives of U a are considered on the same footing as p. The adjustment 
of orders involved in the Leray theory combined with an extra device, allows 
this to be achieved. In order that everything be consistent we expect the order 
of differentiability of the metric to be one greater than that of U a . For the 
Christoffel symbols occur in the evolution equations for U a . Combining these 
two things means that the density should be two times less differentiable than 
the metric. This creates problems with the Einstein equations. For these are 
essentially (and the harmonically reduced equations are precisely) non-linear 
wave equations for the metric. The solution of a system of this kind is only 
one degree more differentiable than the right hand side. This does not fit, 
since the density occurs on the right hand side. The extra device consists 
is differentiating the Einstein equations once more in the direction U a and 
then substituting the evolution equation for p into the result. This gives the 
equation: 

J7 7 V 7 G Q/3 = SirpUcUpV^ 

Note that the right hand side of this contains no derivatives of p and so is not 
worse than the right hand side of the undifferentiated equation. On the other 
hand, the left hand side is, in harmonic coordinates, a third order hyperbolic 
equation in g a p and the solution of an equation of this type has (by the Leray 
theory) two more degrees of differentiability than the right hand side. 

After this intuitive discussion of the Einstein-dust system, let us show how 
it is related to the choice of indices needed to make the Leray hyperbolicity of 
the system manifest. The first thing which needs to be done is to decide which 
equations should be considered as the evolution equations for which variables. 
This has already been done in the above discussion. The relative orders of 
differentiability we have discussed suggest that we choose s(g) = 1, s(U) = 2 
and s(p) = 3. This equalizes the expected differentiability of the different 
variables. In the new variables the evolution equations for g, U and p are of 
order four, three and four respectively. To equalize the orders we can choose 
t{g) = 1, t(U) = 2 and t(p) = 1. The blocks into which the system must 
be split are not just three but fifteen, one for each component. The indices 
s(i) and t(j) are here chosen the same for each component of one geometrical 
object. (Here we have ignored the difficulty that because of the normalization 
condition U a U a = —1 the variables U a are not all independent.) 
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4 Reductions 

In Sect. U we saw an example of a hyperbolic reduction for Einstein's field 
equations which is, at least in the vacuum case, sufficient to obtain local 
existence results and to demonstrate that the local evolution is dominated 
by the light cone structure and the associated concept of the domain of 
dependence. Thus, besides the existence problem, it settles the conceptual 
question whether the evolution process determined by Einstein's equations is 
consistent with the basic tenets of the theory. 

Nevertheless, there are various reasons for considering other types of re- 
ductions. Different matter fields may require different treatments, various 
physical or geometrical considerations may require reductions satisfying cer- 
tain side conditions, the desire to control the gauge for an arbitrarily long 
time may motivate the search for new gauge conditions. In particular, in re- 
cent years various systems of hyperbolic equations deduced from Einstein's 
equations have been put forward with the aim of providing 'good' systems for 
numerical calculations (cf. @,g,0,@,@,@,[f6j). Since it is difficult 
to judge the relative efficiency of such systems by a few abstract arguments, 
detailed numerical calculations are needed to test them and it still remains 
to be seen which of these systems will serve the intended purpose best. 

A general discussion of the problem of finding reductions which are use- 
ful in numerical or analytical studies should also include systems combining 
hyperbolic with elliptic equations (cf. pj|, |77| ], |24| ] for analytic treatments 
of such situations) or with systems of still other types. For simplicity we will 
restrict ourselves to purely hyperbolic reductions. But even then a general 
discussion does not yet appear feasible. Ideally, one would like to exhibit a 
kind of 'hyperbolic skeleton' of the Einstein equations and a complete charac- 
terization of the freedom to fix the gauge from which all hyperbolic reductions 
should be derivable. Instead, there are at present various different methods 
available which have been invented to serve specific needs. 

Therefore, we will present some of these methods without striving for 
completeness. There are various different boundary value problems of inter- 
est for Einstein's equations and tomorrow a new question may lead to a new 
solution of the reduction problem. Our aim is rather to illustrate the enor- 
mous richness of possibilities to adapt the equations to various geometrical 
and physical situations, and to comment on some of the new features which 
may be observed. Apart from the illustrative purpose it should not be forgot- 
ten, though, that each of the reductions outlined here implies, if worked out 
in detail, the local existence of solutions satisfying certain side conditions. 

We begin by recalling in general terms the steps to be considered in a 
reduction procedure. It should be noted that these steps are not independent 
of each other, and that the same is true of the following considerations. 

As a first step we need to choose a representation of the field equations. 
Above we used the standard representation in which Einstein's equations are 
written as a system of second order for the metric coefficients g^ u and we 
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also considered the ADM equations with the first and second fundamental 
forms of the time slicing as the basic variables. Later we shall also employ 
representations involving equations of third order in the metric field, in which 
the metric is expressed in terms of an orthonormal frame field. 

Contrary to what has been claimed by too ambitious or suggested by am- 
biguous formulations, there is no way to get hyperbolic evolution equations 
without fixing a gauge. Hyperbolicity implies uniqueness in the PDE sense 
(in contrast to the notion of geometric uniqueness used for the Einstein equa- 
tions, cf. Sect. |5.2|) . It is necessary to make a choice of precisely four gauge 
source functions. 

Although in the proof of an existence result one has to fix a coordinate 
system, it is of interest to note that there are choices of gauge conditions 
which render the reduced equations hyperbolic for any fixed choice of co- 
ordinates. In Sect. H we considered the functions F v (x^) as gauge source 
functions determining the coordinates. Since these functions can be chosen 
arbitrarily and since the principal part of the reduced equations ( [2.39] ) will 
not be changed if they depend on the metric coefficients, they can be assumed 
to be of the form F v = F"(x»,gx p ). 

The following choice is particularly interesting. Let an affine connection 
r u M a be given on the manifold M = R x S on which we want to construct 
our solution. Since the difference r u \ — T v \, where the first term denotes 
the connection of the prospective solution g^ Vl defines a tensor field, the 
requirement that the equation _T M = g vX F v ^ \ be satisfied, has an invariant 
meaning. This suggests a way to impose a gauge condition which removes the 
freedom to perform diffeomorphisms while leaving all the freedom to perform 
coordinate transformations. In particular, if we assume the connection T v \ 
to be the Levi-Civita conn ection of a metric on M, the condition = 



gvX y v /* A in equation ( 2.39 ) corresponds to the requirement that the identity 
map of M onto itself is a harmonic map from (M, g) to (M,g) (cf. also J37J, 

Again, we do not have a good overview of all the possibilities to impose 
gauge conditions. There exist conditions which work well with quite different 
representations of the field equations. Examples are given by the choice of 
gauge source functions F w (x^) considered above, which work with suitably 
chosen gauge conditions for the frame field also in the frame formalism ( |m| ) , 
or by the gauge, considered also below, in which the shift f3 a and the function 
q = log(a/i _<T ), with h = det(h a b) and a — const. > 0, are prescribed as 
gauge source functions (cf. 0, Q for reductions based on a — | and 
different representations of the Einstein equations). There are other gauge 
conditions which only work for specific representations like some of the ones 
we shall consider in the context of the frame formalism. For some gauge 
conditions, like the ones employing the gauge source functions F^(x^), the 
universal applicability has been shown (cf. our argument in Sect. H) or is easy 
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to see. For others, like the ones using the gauge source function q above, it 
has apparently never been shown. 

The gauge problem does not admit a 'universal solution' which works in all 
possible situations of interest. In fact, choosing the gauge is related to some of 
the most complicated questions of constructing general solutions to Einstein's 
equations. Often one would like to find a system of coordinates which covers 
the complete domain of existence of a solution arising from given initial data. 
If there existed such coordinates x^, we could, in principle, characterize them 
in terms of the associated gauge source functions F v (x^) considered in Sect. ||. 
However, the domain of validity, in particular the 'lifetime', of a coordinate 
system depends on the data, the equations, the type of gauge condition, 
as well as on the given gauge source functions. In practice, gauge source 
functions which ensure that the coordinates exist globally, if there exist any 
at all, have to be identified in the course of constructing the solution. 

Having implemented the gauge conditions, we have to find a hyperbolic 
system of reduced equations from our representation of the Einstein equations. 
As we shall see, there are often various possibilities and the final choice will 
depend on the desired application. With the reduced equations at hand we 
have to arrange initial data which satisfy the constraints and are consistent 
with the gauge conditions. While the second point usually poses no prob- 
lem, the first point involves solving elliptic equations and requires a seperate 
discussion. 

The reduced equations define together with the initial data a well-posed 
initial value problem and we can show, by using standard techniques of the 
theory of partial differential equations (cf. Sect. [O]) , the existence of solutions 
to the reduced problem, work out differentiability properties of solutions etc. 
or start the numerical evolution. 

As the final step one needs to show that the constraints and gauge con- 
ditions are preserved by the evolution defined by the reduced equations. One 
may wonder whether there exists a universal argument which tells us that this 
will be the case for any system of hyperbolic reduced equations deduced from 
the Einstein equations. The example of the spinor equations discussed below 
shows that this cannot be true without restrictions on the matter fields. How- 
ever, even if we ask this question about reduced problems for the Einstein 
vacuum field equations, the answer seems not to be known. The standard 
method here is to use the reduced equations, as in Sect. ||, together with 
some differential identities, to show that the 'constraint quantities' satisfy 
a certain subsidiary system which allows us to argue that these quantities 
vanish. 

If the gauge condition used in our reduced equations has been shown to be 
universally applicable and if there already exists another hyperbolic reduction 
using these gauge conditions and applying to the same geometric and physical 
situation (choice of matter model), it can be argued, invoking the uniqueness 
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property of initial value problems for hyperbolic equations, that our reduced 
equations must preserve the constraints and the gauge conditions. 

In the examples which follow we shall not always produce the complete re- 
duction argument. Often we shall only exhibit a symmetric hyperbolic system 
of reduced equations and remark on some of its properties. 



4.1 Hyperbolic Systems from the ADM Equations 

In the article pl| (c f. also |8(|) the authors derived a system from equations 
( pTC| ), foil), (pT2^) and ( pT23[ ) which seems to be numerically distinctively 



better behaved than any other system derived so far from the ADM equations. 
However, at present there appears to be no clear understanding as to why 
this should be so. The symmetric hyperbolic system to be discussed below 
was found in the course of an attempt to understand whether the system 
considered in jllj is in any sense related to a hyperbolic system. Since other 
hyperbolic systems related to the equations in (llj] have been discussed in Q , 
[fl4| there is now a large family of such systems available. 
In the following we shall consider the fields 

(3 a , g = log(a/ l - ff ), (4.1) 

with h = det(/i Q b), a = const. > 0, as the gauge source functions. The 
density h will be used to rescale the 3-metric and the trace-free part of the 
extrinsic curvature to obtain the densities h a b — h a b, \ab — h~s (x a b — 
5 X h a b)- For simplicity we shall refer to h a b as to the conformal metric or, if 
no confusion can arise, simply as to the metric. The following equations are 
derived from the ADM equations by using the standard rules for conformal 
rescalings. The occurrence of some of the terms in the following equations 
find their proper explanation in the general calculus for densities but we shall 
not discuss this here any further. 

The unknowns in our equations will be the fields 

Kb, V = log h, rja = D a log h, \ = h ab Xab, (4.2) 

7a = 7a -{-X +v)D a log h, Xab, habc = hab,c, (4.3) 

6 



where 7 a is defined for the metric h a b in analogy to (2.32). Note that the power 
of the scaling factor h has been chosen such that we have h — det(h a b) = 1, 
whence = h tC = hh ab h a b,c — h ab h a bc- We denote by D the covariant 
derivative, by j a b c the connection coefficients, and by R the Ricci scalar 
which are defined for the conformal metric h a b by the standard rules. In 
all expressions involving quantities carrying a tilde any index operations are 
performed with the metric h a b- 



Our equations are obtained as follows (cf. also [|ri|). From (2.22) we get 
by a direct calculation 
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_ _ 2 - 

d t h ab - h ab , c (3 C = 2 h c{a (3 C >b) + - h ab P c ,c + 2 a Xab, (4.4) 

d t r ] -r ltC l3 c = 2l3\ c + 2ax- (4.5) 
Taking derivatives on both sides of the second equation gives 



d t Va-Va,cl3 c = 2aD aX + 2 X (ar) a + D a q) + Vc (3 C , a + 2 (3 C , c 



(4.6) 



Equation (2.23) implies together with (2.22) that 



dtX-X.cF = D C D c a-a((l-a)R + aR + X 2 - 3 X + - (h ab T ab - 3 p)), 

with an arbitrary real number a. To replace R in the second term on the 
right hand side, we use the Hamiltonian constraint (2.10), to replace R in the 
third term we use the transformation law of the Ricci scalar under conformal 
rescalings and the expression of the Ricci scalar R in terms of the conformal 
quantities, which give 



U = h ~ \D C (% -~D C log h) - > c b l b - 



1 - 



+h ab j a c d 7 c d b -—D a log hD a \ogh 



Thus we get 



hi -(9tX~ X,c P c ) = -ab a % + {a + a{)-- a)) D a Va (4.7) 
a 2 

+D a D a q + a(% c d Zf d -h ab % c a % d b ) 

+ a + \° + ct2 ) Va V a + {\ + 2 <?) if D a q + D a q D a q 
18 o o 

-0 {(1 - a) Xab X ab + I (1 + 2 a) X 2 - (1 + 2 a) A + | {h ab T ab + (1-2 a) p)}. 

Using the equations above, the definition of and the momentum constraint 
( |2.11| ) as well as its expression in terms of the conformal quantities, 

2 ~ 1 1 - 

D c Xca ~ 3 D aX + - Xca D c log h + - X D a log h = KT IMU h v a , 

we obtain, with arbitrary real number c, by a direct calculation the equation 

- {dt % - %,c P c } = -cD c X ca + 2{\ + \- <J) D a x (4.8) 
a 2 3 

+ {2 + c)kT^ ri l h v a + 2xa c {lc - (77 + 7 -2 a) Vc + D c q} 
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-X ((- + 2 a) (a T) a + D a q) + \ r, a ) ~2x cd lc b d Kb 
3 6 

+ - \h cd h ac , b f3 h id - /3 d ) (h ab lc b d + Kb ld b a ) - h cd , b h ac . d (3» 



(o + °) Vb P ,« + h cd (h bc 13" |0 + hab 13" , c ),d -(1 + 2 a) 13" 



JjQ 



For the rescaled trace free part of the extrinsic curvature we get from (2.23) 
the equation 

2 

dt Xab ~ Xab.c P C ~ 2 Xc(a P° ,b) ~ ^ Xab ,c 

= -ah~^ (Rab- ^Rh ab ) +h~^ (D a D b a- ^D c D c ah ab ) 

-OL (X Xab ~ 2 h cd Xac Xbd)- 

To express the equation in terms of the conformal quantities we use the 
transformation law of the Ricci tensor under conformal rescalings and the 
expression of the conformal Ricci tensor in terms of the connection coefficients 
to get 

Rab - - R Kb = - - h cd Kb,cd + £>( 7b) - g Kb D c 7 C 

+lc d aKd h cf lf e b + 2jd e ch df k{aJb) c f + lKb (7c c dl d -h ef % c d j c d f ) 

+ 4 (D a log h D b log h - I Kb h cd D c log h D d log h). 
36 3 

Thus we obtain 



h*-{d t Xab ~ Xab,c P c } = - h cd kb,cd ~ Dr a %) + \ Kb D c % (4.9) 
a 2 3 

+D a D b q - - Kb D c D c q + D a q D b q - i h ab D c qD c q 
+ (2 cr - -) (?7 (q D b) q - - h ab rf D c q) + {a 2 - - a - — ) (r] a r) b - - Kb Vc rf) 

~lc d aKd K f lf € b~2% e c h d l h e{ alb) C f-\kb {lc C eT-h Cd % e f% f d) 

2 i 

+ 2 Xc(a S° fi) + - Xab S C , c - Jl* (X Xab ~ 2 h cd Xac Xbd)- 

In all the equations expressions like Dajh etc. are defined by the expressions 



which would hold if % denoted a tensor field. Finally, we get from (4.4) that 
dtkbc-kbc,dP d = 2aDcXab (4.10) 
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+2 Xab a(ai] c + D c q) + 4 a j c d (a Xb)d + Kbd P d ,c 



_ 2 ~ 

+2/3 (,a hb)dc + % h d ( b [3 , a )c + g (^ahc /3 A + h a b ft t dc)- 



If c> 0, the system Q, (Q, Q, Q, Q, ( |4~To| ) is symmetric 

hyperbolic. This can be seen as follows, (i) We choose real numbers e, / 
satisfying 



e>0, />0, / a = -2(i+|-<7) 



e=?/-(~«r)(i+|-«7), (4.11) 



2 3 



and multiply some of the equations by overall factors c, e, /, a , to obtain 
them in the form (writing out only their principal parts here) 



dt h a b — h a b^ c f3 c — 



(4.12) 



(4.13) 



{dtVa ~Va,cP c } = e {2 D aX +■■■}, 



(4.14) 



/»* -{d tX - X.c P c } = f{-aD a % + 2^D a Va + ...} 7 (4.15) 
a f 



- {d t la - la,c P c } = -cD c x ca ~ fa D aX 
a 



(4.16) 



1 c 1 - 1 ~ 

h* - {d t x ab - Xab.c P c } = c{-h cd h ab , cd - D (a % } + -h ab D c % + ...}, (4.17) 
a 2 v ' 3 



c 1 ~ 

{9 t /lafcc - h a bc,dP d } = c{-D c Xab + ■■ ■}• 

4a 2 



(4.18) 



(ii) We contract both sides of (|4.14| ) and fl4.16|) with /i ba , both sides of (|4.17|) 
with /i a ( c ft, d ) b , both sides of (|4.18| ) with h^/i e ) b fc/ c , and add on the right 
hand side of the equation obtained in this way from ( 4.16 ) a term of the 



form f h ba h cd D a 

Xcdi which vanishes identically but whose formal occurence 

makes the symmetry manifest. 



The range in which the coefficients a, c, e, /, a are consistent with (4.11) 
can be seen by considering the following cases: 



i) a = «-» a = 



2 3' 
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In this case we have e=(j + §) / and we can choose c, / > arbitrarily. If 
a^Owe have / = -~ (| + § - a) and e = {a 2=1 - 1} (1 + | - <r), which 
gives the following restrictions 

lal 1 c 

" } a< °' 2^TT) <ff< 2 + 3' 

1 c 

wi) < a < 1, o + q <f7 ' 

3 + 2c 1 c 
w ) K a < ^7—. o + o < ° < 



2c ' 2 3 2 (a - 1)' 

We cannot have e > 0, f > with a > 

It is of interest to study the characteristics of the system. Equations (4.12) 
and (4.13) contribute factors n(£) — £ m to the characteristic polynomial. 
To find the other characteristics wc analyse for which £ M ^ the following 
linear system of equations, defined by the principal symbol map, admits non- 
trivial solutions. We set £ a = h ab f> 

n (0la = -c£, b Xba ~ fa£ a X, 
h* Tl(£) Xab = ~ Kbc C C - £,(a 76) + 2 ^ 7c 
»(0 ^abc = 2£cXafc- 

The condition = implies Xab = and x — but there remains a 
fifteen-parameter freedom to choose the remaining unknowns. If 

n(0^0, (4.19) 

we derive from the equations above the further equations (writing £) = 

1 r 



^ 3 Xafc = n(f ) (£ (a 76) ~ g h ab C %), 



whence 



{<?(£, + c ~h cd L U 7a = - £h~ H c % + / a x) L ■ 

2 b n(£) 

The latter equation implies 

{<?(£, + y ^ Cc £4 C > + ^|y^ r £c / a X = 0, (4.20) 
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and, with p a such that p a £ a = 0, 

{5(^) + ^ cd £cMp a 7a = 0, 

Finally we get 

{9{t, + (a - 1) (1 - 2 <t) /i C!i £ c U X ~ « h-i n(0 £ c 7c = 0. (4.21) 



Equations (4.20) and (4.21), are of the form Au = with the unknown it the 



transpose of (x, C 7i) and the matrix A satisfying 

det A = {g(£, + ^ (1 - a) /i cd £ c &} {<?(£, + (2 <r - 1) /i c<i £ c 

If 

S(€,0=0, (4.22) 

it follows from the first of these equations that % = 0. J£ a 1, a ^= ^ 
there remains the two-parameter freedom to choose Xab with £ a Xab — 0. If 
a = 1 or u = | there remains the three-parameter freedom to choose %, Xab 

satisfying £ a Xab — — ^ £& X- If fl^i 7^ the field x»j is known once % has 
been determined. If 

gfot) + Zh ed t o t d = 0, (4.23) 



there is a two-parameter freedom to choose p 1 % as above. This is the only 
freedom unless a=^orer=i + |, conditions which exclude each other 
Further characteristics are given by the equations 



g{t,S) + Y(. 1 - a ) hCd tctd = o, (4-24) 

+ (2 a -1)^^ = 0. (4.25) 



The 'physical' characteristics correspond of course to ( 4.22 ), the two-parameter 
freedom pointed out in the case a ^ 1, a ^ h corresponding to the two po- 
larization states of gravitational waves and the additional freedom in the 
other cases corresponding essentially to a gauge freedom. The timelike char- 



acteristics corresponding to ( f4.19[ ) occur because of the transition from the 



system of second order to a system of first order. These characteristics are 



neither 'physical' nor 'harmful'. Characteristics corresponding to (4.23) are 
spacelike, while the nature of the characteristics corresponding to (4.24) and 
( 4.25| ) depends on the constants a and a. It can happen that one is spacelike 



while the other is timelike or both are spacelike. 

We note here that it is possible to obtain by similar procedures reduced 
equations, for somewhat more complicated unknowns, which have only char- 
acteristics which are timelike or null H] . 

Though the number c > can be chosen arbitrarily (suitably adjusting 
the others) it is not possible to perform the limit c — > while keeping the 
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symmetric hyperbolicity of the system, however, the equations in fil| ] can 
be considered as limit of equations which are algebraically equivalent to our 
systems. We note also that it is not possible to perform a regular limit a — > 
which would make lapse and shift the gauge source functions. In the limit 
as c — > (which we can perform if we do not insist on hyperbolicity) all 
characteristics, with the exception of the gauge dependent characteristics 



corresponding to (4.25), become non-spacelike. 

We fin ally remark on our gauge conditions. Under certain assumptions 
the gauge (^lj) coincid es w ith the gauge with harmonic time coordinate and 



prescribed shift. From (2/7) follows the general relation 



d t a-a, a t3 a = a 2 X -a 3 r°. 



Equation ((15J), written as an equations for h, entails together with (4.1) the 
equation 

d t a-a, a /3 a = 2aa 2 X + a(2ap\ a + d t q-q, a p a ). 

Thus the time coordinate t is harmonic in our gauge, i.e. r° = 0, if a = ^ 
and dtq — q, a P a = ~P a In more general situations the expression for r° 
implied by the equations above contains information about the solution and 
admits no direct conclusion about time-harmonicity in terms of a and f3 a . 

It follows from complete reductions based on the gauge source function 
q (cf. e.g |Q, 0]) that on solutions of Einstein's equations it is possible 
to achieve the corresponding gauge close to some initial surface. However, it 
has apparently never been shown that for prescribed a > 0, [3 a and given 
spacelike hypersurface of some arbitrary Lorentz manifold, coordinates can 
be constructed which realize these gauge source functions. It would be use- 
ful to have a proof of the universal applicability of this type of gauge and 
information about its general behaviour. 



4.2 The Einstein Euler System 

In the following we shall discuss the Einstein-Euler equations, i.e. Einstein's 
equation coupled to the Euler equation for a simple perfect fluid. Its hyper- 
bolicity has been studied by various authors (cf. @, [||, @). Though the 
system is also important in the cosmological context, our main concern in 
analysing the system here is to control the evolution of compact perfect fluid 
bodies, which are considered as models for 'gaseous stars'. In this situation 
arises, besides the need to cast the equations into hyperbolic form, the side 
condition to control the evolution of the timelike boundary along which the 
Einstein-Euler equations go over into the Einstein vacuum field equations. 

The analysis of the evolution of the fields in the neighbourhood of this 
boundary poses the basic difficulty in the discussion of compact fluid bodies. 



In the case of spherical symmetry, this problem was overcome in 61 . If 



in more general situations the coordinates used in the hyperbolic reduction 
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are governed e.g. by wave equations, as is the case in the harmonic gauge, 
there appears to be no way to control the motion of this boundary in these 
coordinates. In J39| a system of equations has been derived from the Einstein- 
Eulcr system which combines hyperbolicity with the Lagrangian description 
of the flow, so that the spatial coordinates are constant along the flow lines. 
Thus the location of the body is known in these coordinates. In the following 
we shall discuss this system and derive the subsidiary system, which was not 
given in |3^| . 



The Basic Equations We shall use a frame formalism in which the in- 
formation on the metric is expressed in terms of an orthonormal frame 
{efc}fc = o,...,3 and all fields, with the possible exception of the frame itself, 
are given in this frame. To make the formalism easily comparable with the 
spin frame formalism which will be used later, we shall use a signature such 
that gik = g(ei,ek) = diag(l, —1,-1,-1). Let V denote the the Levi-Civita 
connection of <? M „. The basic unknowns of our representation of the Einstein 
equations are given by 

e M fc , r t 3 fc , C l jki , matter variables, 

where e M k — < &k , x ^ > are the coefficients of the frame field in some coor- 
dinates x^, Ti 3 k are the connection coefficients, defined by V, eu = 3 k &j 
and satisfying Jj 3 k gji + Ti 3 1 gj k — 0, and C l j k i denotes the conformal Weyl 
tensor in the frame e k . The latter is obtained from the decomposition 

Rijki = Qjki + {g, t [ k S^j - g^ k Sqi}, (4.26) 

of the curvature tensor 

^ 3P1 = e p(-Tq j) ~ e q{-Tp j) — r k j (r p q — T q p ) (4.27) 

I T~t i pfc r-i i pfc 

p k ! q j * q fc ^ p ji 

where we also set Sij — Rij — g gij R, with Rij and R denoting the Ricci 
tensor and the Ricci scalar. We shall need the notation 

Ti j Cfc [ci, 6j] -\- (I^i j Rj i) 6/ , 

A 1 jki = R l jki - C % jki - g l ^ Si]j + g 3 [k Sq \ 
with R 1 jki understood as being given by ( 4.27 ). Furthermore we set 

Fjki = Vi F l jki, with F l jki — C % jki — g % \k Sqj. 

In the equations above we take account of the Einstein-Euler equations 
in the form 

S ik = K{T ik -\g jk T), (4.28) 
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with an energy-momentum tensor of a simple perfect fluid 

T lk = {p + p)U l U k -pg, k . (4.29) 

Here p is the total energy density and p the pressure, as measured by an 
observer moving with the fluid, and U denotes the (future directed) flow 
vector field, which satisfies U% U l = 1. 
We shall need the decomposition 

V J ' T jk = q k + qU k , (4.30) 

Jjk = Vyg fc ], (4-31) 



and the field 
with 



q = U l V lP +{p + p)V l U\ (4.32) 

q k = {p + p) U l V, U k + {U k U* V, - V k }p. (4.33) 

We assume that the fluid is simple, i.e. it consists of only one class of par- 
ticles, and denote by n, s, T the number density of particles, the entropy per 
particle, and the absolute temperature as measured by an observer moving 
with the fluid. We shall assume the first law of equilibrium thermodynamics 
which has the familiar form d e = —pdv + Tds in terms of the volume v = — 
and the energy e = £ per particle. In terms of the variables above, we have 

d p = P-LZdn + nTds. (4.34) 

n 

We assume an equation of state given in the form 

p = f(n,s), (4.35) 

with some suitable non- negative function / of the number density of particles 
and the entropy per particle. Using this in (4.34), we obtain 

V = n- p, T=-—, (4.36) 

on n a s 

as well as the speed of sound v, given by 

v = ^~) s = — r~ ( 4 - 37 

op p+P on 

as known functions of n and s. We require that the specific enthalpy and the 
speed of sound are positive, i. e. 



(4.38) 
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We assume the law of particle conservation 

U i V l n + nV l U l = 0. 



(4.39) 



It implies together with q = and ( 4.34 ) that the flow is adiabatic, i.e. the 
entropy per particle is conserved along the flow lines, 



U l Vi s = 0. 



(4.40) 



The case of an isentropic flow, where the entropy is constant in space and 
time, is of some interest. In this case the equation of state can be given in 
the form 

p = h(p) (4.41) 

with some suitable function h. As a special subcase we shall consider pressure 
free matter ('dust' ), where h = 0. 

We note that |L3f ), j~36|) , ( p^ ), flOg ), and ( p0|) imply 

TP V lP = -( P + p) v 2 V t U\ 

from which we get 

q k = -V fc p + (p+p) {U l V, U k - v 2 U k Vi U 1 }. 
Finally, ( 4.40| ) implies the equation 

C-u Sk = 0, 



(4.42) 



(4.43) 

for Sfe = Vfc s, where Cjj denotes the Lie derivative in the direction of U. 

The Einstein-Eulcr equations are given in our representation by the equa- 
tions z = 0, q = 0, (HI), (HI), and ( ^4C| ), where we denote 
by z the vector-valued quantity 



z = (Ti j, A 1 jki, Fjki, qk)- 



(4.44) 



which we shall refer to (as well as to each of its components) as a 'zero quan- 
tity'. These equations entail furthermore Jjk = and (4.43). After making 
a suitable choice of gauge conditions we shall extract hyperbolic evolution 
equations from this highly overdetermined system . 



Decomposition of Unknowns and Equations From now on we shall 
assume that 

e = U, 

so that we have U l = 8 l q. For the further discussion of the equations we 
decompose the unknowns and the equations. With the vector field U we 
associate 'spatial' tensor fields, i.e. tensor fields Ti^...^ satisfying 



i' o. / i. ••-./;. 
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The subspaces orthogonal to U inherit the metric hij = gij — Ui Uj , and hi 3 
(indices being raised and lowered with <jty) is the orthogonal projector onto 
these subspaces. 

We shall have to consider the projections of various tensor fields with 
respect to U and its orthogonal subspaces. For a given tensor any contraction 
with U will be denoted by replacing the corresponding index by U and the 
projection with respect to hi J will be indicated by a prime, so that for a 
tensor field we write e.g. 

^iUk ~ Trnpq hi U P hk 9 , 

etc. Denoting by eyjy the totally antisymmetric tensor field with 60123 = 1 
and setting ejki = e'ujkU we nave the decomposition 

eijki = 2 {U[i€j] k i - eij[k Uq), 

and the relations 

^ e jpq — ~ ^ e h k [ p h g ], e? ej kq — ~2eh l q . 

Denoting by C*j k[ — \ Cij pq eki pq the dual of the conformal Weyl tensor, 
its {/-electric and U- magnetic parts are given by by Eji = Cy j[/; , Bji = C\ s ^ lJl 
respectively. With the notation Ijk = hjk — Uj Uk, we get the decompositions 

Cijkt = 2 (lj[k Eqi - h[k Eqj) - 2 (E7j- fc Bq p e p i3 + B^ p e p M ), (4.45) 

C* jkl = 2 U[i E j]p e p kl -4 E p[t e j} p [k Uq - 4 U [t B j} [k Uq - B pq e p y e 9 w . (4.46) 
We set 

a l = U k V k U\ Xij = h i k V k U j , X =h ij Xij , 
so that we have 

V j U i = U j a i + X j i , a* = 1^*^*0, Xa = hj l r k ° ,. 

Since U is not required to be hypersurface orthogonal, the field Xij will in 
general not be symmetric. If the tensor field T is spatial, i.e. T — T", we 
define its spatial covariant derivative by T>T — (VT)'. i.e. 

-^ii,.. . ,i p — ^ j -^Ji ,. . . ,j p hi 3 hi 1 3 ... hi p ^ p . 

It follows that 

T>i hjk = 0, T>i e jk i = 0. 



Equation (4.40) then implies Sk = T>k s. 

To decompose the equations, we observe the relations 

q = - h? F> U} , q k = -~ (h» F[ jk - Flj ku ) , 
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and set 

Pi = Fmu 

= V j {Eji - « (I p + I p ) /i^} + I « (p + p ) a, + 2 x fcl e j , (i S fc)j , 

r> — — ^ t ki -pi 
= V k B kj + ej kl (2 X \~ X k ') + k (p + p) X ki ej k \ 

P H = i F (i\u\j) ~ 3 h v hkl F Lui} = E i:j + T> k tj) kl 
-2 a k e kl (i B j)t - 3 X (i 1 E j)t -2 X l (l E j)t + h tJ X kl E kl +2 X E i5 

K 1 

+ 2 (P + P) (X(ij) ~ ^Xhij), 

Qki = 2 e (fe lj F[)ij = £u B k i - T>i E i(k e Z ) %i 

+2 a, e ij (k E t)j - X l (k B l)t -2 X{k l B X)l + X B kl - Xij B pq e pl (k e jq . 
Then we find the splitting 

Ejkl = 2 Uj P [k Uq + h j[k Pq + Q 4 (Uj e* kl - e l j[k Uq ) (4.47) 

-2Pj[k Ui] - Q 0l e l a - -Khj[k qi] - -nqhj[ k Uq. 

The Reduced System i) In the case of pressure free matter the equation 
q k = tells us that e is geodesic and we can assume the frame field to 
be parallelly transported in the direction of eo- Furthermore, we can assume 
the coordinates x Q , a — 1,2,3 to be constant on the flow lines of e and 
the coordinate t — x° to be a parameter on the integral curves of eo. These 
conditions are equivalent to 

r j k = o, e% = = 

which, in turn, imply together with the requirement p = the relation q k = 
if p + p > 0. The unknowns to be determined are given by 

^ ( e ^a, Pb j : Eij , Bij , p), 

where a, b = 1,2,3. The reduced system is given by 

T i a =0, 4 0a = 0, P ij= 0, Q H = 0, q = 0. (4.48) 
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ii) In the general case we cannot assume the frame to be parallelly trans- 
ported because the evolution of U is governed by the Euler equations. We 
assume the vector fields e a , a = 1,2, 3, to be Fermi transported in the direc- 
tion of U and the coordinates to be chosen as before such that 



r Q = o, 



= U» = 5» i 



Our unknowns are given by 

u=(e^ ai r ° a , r a k h Eij, B k i, p, n, s, s k ), 
and the reduced system by 



bOa 



0. 



V 2 A C 0ac - 



v 2 (A° aQb + 



J'ab) 



= o, 



Pij = 0, Qki = 0, q = 0, 
Cu<n=-nx, C v s = Q, C v s k 



(4.49) 



(4.50) 



0, 



where it is assumed that the functions p, u 2 , a, /3 are determined from the 
equation of state ( 4.35 ) according to ( 4.36 ), ( 4.37 ), etc. and p + p > 0. 

Since the functions p, s, n are then determined as functions of the coordi- 
nates x^, we remark that the relation p(x^) = f(n(x^),s(x^)) is satisfied as 
a consequence of the equations for p, s, and n, and the relation (Of]), since 
q = implies dt p = -ji f{n, s). Furthermore, we note that in our formalism 
Cu Sk = d t s k - (r 3 k - r k 3 ) sj = k d t s^ etc. 



When we solve the equations Py = and Qu = 0, the symmetry of 
the fields E^, B k i has to be made explicit. The trace- free condition then 
follows as a consequence of the equations and the fact that the initial data 



are trace- free. With this understanding it is easy to see that the system (4.48) 
is symmetric hyperbolic, the remaining equations only conta ining deriv atives 
in the direction of U. To see that the system consisting of ( 4.49| ) and ( 4.50 ) 



is also symmetric hyperbolic, we write out some of the equations explicitly 
(taking the opportunity to correct some misprints in p9| ). It follows directly 
from the definition that 

J kj = {p + p) {U l (V fc V, Uj - Vj V, U k ) - v 2 Uj V fc V,; U l (4.51) 
W U k Vj V, U l - v 2 V; U l (Vfc Uj - Vj U k ) + Vfc U l V 4 Uj - Vj IP v* u k 
+e p -^r {q^)> v * u 1 (P h v* V, Uj - Uj W V, u k ) 
+(a U k V; U 1 ~PU 1 V* U k ) Sj - {a Uj V, U* -[3U l V l Uj) s k , 
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where we set 

. , , , 

9s v 2 dn d s v 2 d s 



, , ,du 2 n dv 2 dp 2 I dp 



In particular, 
1 



= e (r a ) - v 2 e a (r c c ) + r c (r a c -r Q c a ) (4.52) 
p + p 

P + P r d 2 p^ __ 2 \ c 



and 

1 



, (tt) s — V M c C ^0 a j ^c C s a, 

IS Z O p z ) P+P 

j' ab = e a {r \)~ e b {r Q ° a )~r ° c {r a c b ~r b c a ) (4.53) 



-^r c %(r a \-A^}-^(r ^ 6 -r \ Sa ). 

p + p 



From this follows that the last two equations of ( 4.49Q are given by 



dt r a - v l e c (r a c o) = -To c r a c (4.54) 
p + p / d 2 p^ 



, (TT) S _ ^ r c C 0-To aH ; ^c C 0S a 

v l a p z J P + P 

-v 2 (ifc c o (r a k c — r c k a ) — r a c k r c k o + r c c k r a k o + r c 0ac ) , 
v 1 d t r a \- v 1 e b (r a ) = v 2 (r k b (r k a - r a k ) (4.55) 

— -H) ° c -Ta C 6 + -R° 60a + -H) ° c (-T a C & — A C a ) 

+i , 2 r c %(r a %-A\) + ^(r \ S6 -r \ Sa ) > ) . 

p + p / 



The remaining equations of (4.49) and (4.50) (besides Py = 0, Q^; = 0) 



again only contain derivatives in the direction of U. 

The Derivation of the Subsidiary Equations We have to show that any 
solution to the reduced equations which satisfies the constraints on an initial 
hypersurface, i.e. for which z = on the initial hypersurface, will satisfy z = 
in the domain of dependence of the initial hypersurface with respect to the 
metric supplied by the solution. For this purpose we will derive a system of 
partial differential equations for those components of z which do not vanish 
already because of the reduced equations and the gauge conditions. 

We begin by deriving equations for Fju . There exist two different expres- 
sions for Fki = V J Fjki ■ From the definition of Fju and from the symmetries 
of the tensor field involved follows 
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F H = V J V 1 F ijkl = V J V 1 Cyki - V J V [fc % 

= -R p [i y Cj] P kl + Cij P [i BP k] 13 + 2 Ti p j ■ V p C n ki 

+V [fc V' S l]3 - RJ> [j » fc] 5 pj - R p j 3 {k S l]p + V p p T fc] p „ 

where we took into account that we do not know at this stage whether the 
connection coefficients Fi 3 k supplied by the solution define a torsion free 
connection. From the reduced field equations, the definition of the zero quan- 
tities, and the symmetries of Cijki , it follows that 

F M = -A? {l « C j]pkl + C ijp[l A? k] 13 + \T i p j V p C* ki (4.56) 

+k J k i + A p p 1 k ] S pj - A p j 1 [k S l]p = N(z), 

where N(z) is, as in the following, a generic symbol for a smooth function 
(which may change from equation to equation) of the zero quantities which 
satisfies N(0) = 0. 

On the other hand, because of the reduced equations, equation (4.47) 
takes the form 

F jkt =2Uj P[ k Ui] + hj[ k Pi] + Q t (Uj e l k i - e z j[k Uq)-~K h j[k q t] . (4.57) 

Contracting with V- 5 , decomposing the resulting expression i nto F Lq , F^ and 
equating with the corresponding expressions obtained from fl4.56| ) , we arrive 
at equations of the form 

Cu Pk + \ e k 13 V, Qj = N(z) (4.58) 

CuQk-^ek ii 'DiPj = N{z). (4.59) 

The connection defined by the i~j 3 k and the associated torsion and cur- 
vature tensors satisfy the first Bianchi identity 

(jkl) (jkl) 

where Yl(jkl) denotes the sum over the cyclic permutation of the indices 
jkl. Setting here j = 0, observing that the symmetries of C l j k i, Ski imply 
i kl = S(j/ci) A 1 jki, and taking into account the reduced equations, 
we get from this an equation of the form 

e Q {T a \) =N{z). (4.60) 
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To derive equations for A % jki, we use the second Bianchi identity 



y ' R % mkl — — y ' R % mnj Tk " , 

(Jkl) (jkl) 



(4.61) 



We write 

-R* jfez = A % jki + C jki + E 1 jki + G l jki , 

with E l -jki = g i [k S^j-g^k S£ 4 , G* j fc ; = ± 5 ,g l [fe and S* k = Sjfc S ftfe. 
Using the well known facts that the left and right duals of Cijki and Gijki 
are equal while the left dual of Eijki differs from its right dual by a sign, and 
using the reduced equations, we get 



Vj i?' m fe; — Vj A 1 

{jkl) (jkl) 



mkl 



2 ^jkl q prs ^ ' q p rs -f~ prs) 



^ Vj A 1 mk l + - e^fej P (-Fprs + ^#pr Qs) e 



irs 
m i 



(Jkl) 



whence, by ( [4.57D and ( }4 . 6 1| ) 



2V(*)- 



(4.62) 



i) In the case of pressure free matter we get from ( 4.62] ), by setting j — and 
using ( ff.48 ), an equation of the form 



e Q (A i k ab)=N(z). 



(4.63) 



Since, as remarked earlier, qk = by our gauge conditions, the system of 
equations consisting of ( 4.58| ), ( [f.59] ), ( 4.60 ), and ( 4.63 ) constitutes the de- 
sired 'subsidiary system' for the zero quantities in the pressure free case. 

ii) Using ( 4.49 ), we get in the general case by the analogous procedure 
only an equation of the form 

e (A a icd ) = N(z). (4.64) 
By the antisymmetry of J' ab we know already that 

A° a ob + A° bOa = 0. 

We can express the equations for A [b\o\c], A aoc . in terms of the quantities 



A„ 



be aO 



A* — 

^ab — 



Ha 



cd a0 



b)cd ■> 



A* a = -A c 0ac , (4.65) 
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because 

^° [6|o|c] = e a bc , A abc = -A* ad e bc d + 2 h a [ b A*y 
From ( [4.62 ) we get 

e {A abc) + e c (A a0b ) + e b (A ac0 ) = N{z), 
which implies for the quantities ( |4.65| ) equations 

2e (A:)-e a bc e b (A c ) = N(z), 

and 

e (A* ab ) - e (Q (4 6) ) + h ab h cd e c (A d ) = N{z). 

Since we have 

h ab e a (A b ) = -\ e abc e a (-L- J' hc ) = -\-^~ e ijk V, Vj q k +N(z) = N{z), 
2 P + P 2 p + p 

we can write the second equation in the form 

e (A* ab )-e {a (A b) )=N(z). 



From (4.62) we get furthermore 

e a (A bcd ) + e d (A bac ) + e c (A bda ) = N(z), 
which implies in terms of the quantities ( |4.65| ) an equation of the form 

h ab e a (A* bc ) - e c ab e a (A* b ). 
By a direct calculation we derive from ( |4.31[ ) the equation 

2 Cu J' i} = AV {1 J' Vj] - hi p hj q {A 1 npq + A 1 qnp + A 1 pqn ) U n q t 

-h,P hj* (Ti n p V n q q + T q n ,V„fc +T p n q V n qi)U l 
—2 a,i J'uj + 2 aj J'm, 
which can be rewritten by (4.49) in the form 

(p + v) {eo (A») + 2v 2 e a bc e b (A* c )} = N{z). (4.66) 

From the equations above we obtain the system 

(p + p) {2 v 2 e {AD - v 2 e a bc e b (A c )} = N(z), (4.67) 

(p + p) {e (A,) + v 2 e a bc e b {A* c ) + v 2 h bc e b {A* ca )} = N{z), (4.68) 
cab ( P + P) W 2 eo (A* ab ) - v 2 e (a {A b) )} = N{z), (4.69) 
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where c a b = 1 if a = 6 and c a b — 2 if a ^ b. Finally, we obtain from (4.31), 
(4.4£) the equation 

£uq a = 2v 2 (p + p)A c 0ac . (4.70) 

Equations ( gj| ), ( gj| ), ( ggg ), ( ggg ), ( gjgg ), ( gggj ), ( ggg ), and ( |4~7C| ) con- 
stitute the subsidiary equations for the zero quantities in the general case. 

We note here that in the reduced system it has not been built in explic- 
itly that the energy-momentum tensor has vanishing divergence (cf. equation 
( 4.3C| )). While we assumed the equation q — as part of the reduced equa- 
tions, we verify the vanishing of the quantity q^ by deriving the subsidiary 
equations and using the uniqueness property for these equations. 

In the present formalism the gauge conditions are taken care of by the 
explicit form of some of the unknowns, however the list of constraints is 
much longer than in the previous discussions. We shall not try to demon- 
strate that the constraints are preserved in the specific case of 'floating fluid 
balls'.. Though the construction of data for fluid balls of compact support 
which are embedded in asymptotically vacuum data has been shown and 
their smoothness properties near the boundary have been discussed j70| , the 
evolution in time of these data and the precise smoothness properties of the 
fields near and possible jumps travelling along the boundary have not been 
worked out yet. However, without a precise understanding of the behaviour 
of the solution near the boundary the conservation of the constraints cannot 
be demonstrated. 

Our reduced system also found applications in cosmological context where 
the fluid is spread out, with p + p > 0, over the time slices (cf. Q). 
In this case the desired conclusion follows from the fact that the subsidiary 
systems are symmetric hyperbolic, have right hand sides of the form N(z), 
and the characteristics of the reduced system and the subsidiary system are as 
follows (where we use only the frame components — Cm e M k of the covector 

0- 

(i) In the case of pressure free matter the characteristic polynomial of the 
reduced system is of the form 

o (Zof (Co 2 + \h ab Ca Hb) L (<T Cm & N , 

with positive integers K, L, N and constant factor c, while the characteristic 
polynomial of the subsidiary system only contains the first two factors. Thus 
the characteristics of the subsidiary system are timelike with respect to 
(cf. also the remarks in Sect. |^). 

(ii) In the general case the characteristic polynomial of the reduced equa- 
tions is of the form 

C (£ ) K (Co' + \ £a 6)^ (Co' + ^ h Cd Cc Crf) M (ST Cm Zu)", (4-71) 

with positive integers K, L, M, N and constant factor c, while the character- 
istic polynomial of the subsidiary system is generated by powers of the first 
three factors. 
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We note that the equations (4.54) and (4.55) contribute to the principal 
symbol the third factor which corresponds to the sound cone pertaining to the 
fluid. If one wants to ensure that the sound does not travel with superluminal 
speed one has to require the equation of state (4.35) to be such that v < 1. 
For the question of existence and uniqueness of solutions there is no need 
to impose such a condition, but if v > 1, the domain of dependence with 
respect to <? M „, as it has been defined in Sect. [^, cannot be shown any longer, 
by the arguments given in Sect. ^, to be also a domain of uniqueness. In any 
case, it follows from the characteristic polynomials that in a domain where 
the solution of the reduced system is unique according to those arguments, 
the constraints will be satisfied if they hold on the initial hypersurface. 

We note that the system simplifies considerably in the isentropic case. In 
the reduced system the function p+p then neither occurs in the principal part 
nor in a denominator. In the case of the subsidiary system a more detailed 
discussion is required to understand the consequences of the occurrence of 
the various factors p + p. 

In our procedure the fluid equations serve two purposes, they determine 
the motion of the fluid as well as the evolution of the frame. If we set n = 
in all equations the fluid equations decouple from the geometric equations 
and we obtain a new hyperbolic reduction of the vacuum field equations. In 
this procedure any exotic 'equation of state' may be prescribed as long as it 
ensures a useful, long-lived gauge. 

If the initial data for the Einstein-Euler equations are such that U, p, p, 
n, s, and the equation of state can be smoothly extended through the bound- 
aries of the fluid balls, this suggests using the 'extended fluid' to control the 
evolution of the gauge in the vacuum part of the solution near the boundary. 



4.3 The Initial Boundary Value Problem 

In the previous section we studied a problem involving a distinguished time- 
like hypersurface. Its evolution in time was determined by a physical pro- 
cess. There are also important problems where the Einstein equations are 
solved near timelike hypersurfaces which are prescribed for practical reasons, 
e.g. to perform numerical calculations on finite grids. The underlying initial 
boundary value problem for Einstein's field equations, where initial data are 
prescribed on a (spacelike) hypersurface S and boundary data a (timelike) 
boundary T which intersect at a 2-surface S = T S, has been analysed in 
detail in the article p2[ . The solution to this problem requires a hyperbolic 
reduction which needs to satisfy, beyond the conditions discussed at the be- 
ginning of this section, certain side conditions. In the following we want to 
comment on those aspects of the work in |Q which illustrate the flexibility 
of the field equations in performing reductions and on certain characteris- 
tics of the reduced system. For the full analysis of the initial boundary value 
problem we refer to p^ ]. 
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Since we are dealing with a problem for equations which are essentially 
hyperbolic, the problem can be localized. In suitably adapted coordinates x^, 
defined on some neighbourhood of a point p £ E, the manifold M on which 
the solution is to be determined will then be given in the form M = {x € 
R 4 |:r > 0, x 3 > 0}, the initial hypcrsurface by S = {x £ M\x° = 0} and the 
boundary by T = {x £ M\x 3 = 0}. Clearly, we will have to prescribe Cauchy 
data on S as before, we will have to analyse which kind of boundary data are 
admitted by the equations, and on the edge E = {x £ M\x° = 0, x 3 = 0} 
the data will have to satisfy some consistency conditions, as is always the 
case in initial boundary value problems. 



Maximally Dissipative Initial Boundary Value Problems We have 



seen in Sect. 3.1 that energy estimates provide a basic tool for obtaining 
results about the existence and uniqueness of solutions to symmetric hyper- 
bolic systems. To explain the side conditions which have to be satisfied in a 
hyperbolic reduction of an initial boundary value problem for Einstein's field 
equations, we consider what will happen if we try to obtain energy estimates 
in the present situation. Assume that we are given on M in the coordinates 
x 1 * a linear symmetric hyperbolic system of the form 

A*d„u = Bu + f(x), (4.72) 

for an R^-valued unknown u. The matrices — A^(x), /i = 0,1,2,3, 
are smooth functions on M which take values in the set of symmetric N x 
A^-matrices, there exists a 1-form £ M such that A^ £ M is positive definite, 
B = B(x) is a smooth matrix- valued function and f(x) a smooth R^-valued 
function on M. For convenience we assume that the positivity condition is 
satisfied with £ M = 5° p. 

If we assume that u vanishes for large positive values of x a , a — 1, 2, 3, 
and if the relation 

d^uA^u) = t uKu + 2 t uf with K = B + * B + <9 M A^ , 



implied by (4.72), is integrated over a set M T = {x £ M\0 < x < r}, defined 



by some number r > 0, we obtain the relation 

/ t uA°udS^ [ t uA°udS+ [ { t uKu + 2 t uf}dV + [ t uA 3 udS, 
Js T Js Jm t Jt t 

involving boundary integrals over S T — {x £ M\x° = r} and T T — {x £ 
M\0 < x° < r, x 3 — 0}. Obviously, the structure of the normal matrix A 3 
plays a prominent role here. If the last term on the right hand side is non- 
positive, we can use the equation above to obtain energy estimates for proving 
the existence and uniqueness of solutions. 

By this (and certain considerations which will become clear when we have 
set up our reduced system) we are led to consider the following maximally 
dissipative boundary value problem. 
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We choose g S C°°(S, TL N ) and require as initial condition u(x) — g(x) 
for x £ S. We choose a smooth map Q of T into the set of linear subspaces 
of R N and require as boundary condition u(x) S Q(x) for x £ T. The type 
of map Q admitted here is restricted by the following assumptions. 



(i) The set T is a characteristic of (4.72) of constant multiplicity, i.e. 

dim(kerA 3 (x)) = const. > 0, x £ T. 

(ii) The map Q is chosen such as to ensure the desired non-positivity 

t uA 3 (x)u<0, u£Q(x), x£T. 

(iii) The subspace Q{x), x £ T, is a maximal with (ii), i.e. the dimen- 
sion of Q{x) is equal to number of non-positive eigenvalues of A 3 counting 
multiplicity. 

The specification of Q can be expressed in terms of linear equations. 
Since A 3 is symmetric, we can assume, possibly after a transformation of the 
dependent variable, that at a given point x £ T 



-Ij 
fc 
h 



K > 0, 



where Ij is a j x j unit matrix, Ofe is a k x k zero matrix etc. and j + k + 1 = N . 
Writing u = t (a, b, c) £ R J x R fc x R ; we find that at x the linear subspaces 
admitted as values of Q are neccessarily given by equations of the form = 
c — H a where H = H(x) is a / x j matrix satisfying 

- t aa+ t a t HHa<Q, a£R 3 , i.e. t HH<I j , 

We note that there is no freedom to prescribe data for the component b of u 
associated with the kernel of A 3 . More specifically, if A 3 = on T, energy es- 
timates are obtained without imposing conditions on T and the solutions are 
determined uniquely by the initial condition on S. By subtracting a suitable 
smooth function from u and redefining the function /, we can convert the 
homogeneous problem above to an inhomogeneous problem and vice versa. 
Inhomogeneous maximal dissipative boundary conditions are of the form 

q = c-Ha, (4.73) 

with q = q(x), x £ T, a given R'-valued function representing the free bound- 
ary data on T. 

Maximally dissipative boundary value problems as outlined above have 
been worked out in detail in Q, Q for the linear case and in Q Q for 
quasi- linear problems (see also these articles for further references). If we want 
to make use of this theory to analyse the initial boundary value problem for 
Einstein's field equations we will have to solve two problems which go beyond 
what is known from the standard Cauchy problem. 
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(i) We will have to find a reduction, involving a symmetric hyperbolic 
system, which gives us sufficient information on the normal matrix so that 
we can control the conditions above. 

(ii) To demonstrate the preservation of the constraints we will have to 
discuss an initial boundary value problem for the subsidiary system. This 
should be such as to admit a uniqueness proof. Moreover, there is the problem 
of getting sufficient control on the solution near T. While the initial data for 
the reduced system will of course be arranged such that the constraints are 
satisfied on the initial hypersurface 5, it will a priori not be clear that we will 
have sufficient information on the behaviour of the solution to the reduced 
equations and on the data on T in order to conclude that the constraints will 
be satisfied on T. 

The choice of representation of the field equations, of the gauge conditions 
and the gauge source functions, and, in particular, the choice of the reduced 
equations will largely be dominated by the second problem. 

The Representation of the Einstein Equations In the initial bound- 
ary value problem for Einstein's vacuum field equations was analysed in terms 
of the equations 

2^ = 0, A i jkl =Q, FjM = 0, (4.74) 

of the previous section with everywhere vanishing energy-momentum tensor. 
We shall use these equations together with the conventions and notation 
introduced in the previous section. 

The Gauge Conditions The gauge, which we assume here for simplicity 
extends to all of M, has been chosen as follows. On the initial hypersurface 
x° = and x a , a = 1,2,3, are coordinates with x 3 = on £ and x 3 > 
elsewhere. The timelike unit vector field eo on M is tangent to T, orthogonal 
to the 2-surfaces S c = {x 3 = c = const. > 0} in S, and it points towards M 
on S CiU. The coordinates x^ L satisfy e M = e {x^) — <5 M on M and the sets 
T c = {x 3 = c} are smooth timelike hypersurfaces of M with T = T. The unit 
vector field is normal to the hypersurfaces T c and points towards M on T. 
The vector fields e^, A = 1, 2, are tangent to T c n S and such that they form 
with eo, e3 a smooth orthonormal frame field on S. On the hypersurfaces T c 
these fields are Fermi transported in the direction of eo with respect to the 
Levi-Civita connection D defined by the metric induced on T c . The eu form 
a smooth orthonormal frame field on U . We refer this type of gauge as an 
'adapted gauge'. Notice that it leaves a freedom to choose the timelike vector 
field e on M \ S. 

In analysing the initial boundary value problem it will be necessary to 
distinguish between interior equations on the submanifolds S, T, T c , £, S c . 
Since our frame is adapted to these submanifolds, this can be done by dis- 
tinguishing four groups of indices. They are given, together with the values 
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they take, as follows 

a, c, d, e, f = 0, 1, 2; i, j, k, I, m, n = 0, 1, 2, 3; 
p,q,r,s,t= 1,2,3; A,B,C,D = 1,2. 

We assume the summation convention for each group. 
By our conditions the frame coefficients e M & satisfy 

e f* =(5%, e 3 a = 0, e 3 3 >0 on M, (4.75) 

while the frame connection coefficients satisfy 

r A B = 0. (4.76) 

The fields e a satisfy on T c the equations 

D eo e = r A e A , D eo e A = -gAB r B Q e . (4.77) 

Thus, given the hypersurfaces T c , the evolution of the coordinates x a , a = 
0,1,2, and the frame vector fields e a off S is governed by the coefficients 
r A 

Another part of the connection coefficients defines the intrinsic connec- 
tion D on T c , since D a e c = D Ca e c — r a b c e b . The remaining connection 
coefficients, given by 

Xab = .g(V Ga e 3 ,e b ) = 3 g jb = r a 3 b = r (a 3 b) , (4.78) 

define the second fundamental form of the hypersurfaces T c in the frame e a . 
In the reduced equations, the symmetry of \ab has to be taken into account 
explicitly. A special role is played by mean extrinsic curvature 

X = g ab Xab = 9 jk r j 3 k = V ll e^, (4.79) 

since it can be regarded as the quantity controlling the evolution of the hy- 
persurfaces T c and thus of the coordinate x 3 . 

We now choose two smooth functions F A £ C°° (M) as gauge source func- 
tions. These will occur explicitly in the reduced equations and will play the 
role of connection coefficients for the solution, namely F A — Pq A o- Further- 
more we will choose a function / S C°°(M) which will play the role of the 
mean extrinsic curvature on the hypersurfaces T c . Here the interpretation is 
somewhat more complicated. On T the function \ — /It must be regarded 
as the free datum which, together with certain data on ZJ, indirectly spec- 
ifies the boundary T. However, for x 3 > the function / plays the role of 
a gauge source function which determines the gauge dependent hypersurface 
T c , c > 0. It is a remarkable feature of the Codazzi equations that they admit 
this freedom while at the same time implying hyperbolic equations. 

This example clearly shows the importance of the freedom to dispose of 
the gauge source functions. While we could choose F A = locally (cf. the 
remarks in 0] about certain subtleties arising here) , we need the full freedom 
to make use of the gauge source functions /, since otherwise we could only 
handle restricted types of boundaries. 
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The Reduced Equations Using the gauge conditions above, we extract 



from (4.74) the following reduced system for those components of the un- 
knowns e' 1 k , k, Eij , Bki which are not determined already by the gauge 
conditions and the chosen gauge source functions. Where it has not already 
been done explicitly, it is understood that in the following equations the con- 
nection coefficients fo^o an d X are replaced by the gauge source functions 
F A and / respectively. The torsion free condition gives 

= -T o k p e» k = t e%-(ro%-r p »o)e%-r o V M o- (4.80) 

The Gauss equations with respect to T c provide the equations 

= A B ooa = e (r A B )- e A (F B ) +r c B r A c (4.81) 

—T A b c F c + F b F c g AC + Xo B Xao — Xa B Xoo - C B oa, 

= A B C oa = e (r A B c) + F B r A ° c +r A B F D gen (4.82) 

+-Td B c r A D o + Xo B Xac - XA B Xoc — C B CQA- 
Codazzi's equations with respect to T c imply 

= g ab A 3 abl - D xoi - Di X u - D 2 X u ~ Di(J), (4.83) 
= g ab A 3 ab2 = D X02 - Dx X u - D 2 X22 - D 2 (f), (4.84) 
= A 3 10 i = D xu - Dx xoi - C 3 101 , (4.85) 

= A 3 am + A 3 102 - 2 Do X12 - Di X o2 ~ D 2 X oi - C 3 201 - C 3 W2 , (4.86) 

= A 3 202 = D X22 - D 2 X02 - C 3 202, (4.87) 

where it is understood that the component xoo, which appears only in un- 
differentiated form, is given by xoo = Xu + X22 + /■ The remaining Ricci 
identities give 

= a a B o3 = e (r 3 A B ) + f a r 3 B + r 3 A f c 9bc + r c A B r 3 c (4.88) 
+-^3 A b r 3 3 o + xo A r 3 3 b — r 3 A 3 X ob — rc A b xo c — c A so3, 



o = a a 003 = e (r 3 A o) - e 3 (F A ) + Xo A r 3 3 „ - r 3 A b f b + r B A o r 3 B 

(4.89) 
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+Pa A o r 3 3 o 



r A ^ 

± B XO 



r 3 

^ 3 



s 9 BA Xoo 



F A xoo-C / 



003; 



Zi 3 aos + A 3 03A = 60(^3 3 A ) - e A (r 3 3 ) 



-^3 3 o F B g BA + T 3 3 c Ta 



c 



(4.90) 



= g ab A 3 ab3 = e (r 3 3 ) + g AB e A (r 3 3 B ) - e 3 (/) 

-g ab r 3 3 fe n k a + g ab r b 3 k r 3 k a + g ab r m 3 a (r 3 m b -r b m 3 ). 



(4.91) 



In the previous section we saw how to extract a symmetric hyperbolic 
system from the Bianchi identities. However, for reasons given below, we 
shall not choose that system here. Instead we choose the 'boundary adapted 
system' 

Pu - P22 = Qix - Q22 = 

2Pi 2 =0 2Q 12 = 

Pu + P22 = Q n + Q22 = (4.92) 

Pl3 = 2Q2 Ql3 — — 2P2 

P23 = —?Qi Q23 = oPti 



2 *tfl ^2-3 — 2 

written as a system for the unknown vector u which is the transpose of 

((E-,2Ei2, E + ,Eis, E 23 ), (B-,2Bi2,B + , Bi 3 , B23)). 

Here E± = En ± E 2 2, and B± = Bn ± B22 and it is understood that the 
relations g 1 ^ E^ = and <? y Bij = are used everywhere to replace the fields 
E 33 and B 33 by our unknowns. Written out explicitly, this system takes the 
form (F + A^ 1 ) df, u = b, with 



I" 





" 




A^ 1 " 







A' 1 = 


T A» 



where 



"10000" 







-e^ 3 





e^ 2 


e"i 
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-e"i 


e" 2 
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e" 2 


-c"i 
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e"i 
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-e"i 


-e^ 2 


e"i 









(4.93) 



The reduced system consisting of ( 4.80| ) to (4.92), is symmetric hyperbolic. 
However, beyond that the choice of this particular system was motivated by 
the following specific features. 

(i) The theory of maximally dissipativc initial value problems applies to 
our reduced equations. In equations (4.8C) to (4.91) the derivative d x s, which 
by our gauge conditions occurs only with the directional derivative e 3 , is ap- 



plied to the gauge source functions but not to the unknowns, while (4.93) 
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shows that we have perfect control on the non-trivial part of the normal ma- 



trix arising from (4.92). Our discussion of maximally dissipative initial value 



problems, which led to (4.73), and the form of the matrices (4.93) suggest 
that we can prescribe besides the datum %, which characterizes the boundary, 
precisely two free functions as boundary data on T. This is confirmed by the 
detailed discussion in p2| , though in general a number of technical details 
have to be taken care of . 



(ii) If instead of ( 4.92 ) we had chosen the system = 0, Qui — as 
equations for the electric and magnetic part of the conformal Weyl tensor, 
the theory of maximally dissipative initial value problems would also have 
applied. We would, however, have come to the conclusion that besides the 
mean extrinsic curvature four functions could be prescribed freely on T. This 
apparent contradiction is resolved when one tries to show the preservation 



of the constraints, i.e. that those equations contained in ( 4.74 ) are satisfied 
which are not already solved because of the gauge conditions and the reduced 
equations. In the case of the reduced equations above this can be shown for 
the following reason. The subsidiary system splits in this case into a hierachy 
of symmetric hyperbolic subsystems with the following property. The first 
subsystem has vanishing normal matrix on T. This implies under suitable 
assumptions on the domain of the solution to the reduced equations that all 
unknows in this subsystem must vanish, because the data on S are of course 
arranged such that all constraints are satisfied. Furthermore, it follows for any 
subsystem in the hierarchy that its normal matrix vanishes if the unknowns 
of all previous subsystems in the hierarchy vanish. From this the desired 
conclusion follows in a finite number of steps. If we had considered instead 
the system Pij = 0, Qki = 0, the discussion whether the constraints are 
preserved would have become quite complicated and would have led us in 
the end to the conclusion that only two functions are really free on T while 
the others are subject to restrictions determined by the evolution properties 
of the reduced system. 

We end our discussion of the initial boundary value problem with an 
observation about the characteristics of the reduced system. Equations (4.! 



to (4.91) contribute a factor of the form 

to the characteristic polynomial (using again only the frame components of 
the covector £). The corresponding characteristics are timelike or null with 



respect to g^ v . However, the subsystem (4.92) contributes a factor 

If we denote by or 7 ' the 1-forms dual to the vector fields , so that < c- 7 , efc > 
= 53 k and g^v = ct° er° — a x a x v — a^crl — , the characteristics associated 
with the third factor in the polynomial above can be described as the null 
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hypersurfaces with respect to the metric 

^ = °l ^ 4- ( 4 -94) 

4.4 The Einstein Dirac System 

Apparently, not much has been shown so far about the existence of solutions 
to the Einstein-Dirac system for general data. The initial value problem for 
this system was considered in |)| , but no existence theorem for the evolution 
equations was proved there. J. Isenberg has suggested to us that it should 
be possible to show the well-posedness of the equations by formulating the 
equations as a system of wave equations in a way similar to what was done 
for the Cauchy problem for classical supergravity in ||. However, this idea 
has not been worked out in the literature. 

We shall indicate here how to obtain symmetric hyperbolic evolution equa- 
tions from the Einstein-Dirac system. To avoid lengthy calculations, we shall 
not discuss the complete reduction procedure but only use this system to 
illustrate certain questions arising in the reduction. 

We shall write the equations in terms of the 2-component spin frame for- 
malism, which may be thought of as the spinor version of the frame formalism 
used in the previous sections. The fields and the equations will be expressed 
in terms of a spin frame {t a }a=o,ij which is normalized with respect to the 
antisymmetric bilinear form e in the sense that it satisfies e a b — e(taj L b) with 
eoi = 1. The associated double null frame is given by e aa ' = i a l a < , it satisfies 
e aa ' = e aa > and g{e aa >,ew) = e a bZa'V 

All spinor fields (with the possible exception of the basic spin frame itself 
and the vector fields e aa ' ) will be given with respect to the spin frame above 
and we shall use e a b and e ab , defined by the requirement e a b€ cb = S a c (the 
Kronecker symbol), to move indices according to the rule ui a = e ab u>b, ui a = 

< rb r 

w e bQ . 

We use the covariant derivative operator V acting on spinors, which is 
derived from the Levi-Civita connection of g and satisfies V e a b — 0, to 
define connection coefficients r aa 'bc = ^aa'(6c) by 

^ ' aa' t*b — V 7 e aa r ^b — F aa r 5 L c . 

For any spinor field uj a we have 

(V cc ' V 'dd' — V 'dd' V cc ') LO a — —R a bcc'dd' ^ b — T cc > ee dd' V ee 'W a , 

with vanishing torsion 

= Tbf dd ce> £dd' = e C c> - V CC ' ebb' - [ebb',e C c'} (4.95) 
= -Tfeb' d c e-dc' + rbb> d d e c d' — r cc i d b edb' — r cc > d b< ebd' ~ [ebb',e cc '], 
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and curvature spinor field 

Rabcc'dd' — P-dd'{r c c'ab) ~ e cc ' (/dd'ab) + Rdd'ac Ac' & b (4.96) 

~~\~Red'ab Red d Rcc'ae Rdd 1 b ^ec' ab Rdd' c 
~~ t - de'ab * cc' d' ^ ce' ab 1 dd' c' ^ cc' c£d' ^ ce' ab- 

The latter has the decomposition 

Rabcc'dd' = ~^abcd^c'd' ~ $abc' d' ^cd + A € c > d' {^bd ^ac + € ad ^bc) , (4.97) 

into the conformal Weyl spinor field tfabcd = ^(abcd) as well as the Ricci spinor 
^aba'b' = ^ (ab) (a 1 b' ) = ^aba'b' and the scalar A, which allow us to represent 
the Ricci tensor in the form 

Raa'bb' = 2 <Paba>b> + 6 A tab t-a'b' ■ 

The Bianchi identity reads 

V a '%bcf=V (a f ' QbcW- (4-98) 

The Field Equations The Einstein-Dirac system is specified (cf. @|) by 
a pair of 2-spinor fields <p ai Xa' satisying the Dirac equations 

V a al <f> a = fiXa', V a a 'xa'=M0a, (4.99) 

with a real constant n, and the Einstein equations with energy-momentum 
tensor 

Taa'bb' = "y {<Pa ^bb'4>a> ~ 4>a> ^bb'4>a> + 4>b V a a'<?V ~ 4>b' V ' a a> '4>b (4.100) 
~Xa ^bb'Xa' + Xa' ^bb'Xa ~ Xb ^aa'Xb' + Xb> ^aa'Xb}- 

The Einstein equations then take the form 

A = --^{ct>aX a -k>X a '), (4.101) 

@aba'b> = ^ i>(a ^b)(a'4>b') ~ 4>(a' V^)( a 6 ) (4.102) 
~X(a ^b)(a'Xb') + X(a> ^b')(aXb)}- 

The discussion of this system is complicated by the fact that the Dirac 
equations are of first order while derivatives of the spinor fields also appear 



on the right hand side of (4.102). Consequently, the right hand side of the 
Bianchi identity (4.98) is given by an expression involving the derivatives of 
the spinor fields from zeroth to second order. Therefore we need to derive 
equations for these quantities as well. 
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By taking derivatives of the Dirac equations and commuting derivatives, 
we obtain 

V'a'Vwta = -R h a a a>W0h+^wXa> (4.103) 



V Q u'VccfVwta = -^crjR h a ° a'W <t>h (4.104) 
~R h a a a'bb' ^cc'4>h + bb'V cc'Xa' ~ R fc ° a' cc' ^hb'^a 



-R 



h' 



V a' cc' V bh' 



ba R a a' 



cc' Vbb'Cph 



and similar equations for the derivatives of Xa' ■ It is important here that the 
curvature quantity 

R k a " a'bb' = — $ da'b' —3/1 ta'b'^d 

which occurs in these equations does not contain the conformal Weyl spinor. 
We can use (4.101) and (4.102) to express R h a a a'bv an d its derivative in 
(4.103) and (4.104) in terms of the spinor fields and their derivatives to 
obtain a complete system of equations for cf> a , Vf,f/(/> a , V cc ' Vf,f/</> a and the 
corresponding fields derived from Xa' ■ 

These fields are not quite independent of each other. If we define sym- 
metric fields (f) aC a', 4>abca'V, Xaa'c', Xaba'b'c' by Setting 

a' c' t-7 (a' c') a b' c v7 Co t-j b' c') 

Xa = V a l \ ', Xab =V(a y V b) X , 

we get from the Dirac equations 



Vaa"Pb — Vaba' ~ -^^abXa', 

1 h 2 h 

V cc 'Vfcb/0 a = 4>abcb'c' — ^ £ b'c' $ ' abch <fi + ^ £ c (a ^fc)/ifc'c' <t> 

2 1 2 

+ 2 A (j)( a e b ) c tb'd + g M e a(b Xc)b'c' - 4>a £bc ^b'c' , 

and similar relations for the derivatives of Xa' ■ From these, the equations 
above, and (4.98|), (4.101) and (4.102) we can derive equations of the form 



lvl bca'b'c'i 



a '4>abb' — Mba'b'i V* a '<f>abcb'l 
V a ° Xba'V — Nabb'i V a a Xbca'b'c' = ^abcb' c' 



(4.105) 
(4.106) 



where M^ a>b> , N^ w denote functions of <f) a , Xa', <Paca', Xaa'c', while M^ ca , b , c ,, 
N abcb'c' depend in addition on 4> a bca'b', Xaba'b'c', and ^ a bcd- Note that this 
introduces (or rather makes explicit) further non-linearities. 
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Hyperbolic Equations from the Einstein Dirac System Our field 
equations for the unknowns 



e" fe , 



r, 



aa'bc i 



bed? 



4*abb' j 4 > abcb t c'j \a' j Xba'b' 7 Xbca'b'' 



are now given b y (|4. 95| ) and ( 4.97 ) (wi th the left h and si de understood as 
being give n by (p~96| )) , (pffi , fl4 99|), ( J4.105| ) and Q4.106| ). Here Einstein's 
equations (4.101) and (4.402) are used to express quantities derived from the 
Ricci tensor in terms of the spinor fields and their derivatives. 

When we try to deduce a hyperbolic reduced system from these equations, 
the first two equations, which determine the gauge dependent quantities, will 
require the choice of a gauge, while we expect the remaining equations, which 
are tensorial, to contain subsystems which are hyperbolic irrespective of any 
gauge. This is indeed the case and there are, due to the fact that most of the 
equations are overdetermined, various possibilites to extract such systems. 

In p6[ sytems of spinor equations have been considered which are built 
from systems of the type 

or their complex conjugates, where /? denotes a multi-index of some sort. 
If the components corresponding to different values of the indices b, j3 are 
independent of each other the equations 



1pbl3 = F VI3 

form a symmetric hyperbolic system. Equations ( 4.99| ) are thus symmetric 
hyperbolic as they stand. If symmetries are present which relate the index 
b to a group of unprimed spinor indices comprised by /3, the equations to 
be extracted are slightly different. For instance, we obtain from (4.98) a 
symmetric hyperbolic system (regarding all fields besides the Weyl spinor 
field as given) of form 

V f V&ablf - V f Q'^aWf = ..., 

where the symmetry ^abcd = ^(abcd) is assumed explicitly so that there are 
five complex unknown functions. Equations (4.405) and ( 4. 106| ) can be dealt 
with similarly. 

To compare the characteristics of the system above with previous hyper- 
bolic systems extracted from the Bianchi identity, we set 



1 



eo 



V2 



(eoc + en'), ei 



4 



V2 



(eov + eic), 
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— i 1 
e2 = ~~?^( eo1 ' ~ ei0 ')' e3 = y^ 00 ' ~~ 6l1 '^' 

Then the characteristic polynomial is given up to a positive constant factor 

by 

£„e" fc P X^3 A,T 6^, (4.107) 
which contains the degenerate quadratic form 



The cone {fc pl/ £ p = 0} is the product of a 2-dimensional cone in the plane 
{£m eM 3 — 0} with the real line so that its set of generators is diffeomor- 
phic to 5* 1 x R. The associated characteristics are timelike. The special role 
played here by the vector field e% allows us to adapt the system to situations 
containing a distinguished direction. 

Another method to extract symmetric hyperbolic equations from spinor 
equations has been discussed in |37f| . It is based on the space-spinor formalism 
in which an arbitrary normalized timelike vector field is used to express all 
spinor fields and spinor equations in terms of fields and equations containing 
only unprimed indices. If the fields and equations are then decomposed into 
their irreducible parts (a direct, though somewhat lengthy algebraic proce- 
dure), the equations almost automatically decompose into symmetric hyper- 
bolic propagation equations and constraints. 

For simplicity we choose the timelike vector field to be 

y/2 e = r aa ' e aa > with w = e a e < + t\ ° e v a ' . 

Since T aa / T ba = e a etc., maps generalizing the map uv — ► r a ° uj a i to 
spinors of arbitrary valence are bijective and allow us to obtain faithful rep- 
resentations of all spinor relations in terms of unprimed spinors. Writing 
V a (, = Tb a V aa / = \e a bP + T^ab, we obtain a representation of the co- 
variant derivative operator in terms of the directional derivative operators 
P = r aa V aa / , V a b = 77f, a V w acting in the direction of eo and in di- 



rections orthogonal to eo respectively. In particular, (4.98) splits under the 
operations indicated above into 'constraints' 

vf^ abfg = ..., 

and 'evolution equations' 

P^ ab cd-^{J^bcd)f = ■■■ 

If the latter are multiplied by the binomial coefficients ( a+b ^_ c+d ) , they are 
seen to be symmetric hyperbolic. The characteristic poynomial of this system 
is again of the form 
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however, since there is no privileged spacelike direction singled out here, we 
have a non-degenerate quadratic form 



k^ = (l + c) e%e% 



(4.108) 



with some constant c > 0, as in the case of the system used in the case of the 
Einstein-Euler equations. 

By the method outlined above symmetric hyperbolic systems are also 
obtained for equations (4.105) and ( 4.1061) . For instance, from the first of 
equations ( 4.105 ) we obtain for (j> abc = t c b 4> a bb' — <P(ab)c an equation of the 
form V a d4>abc = ■ ■ where we only indicate the principal part. Using the 
decomposition 



'Jabc 



u abc 



e c(a <Pb) 



<j)t with 4>* abc = 0( o6c ), Cj>* b = cj>fb f 



and the decomposition of V a f,, we get a system of the form 



■ V" 



+ 2V 



be ,/* 
a be 



2V 



be 



%ac 



a + b + c 



Pi 



abc 



a + b + c 



2V d 



(a 4>bc)d ~ n V (ab 0*) 



(o<P|/|6c) 



which is symmetric hyperbolic. 

It is well known that equations for spinor fields of spin y , m > 2 give rise 
to consistency conditions (cf. J72|l ). For instance, the equation 

V a a> <Pabcb>c> = H bca , b , d , (4.109) 

where we consider the right hand side as given, implies the relation 

* abc d + 4 Kbccv(bi $ aM ' c>) = V aa ' H aa 



Jabcb' c' 



ida'b' d ) 



which reduces e.g. in the case of vanishing right hand side to a particular 
relation between the background curvature and the unknown spinor field. 
Depending on the type of equation and the background space-time, such 
consistency conditions may forbid the existence of any solution at all. Nev- 
ertheless, equation (4.109) implies a symmetric hyperbolic system for which 
the existence of solutions is no problem. Difficulties will arise if one wants to 
show that the constraints implied by (4.109) are preserved. 

This example emphasizes the need to show the preservation of the con- 
straints. Because in our case the right hand sides of the equations are given by 
very specific functions of the unknowns themselves, we can expect to obtain 
useful subsidiary equations. 
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There are again various methods to obtain hyperbolic equations for the 
gauge-dependent frame and connection coefficients, which depend in partic- 
ular on the choice of gauge conditions. In JItJ the coordinates and the 
frame field have been subject to wave equations (nonlinear in the case of 
the frame field). Here we shall indicate a gauge considered in |38| (which 
can, of course, also be implemented in the frame formalism considered in the 
previous sections). 

We shall denote by T a 'time flow vector field' and by x^ coordinates 
on some neighbourhood of an initial hypersurface S. We assume T to be 
transverse to S, the 'time coordinate' t = x° to vanish on S, and the relation 
< dx M ,T >= <5 M o to hold on the neighbourhood such that we can write 
T = d t . 

The frame e aa ' is chosen such that the timelike vector field T aa e aa > is 
orthogonal to S. Using the expansion 



1 



= - T aa < T 



2 



r b a ' e ab with e ab = r (b a e a)a /, 



we can write 



with 



T = a t cc e cc , + /3 CC ' e cc , = a r cc e cc , + (3 ab e ab 



r cc /3 cc , =0, P ab = T (a a p b)a ,. (4.110) 

Thus the evolution of the coordinates off S is determined by the fields a^O 
and f3 aa and we can write 

t cc ' e cc , = -(d t -(3 ab e ab ). (4.111) 
a 

Since we have Vt l c = r b c i b , the evolution of the frame is determined 
by the functions 

ri rpaa y 1 

be ^ aa'bc- 



and we can write 



r aa ' r aa , bc = - (r bc - p ae r aebc ), (4.112) 

a 



with r ae bc — T( e a r a ) a , bc . 

We now consider the fields a = a(x tJl ) > 0, (3 ah — I3 ah {x^) (together four 
real functions) as 'coordinate gauge source functions' and the field r bc = 
r bc (x^) (six real functions) as 'frame gauge source functions'. This is feasible, 
because given these functions, we can find smooth coordinates and a frame 
fields close to an initial hypersurface such that the given functions assume 
the meaning given to them above. 

The gauge conditions are then expressed by the requirement that the 



right hand sides of (4.111) and (4.112) are given in terms of the gauge source 
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functions and e a b and r ae bc- Thus it remains to obtain evolution equations 
for e M ab and r aebc . 

Reading the quantity e M aa > for fixed index fi as the expression of the 



differential of x M in the frame e aa ', we can write (4.95) in the form 

Vaa'S^w - Vw^ aa> = 0. 

Contracting this equation with T aa and r c b and symmetrizing, we obtain 
the equation 

= V T cb - ef b , {b V T t d) b ' - e» aa ,V bc T aa \ 
which can be rewritten in the form 

d t e' 1 ab = ■ ■ ■ , 

where the right hand side can be expressed in terms of the gauge source 



functions and their derivatives and the unknowns. By using (4.97) with (4.96) 



on the left hand side, we can derive in a similar way an equation 

dt r ae bc = ■ ■ ■ , 

with the right hand side again being given in terms of the gauge source 
functions and their derivatives and the unknowns. 

Thus we obtain symmetric hyperbolic reduced equations for all unknowns 



except those given by the left hand sides of the gauge conditions (4.111) and 



(4.112). Our procedure applies of course to various other sytems. Our choice 
of gauge is of interest because of the direct relation between the gauge source 
functions and the evolution of the gauge. The causal nature of the evolution 
can be controlled explicitly because the formalism allows us to calculate the 
value of the norm g{T 7 T) = 2a 2 + (3 a bf3 ab - This may prove useful if it is 
desired to control the effect of the choice of gauge source functions on the 
long time evolution of the gauge in numerical calculations of space-times. 



4.5 Remarks on the Structure of the Characteristic Set 

We have seen that for certain reduced systems there occur besides the 'physi- 
cal' characteristics, given by null hypersurfaces, also characteristics which are 
timelike or spacelike with respect to the metric g^. Timelike characteristics, 
which usually occur if a system of first order is deduced from a system of sec- 
ond order, are usually harmless and of no physical significance. The spacelike 
characteristics, which are partly due to the choice of gauge condition and 
partly due to the use made of the constraints, have no physical significance 
either. Though they are associated with non-causal propagation, there is a 
priori nothing bad about them and it rather depends on the applications one 
wants to make whether they are harmful or not. 
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In the characteristic polynomial (4.71) of the reduced equations for the 



Einstein-Euler system there appears a factor £g + \ h ab £ a £b which corre- 
sponds to timelike characteristics. In vacuum the corresponding cone has of 
course no physical meaning, since in general there is no preferred timelikc 
vector field available. In the perfect fluid case there is a distinguished time- 
like vector field present. This, and perhaps the symmetry of the inner cone 
with respect to the fluid vector, has led some people to speculate on the phys- 
ical significance of that cone |n]. However, in some of the later examples of 
hyperbolic equations deduced from the Bianchi identity, which could also be 
used in the fluid case, the structure of the char acter istics is drastically differ- 
ent from the on e obse rved in the fluid case (cf. ( 4.94 ) and also the degenerate 
cone arising in ( 4.107 )). In particular, the factor above does not occur in their 
characteristic polynomials. The large arbitrariness in extracting hyperbolic 
equations, which arises from different use made of the constraints implied by 
the Bianchi identity, suggests that in the case of the Einstein-Euler system 
the only 'physical characteristics' are those associated with the fluid vector, 
the null cone of g^, and the sound cone. 

The null cone of the metric ( 4.94 ) touches the null cone of the metric g 
in the directions of ± but it is spacelike in all other directions. Thus all 
null hypersurfaces of it are spacelike or null for g^ v . Such a cone has the 
effect that the 'domain of uniqueness' defined by the techniques discussed in 
Sect. may decrease. However, as we have seen, is does not prevent us from 
proving useful results. 

There is also no reason to assume that the additional characteristics nec- 
essarily create problems in numerical calculations. In situations where the 
maximal slicing condition can be used, the occurrence of spacelike character- 
istics which are related, as in our examples, in a rigid way with the metric 
should be innocuous. Also, numerical calculations based on equations with 
inner characteristic cones as observed above have been performed without 
difficulties (11, |§). 



5 Local Evolution 



5.1 Local Existence Theorems for the Einstein Equations 

The purpose of this section is to present a local existence theorem for the 
Einstein vacuum equations. By (abstract) vacuum initial data we mean a 
three-dimensional manifold S together with a Riemannian metric h a b and a 
symmetric tensor \ab on S which satisfy the vacuum constraints (see Sect. ^|). 
A corresponding solution of the vacuum Einstein equations is a Lorentzian 
metric g a p on a four-dimensional manifold M and an embedding <f) of S into 
M such that h a b and Xab coincide with the pull-backs via (/) of the induced 
metric on (f>(S) and the second fundamental form of that manifold respectively 
and the Einstein tensor of g a p vanishes. If 4>(S) is a Cauchy surface for the 
space-time (M, g a p) then this space-time is said to be a Cauchy development 
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of the data (S,h a b,Xab)- The basic local existence theorem for the vacuum 
Einstein equations says that every vacuum initial data set has at least one 
Cauchy development. In fact to make this precise it is necessary to fix the 
differentiability properties which are assumed for the data and demanded of 
the solution. For instance, the result holds if the differentiability class for 
both data and solutions is taken to be C°°. Note that there is no need to 
require any further conditions on the spatial dependence of the data. 

The proof of local existence will now be outlined. We follow essentially the 
original method of |l(| except for the fact that we reduce second order equa- 
tions to first order symmetric hyperbolic systems and that we use harmonic 
mappings rather than harmonic coordinates. The use of harmonic mappings, 
as discussed in Sect. |], allows us to work globally in space even if the man- 
ifold S cannot be covered by a single chart. Using harmonic coordinates it 
would be necessary to construct solutions local in space and time and then 
piece them together. Choose a fixed Lorentz metric on S x R, for instance 
the metric product of the metric h a b with —dt 2 . This comparison metric will 
be denoted by g a p. The idea is to look for a solution g a p on an open sub- 
set U of R x S such that the identity is a harmonic map from (U, g a /3) to 
{U,g a p). This is a condition which is defined in a global invariant way. Its 
expression in local coordinates is g^irS 1 — F g L) = where Fg and Fg 
are the Christoffcl symbols of g a/ 3 and g a p respectively. In the terminology 
of Sect. [O] this means that we choose g^F^ as a gauge source function. 

Next consider the question of reduction of nonlinear wave equations to 
symmetric hyperbolic form. This is done as follows. Let g a ^(t,x,u) be func- 
tions of (t, x, u) which for each fixed value of (t, x, u) make up a symmetric 
matrix of Lorentz signature and consider an equation of the form: 

g af} d a d p u + F{t, x, u, Du) = (5.1) 

This has been formulated in a local way but a corresponding class of equations 
can be defined in the case that the unknown u is a section of a fibre bundle. 
As in the above treatment of symmetric hyperbolic equations on a manifold, 
consideration will be restricted to the case of sections of a vector bundle V . 
Choose a fixed connection on V. Then the class of equations to be considered 
is obtained by replacing the partial derivatives in the above equation by 
covariant derivatives defined by the given connection. Let u a — V Q u. Then 



the equation (5.1) can be written as: 



-g 00 V u - 2. 9 0a V a ii = g ab V a u b + F(t, x, u , u„) 
g ab V Q u a = g ab W a u + K b 
V U = u 

Here K b is a term involving the curvature of the connection which is of order 
zero in the unknowns of the system. This is a symmetric hyperbolic system 
for the unknowns u and u a . Since u is allowed to be a section of a vector 



The Cauchy Problem 77 



bundle we are dealing with a system of equations. However the functions g a ° 
must be scalars. The appropriate initial data for the second order equation 
consists of the values of u and <9 t u on the initial hypersurface. From these 
the values of the functions u a on the initial hypersurface may be determined. 
Thus an initial data set for the symmetric hyperbolic system is obtained. It 
satisfies the additional constraint equation V Q u = u a . Applying the existence 
theory for symmetric hyperbolic systems gives a solution (u,u a ). To show 
that the function u obtained in this way is a solution of the original second 
order equation it is necessary to show that the constraint equation is satisfied 
everwhere. That Wqu = uq follows directly from the first order system. It also 
follows from the first order system that Vo(u a — V a it) = 0. Since this is a first 
order homogeneous ODE for u a — V a u, the vanishing of the latter quantity 
for t — implies its vanishing everywhere. 

In the case of the vacuum Einstein equations the bundle V can be taken 
to be the bundle of symmetric covariant second rank tensors. The connection 
can be chosen to be the Levi-Civita connection defined by g a p. This is only 
one possible choice but note that it is important that this connection does 
not depend on the unknown in the equations, in this case the metric g a p- 

Now a proof of local in time existence for the vacuum Einstein equations 
will be presented. Let h a b and Xab be the initial data. In Sect. it was shown 
that the vacuum Einstein equations reduce to a system of nonlinear wave 
equations when harmonic coordinates, or the generalization involving gauge 
source functions, are used. As was already indicated in that section, there is 
no loss of generality in imposing this condition locally in time. If there exists 
a development of particular initial data then there exists a diffeomorphism <fi 
of a neighbourhood of the initial hypersurface such that the pull-back of the 
metric with the given diffeomorphism satisfies the harmonic condition with 
respect to g a p- In fact <f> can be chosen to satisfy some additional conditions. 
The harmonic condition is equivalent to a nonlinear wave equation for <p. 
Solving the local in time Cauchy problem for this wave equation provides 
the desired diffeomorphism. The existence theory for this Cauchy problem 
follows from that for symmetric hyperbolic systems by the reduction to first 
order already presented. The initial data for <j) will be specified as follows. It is 
the identity on the initial hypersurface and the contraction of the derivative 
of 4> with the normal vector with respect to g a p should agree on S with the 
corresponding quantity constructed from the identity on M . Since the vector 
d/dt is the unit normal vector to S with respect to g a p it will also have this 
property with respect to g a p. 

A local solution of the Einstein equations corresponding to prescribed ini- 
tial data can be obtained as follows. Let h a b and \ab denote the components 
of the tensors making up the initial data in a local chart as above. A set of 
initial data for the harmonically reduced vacuum Einstein equations consists 
of values for the whole metric g a p and its time derivative on the initial hy- 
persurface. A data set of this kind can be constructed from h a b and Xab as 
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follows. (The following equations are expressed in local coordinates, but their 
invariant meaning should be clear.) 

9ab — h a b, gOa = 0, #00 = —1 

dtg ab = 2 Xab , d t goa = h bc (h ab . c - (1/2) Va) - h ab r^ d h cd , 

d t gov = -2h ab Xab 

These data are chosen in such a way that the harmonic condition is satisfied 
on the initial hypcrsurface. 

Corresponding to the initial data for the reduced equations there is a 
unique local solution of these equations. It remains to show that it is actually a 
solution of the Einstein equations provided the initial data h ab and \ab satisfy 
the Einstein constraint equations. In order to do this it suffices to show that 
the harmonic conditions are satisfied everywhere since under those conditions 
the reduced equations are equivalent to the Einstein equations. Let A a = 
gP-f^pa^ _ pp ). That the harmonic conditions arc satisfied can be verified 
using the fact that the quantities A a satisfy a linear homogeneous system of 
wave equations. By uniqueness for this system the A a vanish provided the 
initial data A a and d t A a vanish on the initial hypersurface. The first of these 
was built into the construction of the data for the reduced equations. The 
second is a consequence of the combination of the reduced equations with the 
Einstein constraints. 

5.2 Uniqueness 

The argument of the last section gives an existence proof for solutions of the 
vacuum Einstein equations, local in time. It does not immediately say any- 
thing about uniqueness of the space-time constructed. The solution of the 
reduced equations is unique, as a consequence of the uniqueness theorem for 
solutions of symmetric hyperbolic systems. However the freedom to do diffco- 
morphisms has not yet been explored. In fact it is straightforward to obtain a 
statement of uniqueness of the solution corresponding to given abstract initial 
data, up to diffcomorphism. Suppose two solutions g\ and g 2 with the same 
initial data are given. Choose a reference metric g as before and determine 
diffeomorphisms <j>i and <p2 such that the identity is a harmonic map with 
respect to the pairs (g,gi) and (3,32) respectively. The transformed metrics 
satisfy the same system of reduced equations with the same initial data and 
thus must coincide on a neighbourhood of the initial hypersurface. Thus g\ 
and (02 (t>\ 1 )*g2 coincide on a neighbourhood of the initial hypersurface. 
This statement is often referred to as 'geometric uniqueness'. 

For hyperbolic equations it is in general hard to prove theorems about 
global existence of solutions due to the possibility of the formation of singu- 
larities. On the other hand, it is possible to prove global uniqueness theorems. 
Proving global uniqueness for the Einstein equations is more difficult due to 
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difficulties with controlling the freedom to do diffeomorphisms. In this con- 
text we use the word 'global' to mean not just applying to a subset of a 
given solution, but to whole solutions having suitable intrinsic properties. 
The strategy which has just been used to prove local existence using har- 
monic maps cannot be applied directly. For doing so would require a global 
existence theorem for harmonic maps and that we do not have in general. 

A global uniqeness theorem was proved by abstract means by Choquet- 
Bruhat and Geroch who introduced the notion of the maximal Cauchy 
development. If initial data for the Einstein equations coupled to some matter 
fields on a manifold S are prescribed, a development of the data is a solution 
of the Einstein-matter system on a manifold M together with an embedding 
<j> of S into a M which induces the correct initial data and for which the 
image of S is a Cauchy surface. Another development with a solution of the 
Einstein-matter equations on a manifold M' and an embedding <f>' is called 
an extension of the first if there is a diffcomorphism tp from M to an open 
subset U of M' which maps the given metric on M onto the restriction of 
the metric on M' to U and also maps the matter fields on M to those on U 
obtained by restriction from M'. If we are dealing with the vacuum Einstein 
equations then the requirement on the matter fields is absent. In |fl~sf] the 
following theorem was proved for the vacuum case: 

Theorem Let S be an initial data set. Then there exists a development M of 
S which is an extension of every other development of S. This development 
is unique up to isometry. 

The development whose existence and uniqueness is asserted by this theorem 
is called the maximal Cauchy development of the initial data set. Uniqueness 
up to isometry means the following. If we have two developments of the same 
data given by embeddings <f> and <f>' of S into manifolds M and M' respectively 
then there exists a diffeomorphism t/j : M — > M' which is an isometry and 
satisfies <f>' = <j> o ip. The proof of this theorem does not depend strongly on 
the vacuum assumption. One potential problem in extending it to certain 
matter fields is gauge freedom in those fields. This requires the concept of 
extension to be defined in a slightly different way. For instance in the case 
of gauge fields it is necessary to consider not only diffeomorphisms of the 
base manifold, but also automorphisms of the principal bundle entering into 
the definition of the theory. We do not expect that this leads to any essential 
difficulty, but in any concrete example one should pay attention to this point . 

The proof of this theorem applies directly to Zorn's lemma and it is an 
open question whether it is possible to remove the use of the axiom of choice 
from the argument. The maximal Cauchy development is often very useful in 
formulating certain arguments. However it remains very abstract and gives 
the subjective impression of being difficult to pin down. 

It should be emphasized that, despite its global aspects, it would be mis- 
leading to consider the above theorem as a global existence theorem for the 
Einstein equations in any sense. A comparison with ordinary differential equa- 
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tions may help to make this clear. If an ordinary differential equation for 
a function u(t) is given (with smooth coefficients) then the standard local 
existence theorem for ordinary differential equations says that given an ini- 
tial value uq there exists a T > and a unique solution u(t) on the inter- 
val (— T, T) with u(0) = Uo- Now one can ask for the longest time interval 
(— 7\, T 2 ) on which a solution of this kind exists for a fixed initial value uq- 
This is called the maximal interval of existence and has a similar status to 
that of the maximal Cauchy development. The fact that the maximal inter- 
val of existence is well-defined says nothing about the question whether the 
solution exists globally or not, which is the question whether T\ and T 2 are 
infinite or not. The existence of the maximal Cauchy development says noth- 
ing about the global properties of the solution obtained in this way, whereas 
global existence of the solution of an ordinary differential equation does mean 
that the solution has a certain global property, namely that it exists for an 
infinitely long time. 

5.3 Cauchy Stability 

Cauchy stability of the initial value problem for the Einstein equations is the 
statement that, in an appropriate sense, the solution of the Einstein equations 
depends continuously on the initial data. This continuity statement has two 
parts, which can be stated intuitively in the following way. Firstly, if a solution 
corresponding to one initial data set is defined on a suitable closed region, 
then the solution corresponding to any initial data set close enough to the 
original one will be defined on the same region. Closeness is defined in terms of 
Sobolev norms. Care is needed with the interpretation of the phrase 'the same 
region' due to the diffcomorphism invariance. To make it precise, something 
has to be said about how regions of the different spacetimes involved are to 
be compared with each other. Secondly, the solution defined on this common 
region depends continuously on the initial data, where continuity is again 
defined in terms of Sobolev norms. In non-compact situations it is appropriate 
to use local Sobolev norms for this, i.e. the Sobolev norms of restrictions of 
a function to compact sets. 

Rather than make this precise in general we will restrict to one case where 
the formulation of the statement is relatively simple, but which is still general 
enough to give a good idea of the basic concepts. Consider initial data sets 
for the Einstein equations on a compact manifold S. For definiteness let us 
restrict to the vacuum case. As has been discussed above, solutions can be 
constructed by using the harmonically reduced equations. One step in this 
process is to associate to geometric data (h a b, Xab) full data (g a p, d t g a p)- For 
each of these pairs let us choose the topology of the Sobolev space H S (S) 
for the first member and that of H S ^ 1 (S) for the second. Then standard 
properties of Sobolev spaces show that if s is sufficiently large the mapping 
from geometric data to full data is continuous. Suppose now that we have 
one particular solution of the Einstein vacuum equations with data on S. 
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The corresponding solution of the harmonically reduced equations exists on 
some region of the form S x [— T, T]. If s is sufficiently large then there exists 
an open neighbourhood of the given data in H S (S) x H S ~ 1 (S) such that 
for any data in this neighbourhood there exists a corresponding solution in 
H s (Sx [— T, T]). Moreover the mapping from data to solutions defined on this 
neighbourhood is continuous (in fact differentiable). This has been proved by 
Choquet-Bruhat fl§]. 

The theorem concerning a compact initial hypersurface can also be mod- 
ified to give a local statement of the following type. Let initial data for the 
Einstein equations be given on some manifold S and suppose that a corre- 
sponding solution is given on a manifold M. There is a neighbourhood U of 
S where harmonic reduction is possible globally. This identifies U with an 
open subset of S x R. This contains a set of the form V x [— T, T] (for some 
open subset V of S and some T > 0) which contains any given point of the 
initial hypersurface. If we cut off the the initial data for the harmonically 
reduced equations and use the domain of dependence, we can use the above 
statement for a compact initial hypersurface to get continuous dependence on 
initial data for a possibly smaller set V x [— T', T']. Summing up, each point 
sufficiently close to the initial hypersurface has a neighbourhood W\ with 
compact closure such that there is an open subset Wi of the initial hyper- 
surface with compact closure such that the following properties hold. If the 
restriction of an initial data set for the Einstein equations is sufficiently close 
to that of the original data set in B. s {W-i) then there exists a correspond- 
ing solution of class H s on a neighbourhood of W\. Moreover, the resulting 
mapping from H s (W2) to H S (W\) is continuous. 

5.4 Matter Models 

To specify a matter model in general relativity three elements are required. 
The first is a set of tensors (or perhaps other geometrical objects) on space- 
time which describe the matter fields. The second is the equations of motion 
which are to be satisifed by these fields. The third is the expression for the 
energy-momentum tensor in terms of the matter fields which is to be used to 
couple the matter to the Einstein equations. Note that in general both the 
matter field equations and the expression for the energy-momentum tensor 
involve the space-time metric. Thus it is impossible to consider matter in 
isolation from the space-time metric. In solving the Cauchy problem it is 
necessary to deal with the coupled system consisting of the Einstein equations 
and the equations of motion for the matter fields. 

There are two broad classes of matter models which are considered in gen- 
eral relativity, the field theoretical and phenomenological matter models. The 
distinction between these is not sharply defined but is useful in order to struc- 
ture the different models. The intuitive idea is that the field theoretic matter 
models correspond to a fundamental description while the phenomenological 
models represent an effective description of matter which may be useful in 
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certain situations. Within the context of classical general relativity, which is 
the context of this article, the pretension of the field theoretic matter models 
to be more fundamental is not well founded since on a fundamental level the 
quantum mechanical nature of matter should be taken into account. 

Before going further, a general remark on the Einstein-matter equations 
is in order. Suppose that in any given coordinate system the matter equa- 
tions can be written in symmetric hyperbolic form in terms of a variable u. 
Consider the system obtained by coupling the harmonically reduced Einstein 
equations, written in first order symmetric hyperbolic form, to the sytem 
for u. If the coupling is only by terms of order zero then the combined sys- 
tem is symmetric hyperbolic and a local existence theorem for the reduced 
Einstein-matter system is obtained. The condition for this to happen is that 
the equations for the matter fields contain at most first derivatives of the 
metric (in practice the Christoffel symbols) and that the energy-momentum 
tensor contains no derivatives of u. When these conditions are satisfied, local 
existence for the Einstein-matter equations (not just the reduced equations) 
can be proved using the same strategy as we presented in the vacuum case. 
The fact, which should hold for any physically reasonable matter model, that 
the energy-momentum tensor is divergence free as a consequence of the mat- 
ter field equations, can be is used derive the equation which allows it to be 
proved that the harmonic condition propagates. 

It would be unreasonable to try and describe here all the matter models 
which have ever been used in general relativity. We will, however, attempt to 
give a sufficiently wide variety of examples to illustrate most of the important 
features to be expected in general. We start with the field theoretic models. 

The simplest case is where the matter field is a single real- valued function 
<fi. The equations of motion are: 

V a V a <f> = m 2 4> + V'((l)) 

Here to is a constant and V is a smooth function which is 0(<f) 3 ) for <j) close to 
zero. A typical example would be V(cf)) — (j) 4 . The energy-momentum tensor 
is: 

T af} = VaWrf - [(1/2)(VW 7 $ + m 2 <j> 2 + 2V(4)]g a p 

The equation for <f> is a nonlinear wave equation and so may be reduced 
to a symmetric hyperbolic system. When it is coupled to the harmonically 
reduced Einstein equations via the energy-momentum tensor above and the 
whole system reduced to first order there is no coupling in the principal part. 
As mentioned above this is enough to allow a local existence theorem to 
be proved. Note that the splitting off of the nonlinear term V and the sign 
condition following from the form to 2 of the coefficient in the linear term 
are irrelevant for the local well-posedncss of the equations. On the one hand 
they are motivated by considerations of the physical interpretation of the 
equations. On the other hand they have an important influence on the global 
behaviour of solutions. This matter model is often referred to as 'the scalar 
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field' , although when used without qualification this often means the special 
case m = 0, V = 0. 

The scalar field can be thought of as a mapping into the real line. It 
can be generalized by considering mappings into a Riemannian manifold N. 
An equivalent of the massless scalar field with vanishing potential is the 
nonlinear er-model or wave map as it is often known to physicists and math- 
ematicians respectively. It is a mapping from space-time into a manifold N 
with Riemannian metric h called the target manifold. The field equations 
and energy-momentum tensor have a coordinate-invariant meaning but we 
will content ourselves with giving the expressions in coodinate systems on M 
and N. The field equations are 

v Q v Q ^ A + r£ c (^)v a< £ B w c = o 

where r^ c are the Christoffel symbols of h in some coordinate system. The 
energy-momentum tensor is: 

T a p = h AB [V a ^ A V P 4> B - (1/2)(V 7 ^V^ b ) 3q/3 ] 

The special case where N is the complex plane with the flat Euclidean metric 
corresponds to the complex scalar field. In contrast to the case of the scalar 
field, the wave map does not allow the addition of a mass term or a potential 
term in any obvious way. If N has the structure of a vector space there is an 
obvious way of defining a mass term and further structure on N may lead 
to natural ways of defining a potential. These features occur in the case of 
Higgs fields. It may be noted that the wave maps considered here, which are 
sometimes also called hyperbolic harmonic maps, are related mathematically 
to the harmonic gauge discussed in Sect. £| The role of the connection P in 
Sect. ^ is played here by the Levi-Civita connection of the target manifold. 

One of the most familiar matter models in general relativity is the Maxwell 
field. This is described by an antisymmetric tensor F a p. The equations of mo- 
tion for the source- free Maxwell field are V a F af3 = and V Q F ( g 7 + V ' pF ia + 
V~/F a p = and the energy-momentum tensor is 

T a p = F a lF M - {l/A)F^F lS g aP 

The second set of Maxwell equations can be solved locally by writing F a p — 
V ' a Ap — V pA a for a potential A a . Then the other equations can be regarded 
as second order equations for A a . Note, however, that if space-time has a 
non-trivial topology then it may be impossible to find a global potential 
which reproduces a given field F a p. In the same way that the scalar field 
can be generalized to get wave maps, the Maxwell field can be generalized to 
get Yang-Mills fields. We will not give the global description of these fields 
involving principal fibre bundles but only give expressions in local coordinates 
and a local gauge. The model is defined by the choice of a Lie algebra (which 
we describe via a basis) and a positive definite quadratic form on the Lie 
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algebra with components h IJ in this basis. Let Cj K be the structure constants 
in this basis. The basic matter field is a one-form with values in the Lie 
algebra and the field strength is defined by 

Fip = V«4 - V p Ai + C^AiAf 

The field equations are 

V a F Ia ? + Cj K A^F Ka ^ = 

In the special case where the Lie algebra is one-dimensional (and hence 
Abelian) the Yang-Mills equations reduce to the Maxwell equations. Note 
however that the Yang-Mills field cannot be described by the field strength 
alone. The description in terms of a potential is indispensible. 

A complication which arises when studying the initial value problem for 
the Yang-Mills or Einstein- Yang-Mills systems is that of gauge invariance. 
Although the potential is required it is not uniquely determined. Gauge trans- 
formations of the form: 

K ~ Ai + {g- l V a gy 

leave F^g invariant. Here g is a function taking values in a Lie group with 
the given Lie algebra and the expression <? Vg can naturally be identified 
with a Lie-algebra-valued one-form. Fields related by a gauge transformation 
describe the same physical system. The ambiguity here is similar to that of 
the ambiguity of different coordinate systems in the case of the Einstein equa- 
tions. It can be solved in an analogous way by the use of the Lorentz gauge. 
This is similar to the harmonic coordinate condition and reduces the Yang- 
Mills equations on any background to a system of nonlinear wave equations 
which can, if desired, be reduced to a symmetric hyperbolic system. Combin- 
ing harmonic coordinates and Lorentz gauge produces a reduced Einstein- 
Yang-Mills system which can be handled by the same sort of techniques as 
the reduced vacuum Einstein equations. Of course there are a number of 
steps which have to be checked, such as the propagation of Lorentz gauge. 
All these comments apply equally well to the Einstein- Yang-Mills-Higgs sys- 
tem obtained by coupling the Yang-Mills field to a Higgs field. (Now gauge 
transformations for the Higgs field must also be specified.) 

A field theoretic matter model whose Cauchy problem does not fit easily 
into the above framework is that given by the Dirac equation. It has been 
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Probably the best known phenomenological matter model in general rel- 
ativity is the perfect fluid. This has already been discussed at some length in 



Sect. 4.2. Here we mention some complementary aspects. Recall that the ba- 
sic matter fields are the energy density p, a non-negative real- valued function, 
and the four- velocity U a , a unit timelike vector field. The energy-momentum 
tensor is given by 

T al3 = {p + p)U a U fi +pg a0 
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in the signature used here, namely (— , +, +, +). The pressure p is given in the 
isentropic case in terms of an equation of state p = h(p). The usual assump- 
tions on this equation of state is that it is a non-negative continuous function 
defined on an interval [po,oo) which is positive for p > 0. Also it should be 
C 1 for p > with positive derivative h'(p) = dp/dp there. The equations of 
motion of the fluid, the Euler equations, are given by the condition that the 
energy-momentum tensor should be divergence-free. The Euler equations can 
be written as a well-posed symmetric hyperbolic system in terms of the basic 
variables provided we restrict to cases where p > C > for some C > 0. 
Here the condition h' > is crucial. In the study of spatially homogeneous 
cosmological models the equation of state p = kp with k < is sometimes 
considered. We emphasize that a fluid with an equation of state of this kind 
cannot be expected to have a well-posed initial value problem. This has to do 
with the fact that the speed of sound, which is the square root of h', is imag- 
inary in that case. In Sect. 5.5 we prove a related but simpler result, namely 
that in special relativity the Euler equations with this equation of state, lin- 
earized about a constant state, have an ill-posed initial value problem. The 
system obtained by coupling the Euler equations to the harmonically reduced 
Einstein equations can be written as a symmetric hyperbolic system in the 
case that the Euler equations can be written symmetric hyperbolic in terms 
of the basic variables, as stated above. 



In the case where the density is everywhere positive, writing the Euler 
equations in symmetric hyperbolic form is not trivial. One approach is to 
take p and the spatial components U l of the velocity as variables and to 
express U in terms of the U l and the metric via the normalization condition 
U a U a = —1 (see 0). Another possibility (see jlH) is to consider the Euler 
equations as evolution equations for p and U a and treat the normalization 
condition as a constraint, whose propagation must be demonstrated. The 
treatments above are limited to the isentropic Euler equations. For a fluid 
which is not isentropic the conservation equation for the energy-momentum 
tensor must be supplemented by the equation of conservation of entropy 
ii Q V Q s = 0. The equation of state can then be written in the form p — h(p, s) 
or, equivalently, in the form p = f(n,s), where n is the number density of 
particles. We are not aware that a local existence theorem non-isentropic 
Euler equations has been proved by a generalization of the method for the 
isentropic case just outlined, although there is no reason to suppose that 
it cannot be done, provided the condition dh/dp > is satisfied. As we 
saw in Sect. 4.2, the Einstein equations coupled to the non-isentropic Euler 
equations can be brought into symmetric hyperbolic form by introducing 
additional variables in a suitable way, following |3^]. The Cauchy problem 
for the general (i.e. not necessarily isentropic) Euler equations coupled to the 
Einstein equations had much earlier been solved by other means by Choquet- 
Bruhat using the theory of Leray hyperbolic systems ]l7| . 
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If it is desired to show the existence of solutions of the Einstein-Euler 
system representing dynamical fluid bodies (such as oscillating stars) then 
problems arise. Either the density must become zero somewhere, in which 
case the Euler equations as written in the usual variables fail to be symmetric 
hyperbolic there, or at least the density must come arbitrarily close to zero 
at infinity, which is almost as bad. Treating this situation as a pure initial 
value problem is only possible under restrictive circumstances JT6| and the 
general problem is still open. It would seem more promising to try to use the 
theory of initial boundary value problems, explicitly taking account of the 
boundary of the fluid. So far this has only been achieved in the spherically 
symmetric case with po > f6lj . 

One case of a fluid which is frequently considered in general relativity is 
dust. This is defined by the condition that the presure should be identically 
zero. Since in that case hi = the straightforward method of writing the fluid 
equations as a symmetric hyperbolic system does not work. The symmetric 
hyperbolic system of [^| discussed in Sect. [4.2| al so covers the dust case as 
does the existence theorem of Choquct-Bruhat |l7f using a Leray hyperbolic 
system. 

The next phenomenological matter model we will consider comes from 
kinetic theory. It does not quite fit into the framework we have used so far to 
describe matter models since the fundamental matter field is a non-negative 
function / on the cotangent bundle of space-time. (Often the case of particles 
with a fixed mass is considered in which case it is defined on the subset of 
the cotangent bundle defined by the condition g a ^p a P0 = —1, known as 
the mass shell.) The idea is that the matter consists of particles which are 
described statistically with respect to their position and momentum. The 
function / represents the density of particles. The geodesic flow of the space- 
time metric defines a vector field (Liouville vector field) on the cotangent 
bundle. Call it L. This flow describes the evolution of individual test particles. 
The equation of motion for the particles is Lf = Q(f) where Q(f) is an 
integral expression which is quadratic in its argument. This is the Boltzmann 
equation. It describes collisions between the particles in a statistical way. The 
case Q = is the collisionless case, where the equation Lf = obtained is 
often called the Vlasov equation. The energy-momentum tensor is defined by: 

T a0 = J fp a p du;(p) 

where duj(p) represents a natural measure on the cotangent space or mass 
shell, depending on the case being considered. 

The coupled Einstcin-Boltzmann system in harmonic coodinates cannot 
be a hyperbolic system is any usual sense for the simple reason that it is not 
even a system of differential equations, due to the integrals occurring. Never- 
theless, the techniques used to prove existence and uniqueness for hyperbolic 
equations can be adapted to prove local existence for the Einstein-Boltzmann 
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system |jj . In contrast to the case of the perfect fluid nothing particular hap- 
pens when the energy density vanishes so that there is no problem in de- 
scribing isolated concentrations of matter. The Einstein- Vlasov system can 
be used to describe globular clusters and galaxies. For the analogue of this 
in Newtonian theory see . Literature on the relativistic case can be found 
in©. 

Another kind of phcnomenological matter model is elasticity theory. Apart 
from its abstract interest, self-gravitating relativistic elasticity is of interest 
for describing the solid crust of neutron stars. As in the case of a perfect fluid 
the field equations are equivalent to the equation that the divergence of the 
energy-momentum tensor is zero. What is different is the nature of the matter 
variables and the way they enter into the definition of the energy-momentum 
tensor. This is complicated and will not be treated here. The Cauchy problem 
for the Einstein equations coupled to elasticity theory has been discussed by 
Choquet-Bruhat and Lamoureux-Brousse . A local existence theorem for 
the equations of relativistic elasticity has been proved by Pichon j73| for data 
belonging to Gevrey classes. These are classes of functions which are more 
special that C°° functions in that the growth of their Taylor coefficients is 
limited. However they do not have the property of analytic functions, that 
fixing the function on a small open set determines it everywhere. 

In non-relativistic physics, the Euler equations are an approximation to 
the Navier-Stokes equations where viscosity and heat conduction are ne- 
glected. The Navier-Stokes equations are dissipative with no limit to the 
speed at which effects can propagate. Mathematically this has the effect that 
the equations are parabolic with no finite domain of dependence. It is prob- 
lematic to find an analogue of the Navier-Stokes equations in general rel- 
ativity which takes account of the effects of diffusion and heat conduction. 
One possibility is to start from the Boltzmann equation, which does have a 
finite domain of dependence and try to do an expansion in the limit where 
the collision term is large. This is analogous to the Hilbert and Chapman- 
Enskog expansions in non-relativistic physics. The first attempts to do this led 
to equations which probably have no well-posed Cauchy problem (Landau- 
Lifschitz and Eckart models). Hiscock and Lindblom |5J| have shown that 
the linearization of these equations about an equilibrium state have solutions 
which grow at arbitrarily large exponential rates. In response to this other 
classes of models were developed where the fluid equations are symmetric 
hyperbolic. (For information on this see 0]). These models do have a well- 
posed Cauchy problem and the main difficulty seems to be to decide between 
the many possible models. Since the variables used to formulate the sym- 
metric hyperbolic system for the fluid are not differentiated in forming the 
energy-momentum tensor, the system obtained by coupling these fluids to 
the Einstein equations can be written in symmetric hyperbolic form. 

More general matter models can be obtained from those already men- 
tioned by combining different types of matter field. For instance there is the 
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charged scalar field which combines a scalar field and a Maxwell field. This 
often produces no extra difficulties at all for the local in time Cauchy problem 
due to the fact, mentioned above, that the system obtained by taking two 
symmetric hyperbolic systems together is also symmetric hyperbolic, pro- 
vided the coupling between the two systems is via terms of order zero. It 
is also routine to allow charged particles in a kinetic model, obtaining the 
Einstein-Maxwell-Boltzmann system [Q. Charged fluids are more compli- 
cated. The local in time Cauchy problem has been treated in two cases, those 
of zero conductivity and infinite conductivity. The latter model is also known 
as magnetohydrodynamics. It has been shown to be have a well-posed initial 
value problem only for data in Gevrey classes [[37J , J6^] . 



5.5 An Example of an 111- Posed Initial Value Problem 

It may be hard to appreciate the significance of a system of equations having 
a well-posed initial value problem since most examples which come up in 
practice do have this property. In this section we present an ill-posed example 
which is close to examples which relativists are familiar with. Consider the 
special relativistic Euler equations with equation of state p = kp where k < 0. 
A special solution is given by constant density and zero spatial velocity. Now 
consider the equations obtained by linearizing the Euler equations about this 
background solution. The unknowns in the linearized system will be denoted 
by adding a tilde to the corresponding unknowns in the nonlinear system. 
For simplicity we take the background density to be unity. The linearized 
equations are: 

d t p = -2(1 + k)d a u a 

It will be shown that given any T > there exist periodic initial data (po, Uq) 
of class C°° such that it is not true that there is a unique corresponding 
solution, periodic in the space coordinates, on the time interval [0, T). The 
condition of periodicity here does not play an essential role. It is adopted for 
convenience. Instead of thinking of smooth periodic functions on R 3 with can 
equally well think in terms of smooth functions on a torus T 3 . The space of 
smooth functions on T 3 can be made into a topological vector space X in a 
standard way. Convergence of functions in the sense of this topology means 
uniform convergence of the functions and their derivatives of all orders. In a 
similar way the space of smooth functions on [0, T] x T 3 can be made into 
a topological vector space Y. These are Frechet spaces The unknown in 
the linearized Euler equations can be thought of as an element of X 4 and 
the data as an element of Y 4 . Let Z be the closed linear subspace of X 4 
consisting of solutions of the linearized Euler equations. Consider the linear 
mapping L : Z — ► Y 4 defined by restricting solutions to t = 0. It is continuous 
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with respect to the relevant topologies. If existence and uniqueness held for 
all initial data then this linear map would be invertible. Hence, by the open 
mapping theorem it would have continuous inverse. What this means 
concretely is that given any sequence (po,m Uq n ) which is uniformly bounded 
together with each of its derivatives, the corresponding sequence of solutions 
(which exists by assumption) must also be bounded together with each of each 
of its derivatives. Thus to prove the desired theorem it is enough to exhibit 
a uniformly bounded sequence of initial data and a corresponding sequence 
of solutions which is not uniformly bounded. This can be done explicitly as 
follows: 

Po. n = sinnx 

<,n = 

p n = sin nx cosh(nV—kt) 

_-, \J—kn . , — - 

u„ = — r sin nx Sinn inv — kt) 

2(1 + jfe) v ' 

«n = fi n = 

From these explicit formulae we get an idea what is going wrong. Fourier 
modes of increasing frequencies of the initial data grow at increasing expo- 
nential rates. In the case of a fluid where the equation of state has a positive 
value of k the hyperbolic functions in the above formulae are replaced by 
trigonometric ones and the problem does not arise. 

In fact in the above example the density perturbation satisfies a second 
order equation which is elliptic. After a rescaling of the time coordinate it 
reduces to the Laplace equation. The computation which has just been done 
should be compared with the remarks on the Cauchy problem for the Laplace 
equation on p. 229 of ^7j, Vol. 2. The corresponding example for the Laplace 
equation goes back to Hadamard [ff9| . 

Solutions of the Einstein equations coupled to a fluid with an equation of 
state of the type considered in this section have been considered in the context 
of inflationary models J93|. While this is unproblematic for spatially 

homogeneous models, the above ill-posedness result suggests strongly that 
this kind of model cannot give reasonable results in the inhomogeneous case. 

Some general comments on well-posedness and stability will now be made. 
Suppose a solution of a system of evolution equations is given. Assume for 
simplicity that this solution is time- independent, although a similar discus- 
sion could be carried out more generally. The solution is called stable if in 
order to ensure that a solution stays close to the original solution to any 
desired accuracy, it is enough to require it to be sufficiently close at one time. 
Closeness is measured in some appropriate norm. Well-posedness has no in- 
fluence on stability in this sense. Already for ordinary differential equations 
with smooth coefficients, which always have a well-posed initial value prob- 
lem, stability does not in general hold. Solutions of the linearized sytem about 
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the given solution can grow exponentially. However there is a constant k such 
that no linearized solution can grow faster than Ce kt . This is a rather general 
feature of well-posed evolution equations. For instance any linear symmetric 
hyperbolic system allows an exponential bound with some constant k inde- 
pendent of the solution. If no such bound is possible, then it is said that there 
is a violent instability. (Cf. (69j, Sect. 4.4 for this terminology.) The presence 
of a violent instability is closely related to ill-posedness, as can be seen in the 
above example. 

5.6 Symmetries 

If an initial data set for the Einstein-matter equations possesses symmetries, 
then we can expect these to be inherited by the corresponding solutions. This 
will be discussed here in the case of the vacuum Einstein equations. There is 
nothing in the argument which obviously makes essential use of the vacuum 
condition and it should extend to reasonable types of matter. It makes use 
of the maximal Cauchy development and so any restrictions on the matter 
model which might come up there would appear again in the present context. 

The following only covers symmetries of spacetime which leave a given 
Cauchy surface invariant. It is based on group actions rather than Killing 
vectors. A more extensive discussion of symmetries of spacetime and their 
relations to the Cauchy problem can be found in Sect. 2.1 of p6fl . 

By a symmetry of an initial data set (S, /i 6, Xab) for the vacuum Einstein 
equations we mean a diffcomorphism ip : S — ► S which leaves h a b and Xab 
invariant. Let <fi be the embedding of S into its maximal Cauchy development. 
Then <p — <poip also satisfies the properties of the embedding in the definition 
of the maximal Cauchy devlopment. Hence, by uniqueness up to isometry, 
there exists an isometry ip of the maximal Cauchy development onto itself 
such that ip o <p = (f>. This means that ip o <p = cp o ip. Thus ip is an isometry 
of M whose restriction to 4>(S) is equal to ip. We see that a symmetry of 
the initial data extends to a symmetry of the solution. Next we wish to 
show that this extension is unique. Since a general theorem of Lorentzian 
(or Riemannian) geometry says that two isometries which agree on a open 
set agree everywhere it suffices to show that any two isometries ip with the 
properties described agree on a neighbourhood of ^(S 1 ). Let p be a point of 
4>(S). A neighbourhood of p can be covered with Gauss coordinates based 
on <p(S). An isometry preserves geodesies and orthogonality. Hence if, when 
expressed in Gauss coordinates, it is the identity for t = it must be the 
identity everywhere. This completes the proof of the uniqueness of ip. 

Now consider the situation where a Lie group G acts on S in such a 
way that each transformation tp g of S corresponding to an element of the 
group is a symmetry of the initial data. Let H be the isometry group of 
the maximal Cauchy development and H$ the group of all isometries of the 
maximal Cauchy development which leave (p(S) invariant. The group H$ is 
a closed subgroup of the Lie group H and thus is itself a Lie group. Each 



The Cauchy Problem 91 



tpg is the restriction of a unique element ip g of H$- Using uniqueness again 
we must have ip g h = ipgiph for all elements g and h of G. Thus we obtain a 
homomorphism from G to Hs- This shows that there exists an action of the 
group G on M by isometries which extends the action on 4>{S) arising from 
the original action on S by means of the identification using <fi. However this 
argument does not show that the resulting action of G is smooth. To show 
this consider first the group Hi of all symmetries of the initial data. It is a 
closed subgroup of the isometry group of h and therefore has the structure 
of a Lie group. The above considerations show that restriction defines an 
isomorphism of groups from Hs to Hj. (We identify S with 4>(S) here.) If 
we knew that this mapping was continuous a general theorem on Lie groups 
[ j92| would show that it is also an isomorphism of Lie groups. The continuity 
can be seen by noting that the topology of an isometry group coming from 
its Lie group structure coincides with the compact open topology ||62)| . The 
continuity of the restriction mapping in the compact open topology follows 
immediately from the definitions. We conclude that as Lie groups Hs and 
Hi can be identified. 

Now we come back to the action of G. The action of G on initial data 
is a smooth mapping G x S — > S. It is the composition of a smooth homo- 
morphism from G to Hi with the action of Hi on S. By the comments of 
the last paragraph this can be identified with the composition of a smooth 
homomorphism from G to Hs with the action of Hs- In this way we obtain 
a smooth action of G on M which leaves S invariant and restricts to the 
original action on the initial data. It is the action of G which we previously 
considered. 

6 Outlook 

This article is intended to be an informative tour through its subject, rather 
than an exhaustive account. The latter would in any case be impossible in 
an article of this length, given the amount of literature which now exists. 
The aim of this section is to mention a few of the important things which 
have been left out, and to direct the reader to useful sources of information 
concerning these. A good starting point is the review article of Choquct- 
Bruhat and York |^3| which is still very useful. (See also There is an 

extensive treatment of the constraints in reference ]23| . For a selection of 
newer results on the constraints, see e.g. [||, ||, p7{, |5£{], p9fl and ]2C| . 

The most obvious omission of the present article is the lack of statements 
on the global Cauchy problem. A review with pointers to further sources can 
. This material is too recent to be discussed in Q and a 
lot has happened since then. It is natural that once some of the basic local 
questions had been solved attention turned to global issues. The latter are 
now central to present research on the Cauchy problem. Most of the existing 
results concern spacetimes with high symmetry, although in the meantime 



be found in |78 
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there are also a few theorems on space-times without symmetries which are 
small but finite perturbations of space-times of special types. This is an area 
which is developing vigorously at the moment. Often the statements obtained 
about global existence of solutions are accompanied by information on the 
global qualitative properties of the solutions. A discussion of the asymptotic 
behaviour of a particular class of solutions and its relevance for the modelling 
of certain physical systems can be found in |l0| and JOJ . 

In the introduction we mentioned the possible applications of ideas con- 
nected with the Cauchy problem to analytical and numerical approximations. 
Up to now progress on establishing an effective interaction between theoret- 
ical developments and the applications of approximate methods to concrete 
physical problems such as the generation of gravitational waves has been lim- 
ited. For instance, little has been done on the question of proving theorems 



on analytical approximations in general relativity since E8|, 75 , and fl? 



There is no shortage of things to be done. For instance one tempting goal 
would be a precise formulation and justification of the quadrupole formula. 

As for the link to numerical relativity, the discussions in Sect. ||, apart 
from their interest for purely analytical reasons, could potentially be ex- 
ploited for improving numerical codes. Many new hyperbolic reductions have 
been suggested recently with the aim of providing equations which would 
ensure a stable time evolution. We have added a few more. We also pointed 
out various different gauge conditions. Their usefulness for stable long-time 
numerical calculations still has to be explored. We have seen that different 
representations of the field equations and different formalisms allow us to em- 
ploy different gauge conditions. The possibilities of controlling the lifetime of 
a gauge by a judicious choice of gauge source functions have neither been 
investigated analytically nor numerically in a systematic way. 

One of the main interests in the analysis of the initial boundary problem 
lies in the fact that many approaches to numerical relativity require the in- 
troduction of timclike boundary hypersurfaces which reduce the calculations 
to spatially finite grids. A good analytical understanding of the initial bound- 
ary value problem does not guarantee the stability of long-time evolutions for 
problems with timelike boundaries, but the latter are likely to fail without 
a proper understanding of the analytical background. In this respect we also 
consider the analytical investigations of the conformal Einstein equations (cf. 
p0| , p[ ) as an opportunity for numerical relativity. They allow us in princi- 
ple to calculate entire spacetimes, including their asymptotic behaviour, on 
finite grids without the need to introduce artificial boundaries in the physical 
spacetime. 

Another place where the theory of the Cauchy problem could have some- 
thing to contribute is the application of numerical methods developed for 
systems of conservation laws (or, more generally, systems of balance laws) 
to the Einstein equations, a procedure which has been popular recently. In a 
certain sense this means accepting an analogy between the Einstein equations 
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and the Euler equations. At the moment this amounts to no more than the 
fact that both systems are quasi-linear hyperbolic and can be formulated as 
systems of balance laws. Unlike the Euler equations, the Einstein equations 
do not appear to admit a preferred formulation of this kind. There are differ- 
ent alternatives and no known criterion for choosing between them. It is very 
tempting to import the vast amount of knowledge which has been accumu- 
lated concerning the numerical solution of the Euler equations into general 
relativity. On the other hand it is not at all clear how many of these tech- 
niques are really advantageous for, say, the vacuum Einstein equations. One 
might hope that analytical theory could throw some light on these questions. 

There can be little doubt that increased cooperation between people work- 
ing on analytical and numerical aspects of the evolution of solutions of the 
Einstein equations would lead to many new insights. As both fields progress 
cases where a rewarding collaboration would be possible are bound to present 
themselves. It suffices for someone to show the initiative required to profit 
from this situation. 

An open problem which has been touched on already is that of the exis- 
tence of solutions of the Einstein-Euler system describing fluid bodies. This 
kind of free boundary problem is poorly understood even in classical physics, 



although there has been significant progress recently 1 95 . The central impor- 
tance of this becomes clear when it is borne in mind that the usual applica- 
tions of gravitational theory in astrophysics (except for cosmology) concern 
self-gravitating bodies. 

We hope to have succeeded in showing in this article that the study of the 
Cauchy problem for the Einstein equations is a field which presents a variety 
of fascinating challenges. Perhaps, with luck, it will stimulate others to help 
tackle them. 
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